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ABSTRACT 

An  analysis  was  done  of  central  Arctic  Ocean  acoustic  data, 
to  determine  the  temporal  and  spatial  characteristics  o-f 
transient  noise  events.   Digital  ambient  noise  data  from 
the  FRAM  IV  experiment  of  April  1982  were  searched  for 
ambient  noise  transients  using  a  detection  program.   The 
time  series  of  the  resulting  detections  were  examined 
visually  to  categorize  each  detection  as  a  transient, 
artifact  or  false  alarm.   The  transient  events  were  located 
in  space  using  time  delays  between  signal  arrival  at 
different  hydrophones.   The  cross  shape  of  the  FRAM  IV 
horizontal  array  permitted  location  in  both  bearing  and 
range.   The  source  strength  of  each  event  was  calculated 
using  a  simole  dipole  source  model.   Refraction  and 
scatterring  of  the  acoustic  path  in  the  Arctic  Ocean  was 
taken  into  account. 

The  overall  number  of  events  detected,  and  hence  their 
interarrival  times  and  spatial  density,  v^ere  all  affected 
by  the  background  ambient  noise  level.   The  detection 
program  used  the  same  threshold  si  gnal -tc— noi  se  le-/el  for 
all  data  tapes,  so  when  ambient  noise  levels  were  low  more 
detections  occurred.   The  mean  interarrival  time  between 
events  was  100  seconds.   The  interarrival  time  fit  a  J 
shaped  gamma  probability  distribution.   The  number  of 
events  detected  per  a.rsa    decreased  with  range  from  the 
Arr 3.y    center.   Half  of  the  events  occurred  within  3000 
meters  of  the  array.   In  this  area  there  were  0.3  events 
per  square  kilometer  per  hour.   The  event  population  showed 
no  predcminant  angular  dependence.   The  strengths 
calculated  using  the  simple  dipole  model  had  a  mean  o-r  430 
kN  overall  and  260  kN  during  quiet  times.   Stronger  events 
occur r^d  during  times  '-^ith  high  ambient  ncise  levels. 

Theis  Supervisor:   Dr.  Ira  Dyer 

Title:   Professor  of  Ocean  Engineering 


TABLE  OF  CONTENTS 

Abstract  2 

Table  o-f  Contents  3 

List  D-f  Figures  4 

List  o-f  Tables  7 

Chapter  1  -  Introduction  S 

Thesis  Motivation  S 

Thesis  Contents  12 

Chapter  2  —  Detection  o-f  Noise  Events  13 

Data  Collection  13 

Event  Detector  Program  16 

Visual  Con-f  irmation  2S 

Chapter  3  —  Location  o-f  Noise  Events  35 

Event  Location  Program  35 

E-f-fects  o-f  Refraction  on  Location  47 

Chapter  4  —  Strength  o-f  Noise  Events  56 

Acoustic  Source  Model  56 

E-f-fects  G-?  Re-fraction  on  Transmission  Lc<b5  5'^' 

Strength  o-f  Background  Noise  71 

Chapter  5  -  Analysis  at  Noise  Evants  76 

Detection  Analysis  76 

Temporal  Analysis  79 

Spatial  Analysis  91 

Strength  Analysis  9S 

Chapter  6  —  Summary  and  Thoughts  108 

Re-ferences  110 

Ackrow 1 edgements  112 

HDcendix  A  —  Event  Detector  Program  il4 

Appendix  B  -  Location  Program  133 

Appendix  C  -  Event  Summary  144 

ADoendix  D  -  Re-*^ractian  Path  Calculations  164 


-4- 


LIST  OF  FIGURES 


Figure  1-1:   Location  o^^  the  FRAM  IV  Arctic  11 

experiment  conducted  in  the  spring  o-f  1982. 

Figure  2-1:   FRAM  IV  horizontal  hydrophone  array.        14 

Figure  2-2:   Schematic  of  the  recording  system  15 

used  -for  FRAM  IV  data  collection. 

Figure  2—3:   Block  diagram  o-f  major  modules  o-f  the      18 
detect  ioTi    program. 

Figure  2—4:   Frequency  response  o-f  the  Parks—  19 

McClellan  digital  bandpass  filter  used  in 
detect  ion . 

Figure  2—5:   Diagram  of  event  detection  module  21 

flow  and  decision  making. 

Figure  2—6:   Sketch  showing  point  of  measurement        33 

for  event  time  delays.   The  event  is  timed  at  its 

zero  crossing  between  the  largest  pair  of  positive 

and  negative  peaks.    (The  measurement/analysis 

system  has  a  polarity  of  negative  voltage  for 

positive  pressure.) 

Figure  3—1:   Coordinate  system  used  for  36 

calculation  of  event  location. 

Figure  3—2:   Least  squares  fit  of  time  delays  and       37 
test  location  slant  range. 

Figure  3—3:   Least  squares  fit  to  a  false  location      41 
180°  away. 

Figure  3—4:   Areas  covered  by  locatioTt  ,    -far  locate  42 

and  tine  locate . 

Figure  3—5:   Noise  events  located  within  a  2  km         45 
square  surrounding  the  array  origin. 

Figure  3—6:   Position  of  all  noise  events  located.       46 

Figure  3-7:   Sound  velocity  profile  used  in  48 

predicting  refractive  paths. 


Figure  3-8:   Scheme  used  to  compute  re-fractive  49 

propagation  paths.   Source  is  located  at  z  =  Om 
and  hydrophone  at  z  =  93m.   Sound  speed  gradients 
change  at  30m,  254m  and  362m. 

Figure  3-9:   Time  di-f-ference  between  the  slant  52 

range  acoustic  path  and  the  re-fractive  path  as  a 
function  o-F  horizontal  range. 

Figure  3-10:  Scatter  plot  o-f  standard  deviation         53 
and  range  -for  all  located  noise  events- 
Figure  4—1:   Assumed  dipole  source  orientation  and      57 
de-finition  o-f  launch  angle,   Q  . 

Figure  4-2:   Spherical  spreading.  60 

Figure  4-3:   Spreading  loss  in  re-fraction.  62 

Figure  4—4:   Specular  and  non— specular  re-flections      65 
■from  the  ice/water  inter-face. 

Figure  4—5:   Probability  of  a  ray  bundle  crossing       67 
a  certain  depth  at  a  given  horizontal  range. 

Figure  4—6:   Spherical  (heavy  line)  and  refractive      70 

spreading  loss,  G(r) ,  as  a  function  of  horizontal 

range. 

Figure  4—7:   Composite  central  Arctic  ambient  72 

noise  spectrum  observed  during  the  FRAM  IV 
ex per i  ment . 

Figure  5—1:   Normalized  ambient  noise  pressure,  77 

number  of  unique  events  per    tape,  and  number  of 
false  alarms  per  tape  for  each  data  tape  examined. 

Figure  5—2:   Average  number  of  false  alarms  and         78 
unique  events  per  data  tape  for  four  ranges  of 
ambient  noise  pressure. 

Figure  5—3:   Number  of  events  per  interarrival  SO 

t  i  me  bin. 

Figure  5-4:   Hal f -gaussi an  distribution  compared        S4 
with  experimental  values. 

Figure  5-5:   Exponential  distribution  compared  B5 

with  experimental  values. 


-6- 


Figure  5-6:   Semi  log  plot  of  exponential  86 

distribution  and  experimental  values  against  bin 
number. 

Figure  5-7:   J  shaped  distribution  compared  with        88 
experimental  values. 

Figure  5-Q:   Semi  log  plot  o-f  J  shaped  89 

distribution  and  experimental  values  against  bin 
number . 

Figure  5-9:   Scheme  for  annul i  each  of  area  equal       92 
to  30  km'^. 

Figure  5—10:  Number  of  events  per  radius  annulus.        93 

Figure  5—11:  Number  of  events  per  30°  sector.  96 

Angles  Are    measured  from  the  northern  leg  of  the 
hydrophone  array. 

Figure  5-12:  Number  of  events  per  30°  sector.  97 

Radius  gives  the  number  of  events,  while  angle 
indicates  the  sector  measured  from  the  northern 
leg  of  the  array.   Each  ring  represents  10  events. 

Figure  5—13:  Number  of  events  per  10°  sector.  99 

Radius  gives  the  number  of  events,  while  angle 
indicates  the  sector  measured  from  the  northern 
leg  of  the  array. 

Figure  5—14:  Mean  hydrophone  peak  pressure  100 

measured  for  events  between  100  m  and  20,000  m, 
plotted  against  horizontal  range  from  the  FRAM  IV 
arraty    origin. 

Figure  5—15:  Distribution  of  strength  for  a  103 

population  of  noise  events. 

Figure  5—16:  Dipole  strength  for  events  between        104 
300  m  and  20,000  m,  plotted  against  horizontal 
range  from  the  Array    origin. 

Figure  5-17:  10  log  of  Dipole  strength  (dB  re  1  N)      106 
versus  10  log  of  horizontal  range  from  the  array 
origin  (dB  re  1  m)  for  events  occurring  during  an 
ambient  noise  level  of  0.01  to  0.02  Pa. 

Figure  A— 1:   Flow  chart  of  the  hdetect    program.         119 


LIST  OF  TABLES 


Table  2-1:   Determination  of  the  Best  RATIO  25 

Table  2-2:   Detection  Statistics  for  Two  Visual  31 

Con-f  i  rmati  on  Methods 

Table  3-1:   Results  o-f  locat ion    Program  Accuracy        55 
Test 

Table  4-1:   20  to  80  Hz  Band  Ambient  Noise  rms  75 

Pressure 

Table  3—1:   Number  o-r  Events  per  Inter  am  val  Time      31 
Bin 

Table  5—2:   Comparison  o-f  Distribution  Functions        90 

Table  5-3:   Background  Noise  Level  Dependence  of        91 
Interarrival  Time 

Table  5-4:   Number  o-f  Events  per  Annul  us  94 

Table  5—5:   Average  Number  o-f  Events  per  Annul  us         95 
-for  4  Amoient  Noise  Levels 


Table  5—6:   Mean  Hydrophone  Peak  Pressure  for  4 
Ambient  Noise  Levels 


1  n  1 


Table  5-"?':  Strength  distribution  -for  a  Fcoulatian  lOv 
of  Events 

Table  5-B:  Dioole  Strength  versus  Ambient  Noise  102 
Levels 

Table  C— 1:  Event  Location  Summary  145 

Table  C— 2:  Tape  Summary  152 

Table  C-3:  Event  Interarrival  Time  Summary  153 

Table  C-4:  Event  Strength  Summary  159 

Table  D— 1;  Angles,  Ranges  ana  Times  for  165 

Ref  r a c t  i  ve  Prop  ag a. t  i  on  Paths 

Table  D— 2:  Soreadinq  Loss  Function,  G(r)  IGO 


-a- 

CHAPTER  1 
INTRODUCTION 

Thesis  Motivation 

This  thesis  investigates  the  spatial  distribution, 
strength  and  rate  o-f  occurrence  o-F  low  frequency  central 
Arctic  noise  events.   During  the  last  30  years  there  has 
been  a  growing  commercial,  military  and  academic  interest 
in  the  Arctic  region. 

In  this  relatively  unexplored  area  there  is  increasing 
evidence  of  rich  mineral  and  petroleum  resources.   Research 
into  methods  of  locating  these  assets  and  constructing 
facilities  to  exploit  them  have  received  much  attention. 
These  facilities  must  be  able  to  withstand  the  harsh  Arctic 
surroundings.   The  study  of  Arctic  acoustics  helps  in 
understanding  the  Arctic  environment  and  climate.   It  has 
been  shown  that  there  is  a  direct  correlation  between  iO  to 
20  Hs  ambient  sound  pressure  and  environmental  stresses  and 
sTiomentsC  103 .   The  ability  to  use  acoustic  noise  levels  as 
an  environmental  predictor  would  be  a  useful  tool  in  the 
protection  of  commercial  Arctic  facilities. 

The  Arctic  ocean  serves  as  a  military  arena  for 
several  submarine  fleets.   The  under ice  environment  makes 
detection  difficult,  increasing  the  stategic  role  of  those 
fleets.   Because  of  the  sound  velocity  profile  of  the 
central  Arctic  there  is  a  surface  duct  which  channels  sound 
for  long  distances.   But,  the  underice  profile  scatters 


-9- 

BOLind  energy,  e-f -f  ecti  vel  y  filtering  out  high 
■FrequenciesC  1  ID  .   The  result  is  that  only  low  frequency 
signals  travel  far  in  the  Arctic,  and  therefore  the  low 
frequency  range  is  the  best  for  detecting  adversary 
submarines.   The  importance  of  understanding  the  low 
frequency  ambient  noise  field  becomes  apparent.   The  actual 
central  Arctic  ambient  noise  level  is  at  times  much  quieter 
than  the  open  ccean,  but  it  contains  unpredictable 
transient  noise  events  which  interfere  with  conventional 
detection  schemes.   It  has  been  hypothesized  that  the 
background  ambient  noise  is  the  summation  of  these 
transients  from  throughout  the  Arctic  basinC7Il.   Analysis 
of  the  spatial  and  temporal  distribution  of  these 
transients  is  a  logical  next  step  in  understanding  low 
frequency  noise,  and  improving  our  submarine  detection 
capabi 1 i  ty . 

The  academic  challenge  of  the  Arctic  lies  in  the 
sparseness  of  field  data.   The  Arctic  cannot  be  casual iy 
sampled.   Even  simple  experiments  require  expensive 
expeditions.   The  harsh  environment  takes  its  toll  on 
researchers  and  equipment,  and  reduces  the  amount  of  usable 
data.   Hence,  the  study  of  the  Arctic  is  like  a  jigsaw 
puzzle  with  few  pieces  present.   The  total  picture  remains 
a  stimulating  mystery. 

This  study  analyzes  data  collected  during  the  FRAM  IV 
experiment  by  Massachusetts  Institute  of  Technology  and 
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Woods  Hole  Oceanographic  Institute  personnel.   The  FRAM  IV 
ice  camp  was  located  in  the  Barents  (Nansen)  Abyssal  Plain 
at  approKimately  34°  M  by  15°  E,  as  shown  in  Figure  1-1. 
The  ice  was  3  meter  thick  multi-year  pack  ice.   The  ice 
activity  was  low;  there  was  no  ice  ridging  or  lead 
formation  around  the  camp  during  the  experiment. 

The  FRAM  IV  ice  camp  was  set  up  from  25  March  to  11 
May  1982.   This  study  analyzes  data  taken  bet>jeen  March 
27th  and  April  22nd.   The  weather  was  mild,  with 
temperature  ranged  from  -35°  to  —4°  C,  and  wind  speed  from 
1  to  23  knots. 

The  ambient  noise  was  sensed  with  a  large  horizontal 
hydrophone  array  which  consisted  of  two  non— uni f ormly 
spaced  line  arrays,  crossing  at  right  angles.   The  data 
were  digitally  recorded  on  a  multichannel  system. 
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Figure  1-1   Location  o-f  the  FRAM  IV  Arctic  experiment 
conducted  in  the  spring  o-f  1982.  llil 
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Thesis  Contents 

My  work  in  this  airsA    began  long  after  the  FRAM  IV 
ambient  noise  data  were  collected.   The  -first  step  was 
■finding  the  events  in  the  raw  data.   The  data  consisted  of 
magnetic  tapes  each  containing  20  minutes  of  digitized 
noise  levels.   A  program  was  written  which  searched  the 
ambient  noise  tapes  for  possible  events.   Chapter  2 
di scribes  the  automated  and  manual  techniques  used  to 
accomplish  this  detection. 

These  events  were  then  located  in  space  using  the 
difference  in  arrival  time  between  hydrophones.   This  was 
also  done  with  a  computer  program.   The  program  plotted  the 
arrival  time  delays  against  range  to  a  trial  location,  did 
a  least  squares  fit,  and  chose  the  location  with  the  best 
fit.   This  is  covered  in  Chapter  3. 

The  peak  voltages  for  each  event  were  used  with 
the  dipole  source  model  to  predict  peak  source  strength. 
The  background  ambient  noise  strength  was  also  determined. 
These  strength  calculations  aire    found  in  Chapter  4. 

In  Chapter  5  the  distribution  of  event  interarrival 
time  was  determined.   The  event  locations  and  strengths 
were  analyzed,  and  a  spatial  density  found. 

Chapter  6  summarizes  the  key  results  of  this  study. 
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CHAPTER  2 
DETECTION  OF  NOISE  EVENTS 

Data  Collection 

Twenty  nine  FRAM  IV  ambient  noise  tapes  were  searched 
in  order  to  find  a  population  o-f  noise  events  -for  this 
study.   The  specific  tapes  were  chosen  -from  the  possible  67 
in  order  to  cover  the  entire  range  o-f  days  of  the  FRAM  IV 
experiment.   However,  there  were  several  days  when  no 
ambient  noise  tapes  were  recorded.   To  help  fill  these  gaps 
five  reverberation  tapes  were  also  searched.   These  tapes 
were  recorded  prior  to  the  reverberation  shot  being  fired, 
ar    they  were  recorded  so  late  in  the  experiment  (SO  minutes 
after  the  shot)  that  reverberations  were  no  longer  present.- 

The  FRAM  IV  experiemnt  used  the  horizontal  array  of 
omnidirectional  hydrophones  pictured  in  Figure  2-1.   The 
hydrophones  wer=  suspended  from  the  ice  into  the  water  to  a 
depth  of  93  meters  below  the  air/ice  interface.   The  two 
crossed  lines  of  the  array  allowed  the  possibility  cf 
localising  events  in  space.   Although  26  hydrophones  are 
shown  in  Figure  2-1,  only  24  at  a  time  could  be  used  to 
record  data.   In  most  cases  a  few  of  the  recording  channels 
were  used  for  other  sensors  (geophones  or  hydrophones  used 
in  a  vertical  array).   Most  of  the  time  19  to  21  horizontal 
array  h ydrcuhone  data  were  recorded. 

The  FRAM  IV  ambient  noise  tapes  were  recorded 
digitally.   Figure  2-2  shows  a  schematic  of  the  system  used 
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FRAM  IV  DATA  ACQUISITION 
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-for  collecting  the  data.   The  input  -from  each  hydrophone 
was  put  through  a  gain  ranging  amplifier  and  a  low  pass  80 
Hz  filter.   It  was  sampled  at  250  Hz  and  recorded  on  20 
minute  magnetic  tapeC4;].   The  24  channel  recorder  had  a  120 
dB  dynamic  range[!12D. 

Event  Detector  Program 

The  event  detector  program  was  written  to  take  digitel 
data  -from  a  FRAM  IV  tape  and  determine  where  in  that  tape 
noise  events  occurred.   The  program  was  originally  written 
to  take  data  directly  -from  a  tape  drive,  but  subsequently 
modified  to  take  the  data  from  a  file.   The  t'ramread 
program,  with  a  -head  switch  is  used  to  read  the  tape  into 
the  file.   This  will  eliminate  any  headers  and  then  read 
the  digital  data  straight  into  a  file.   A  FRAM  II  tape  may 
be  read  into  a  similar  file  using  t'ramread    and  the  switches 
—head  and  — fram2.   The  tramread    program  was  written  by 
G.  Duckworth,  and  is  available  to  .the  Arctic  Acoustics 
Program  at  MIT. 

The  event  detection  programs  source  codes,  flow  chart, 
and  a  short  users  guide  are    found  in  Appendix  A.   The  event 
detection  program  which  reads  from  a  file  is  called 
hdetect .        The  detect  ion    program  reads  from  a  tape  drive. 
Both  prograrr^B  are  written  in  the  c  programming  language  for 
a  UNIX  operating  system. 

The  event  detection  program  follows  the  block  diagram 
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o-f  Figure  2-3.   The  program  i  ni  ti  al  i  zati  an  portion  defines 
variables  and  constants,  zeros  -flags,  and  requests  user 
input  such  as  tape  number,  date,  time  and 

channels (hydrophones)  to  be  used,  as  well  as,  the  name  ai 
input  and  output  -files.   After  this  information  is 
requested  from  the  user,  the  program  no  longer  requires 
attenti  on. 

The  event  detection  program  then  reads  in  a  file  of 
data,  filters  the  data,  squares  each  data  point,  and  takes 
the  square  root.   The  filter  was  a  Parks-McCl el  1  an  digital 
20  to  SO  Hz  bandpass  filter.   Its  frequency  response  is 
shown  in  Figure  2—4.   The  range  of  this  filter  was  chosen 
to  avoid  the  Nyquist  frequency  (125  Hz)  and  hydrophone 
cable  strum  (1-20  Hz) ,  and  to  be  compatible  with  the  analog 
30  Hz  low  pass  filter  the  data  went  through  before  being 
recorded.   The  data  -vjere  squared  and  then  rooted  to  ensure 
positive  peak  values  for  all  data  points. 

The  -ie>;t  portion  of  the  program  used  a  threshold 
detection  scheme  to  check  each  channel  for  possible  noise 
events.   For  a  particular  channel  a  short  average  of  the 
four  most  recent  data  points  was  compared  to  a  long  average 
of  64  recent  data  points.   If  the  ratio  of  short  average  to 
long  average  was  over  a  certain  value  that  channel  would  be 
flagged  -for  a  possible  event.   The  time  of  the  flag  and  the 
value  of  the  short  average  were  also  recorded.   All  other 
channels  averaoes  were  taken  similarly.   If  at  least  507.  of 
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Figure  2—3   Block  diagram  o-f  major  modulues  o-f  the 
detect  ion    program. 
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Figure  2—4   Frequency  response  o-f  the  Parks-McCl  el  1  an 
digital  bandpass  -filter  used  in  detect  ion . 
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the  channels  were  flagged,  an  event  would  be  declared. 
Then  the  program  would  shift  to  the  next  time  increment  of 
data  and  the  process  would  be  repeated. 

A  more  detailed  diagram  of  the  event  detection  module 
is  seen  in  Figure  2-5.   There  are  four  submodules:  reset 
flag,  set  flag,  new  event,  and  deactivate  old  event 
modules.   The  reset  and  set  flag  modules  deal  with  the 
channel  flags  which  trip  when  a  particular  channel 
experiences  a  large  si gnal — to— noi se  ratio  (i.e.  the  RATIO 
of  short  average  to  long  average  exceeds  a  certain  level). 
The  new  event  and  deactivate  old  event  modules  deal  with  an 
active  event  matrix  which  identifies  active  events,  and 
stores  channel  flag  time  and  amplitude  for  each  declared 
event . 

The  reset  flag  module  resets  the  channel  flag  if  it 
has  been  more  than  0.3  seconds  since  the  channel  tripped. 
Spurious  peaks  on  a  channel  might  flag  a  channel 
prematurely.   This  reset  module  prevents  a  number  of 
channels  with  spurious  peaks  over  a  long  time  period  from 
being  falsely  declared  an  event.   The  value  of  this  RESET 
DELAY  was  determined  by  examining  known  events  and  noting 
that  about  half  the  channels  tripped  within  a  0.3  second 
peri  ad . 

The  set  f 1 aq  module  determines  if  the  short  average  to 
long  average  RATIO  has  been  exceeded  and,  if  it  has,  the 
module  1)  checks  to  see  if  the  channel  flag  is  already 
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Figure    2—5       Diagram    o-f    event    detection    module    -flow    and 
decision    making. 


triDped,  2)  checks  to  see  i -f  this  detection  is  part  o-f  a 
recently  declared  event,  and  3)  adds  its  detection 
i  n-f  ormat  1  on  to  the  active  event  or  the  channel  + 1  ag  ,  or 
discards  the  in-formation  depending  on  the  circumstances. 
The  short  average  length  g+  4  <=  0.016  sec)  was  chosen  so 
that  it  matched  the  length  o-f  the  signal  (0.02  sec).   This 
provides  the  maximum  signal  level  since  this  is  long  enough 
to  get  all  o-f  the  signal  and  short  enough  not  ~o  a /erage  it 
f.'^ith  the  lower  surrounding  background  noise  level.   The 
long  average  length  o-f  64  was  chosen  because  a  length  ratio 
o-f  15  to  1  had  been  suggested  by  KellyC9Il  -for  the  L.arge 
Aperture  Seizmic  Array  (LABA).   The  detection  and 
localization  schemes  used  by  this  large  horizontal  e^rrs^y 
i-vere  directly  applicable  to  the  FRAM  IV  hydrophona  array 
data.   The  choice  o-f  detection  RATIO  was  done  through  a 
series  o-f  tesrs.  and  the  selection  was  made  by  balancing 
detection  rate  and  -false  alarm  rate. 

when  an  event  is  declared  the  in-formation  in  the 
channel  -flags  is  trans-f erred  to  the  active  event  matrix, 
and  the  channel  -flags  a.rs    cleared.   So  if  a  channel  has  a 
detection  and  its  channel  -flag  is  not  alreadv  tripped,  the 
set  flag  module  must  first  see  if  the  detection  belongs  to 
a  recently  declared  event.   If  the  active  event  alreadv  has 
that  channel  flagged,  the  information  is  replaced  bv  the 
new  detection  only  if  their  times  differ  by  less  than  0.02 
seconds  and  the  new  detection  amplitude  is  or eater.   This 


means  that  the  existing  data  can  only  be  replaced  by  a 
detection  of  the  same  signal  having  a  higher  peak  value. 
I-f  the  active  event  does  not  have  that  channel  flagged,  the 
detection  information  is  entered  in  the  active  event 
matrix,  and  the  channel  flag  is  not  tripped.   The  oldest 
active  events  atre    checked  first,  and  the  detection 
information  entered  in  only  one  event.   If  all  of  the 
active  events  already  have  input  for  tnis  channel,  anc  it 
has  been  more  than  0.02  sec  since  the  most  recent  input, 
the  detection  is  considered  potentially  part  of  an 
undeclared  event,  and  its  channel  flag  is  tripped. 

If  a  detection  is  made  on  a  channel  and  its  channel 
flag  IS  already  tripped,  the  information  in  the  channel 
flag  will  be  replaced  with  the  new  detection  only  if  the 
new  detection  amplitude  is  greater. 

The  new  event  module  checks  to  see  if  at  least  50"';  of 
the  channel  flags  have  been  tripped,  and  if  so,  declares  a 
new  event.   The  50'/.    mutual  occurrence  criterion  was  used  in 
the  LASA  program  with  good  resultsC93.   All  the  channel 
flag  information  is  transferred  to  the  active  event  matrix, 
and  the  channel  flags  are  reset. 

The  deactivate  old  event  module  was  used  to  remove 
events  which  were  past.   A  set  EVENT  DELAY  time  after  an 
e-.  ent  is  declared,  it  is  written  to  tne  outout  file  SL~d 
erased  from  the  active  event  matrix.   This  prevents 
spurious  peaks  from  being  added  to  a  event  long  past.   The 
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EVENT  DELAY  time  of  0.5  second  was  chosen  because 
inspection  o-f  the  known  events  revealed  that  most  channels 
tripped  within  a  0.8  second  period.   The  RESET  DELAY  (0.03 
sec)  plus  the  EVENT  DELAY  result  in  this  0.8  second  look 
-for  each  event. 

The  output  o-f  the  event  detection  program  is  an  output 
■file  which  contains  the  time  each  event  was  declared,  which 
channels  were  -flagged  and  the  tiiXie  delay  and  peak  amplitude 
-for  each  channel.   The  time  delays  Are    relative  to  the 
earliest  channel  -flag  time,  so  one  o-f  the  channels  always 
has  a  zero  time  delay. 

The  -first  version  o-f  the  detect  ion    program  was  written 
to  take  a  new  short  average  and  long  average  at  every  data 
point  (every  0.004  seconds).   This  program  took  4  to  8 
hours  to  search  a  20  minute  data  tape.   A  concession  to 
speed  was  made  and  the  program  changed  to  compute  averages 
at  every  -fourth  sample  point  (every  0.016  seconds).   This 
reduced  the  accuracy  o-f  the  time  delays  and  the  ability  ot 
the  program  to  pick  up  events.   The  RATIO  had  to  be  lowered 
in  order  to  get  the  same  detections  which  were  obtained 
previously. 

Studies  to  -find  the  best  signal-to-noise  RATIO  were 
conducted  several  times.   Development  of  the  LASA  detection 
system  had  revealed  that  a  7  dB  si  gnal  — to- noi  se  ratio  wa.s 
needed  -for  75/1  detecti  onC8  3 .   This  equates  to  a  5  to  1 
ratio  o-f  signal  power  to  noise  power.   Since  I  was  working 


with  presBLire  vice  power  I  used  2.2  as  my  starting  RATIO. 
This  RATIO  engul-fed  me  in  -false  alarms.   A  quick  study  was 
dene  around  the  2.4  level.   The  -first  10  minutes  of  tape 
4013  were  run  at  RATIOs  of  2.4,  2.45  and  2.48.   This  tape 
had  been  visually  examined  in  detail  previously,  so  the 
events  were  known.   The  results  Are    shown  in  Table  2—1. 
Also  shown  in  Table  2-1  Are    the  results  of  a  second  study, 
done  after  the  program  had  been  changed  to  average  less 
often. 

Table  2-1   Determination  of  the  Best  RATIO 

Average  taken  at  every  data  point 

RATIO  detection  rate  false  alarm  rate 

2 . 4  947.  257. 

2.45  757.  137. 

2.48  717.  67. 

Average  taken  at  every  fourth  data  point 

RATIO  detection  rate  false  alarm  rate 

2.3  767.  467. 
2 .  38                         717.  297. 

2.4  597.  337. 

2.5  597.  337. 


aetec^ion  rate  =  #  event  detections  /  #  of  known  events 
FA  rate  =  #  non-event  detections  /  total  detections 
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Natice  how  the  detection  rate  has  decreased  and  the 
"fslse  alarm  rata  increased  as  a  result  a-f  only  averaging  at 
every  -fourth  data  point.   Averaging  less  -frequently  means 
there  is  a  smaller  probability  that  the  data  paints  to  be 
averaged  will  all  lie  near  the  peak  amplitude  of  the 
signal.   The  signal  level  is  generally  lower  than  that 
detected  when  averaging  every  data  point,  and  a  lower 
si gnal — to-noi se  RATIO  must  be  usee  to  detect  the  same 
events.   But  when  the  signal-to-noise  RATIO  is  lowered  the 
-false  alarm  rate  increases. 

The  RATIO  o-f  2.3S  was  settled  on.   This  is  a 
compromise  which  gives  a  detection  rate  which  -finds  most 
high  and  medium  strength  events,  and  which  has  a  tolerable 
-false  alarm  rate.   Because  the  detection  rate  is  less  than 
1007.     (71"'.)  there  were  events  present  which  could  be  seen 
visually,  but  were  not  picked  up  by  the  detect ioti     progr^jTi, 
The  RATIO  could  have  been  adjusted  to  detect  all  events 
seen  visually,  but  at  the  cost  o-f  a  multitude  c-f  -false 
alarms.   The  RATIO  was  kept  at  2.38  and  used  for  the 
detection  o-f  all  data  tapes. 

The  -final  version  o-f  the  hdetect    program  read  digital 
data  from  a  framread  file,  detected  possible  events  using 
the  less  frequent  averaging  scheme,  and  supplied  the  event 
tine  and  channel  time  del avs  and  amplitudes  ta  an  output 
file.   Once  the  RATIO  had  been  satisfactorily  set  the 
program  was  used  to  search  the  FRAM  IV  ambient  noise  tapes 
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■For  events,  and  no  further  program  development  was  done. 
Weaknesses  in  the  program  were  subsequently  discovered,  but 
have  not  been  corrected. 

The  biggest  problem  is  the  Accuracy    o-f  the  event  time 
(the  time  when  50/C  o-f  the  channels  have  been  flagged).   An 
event  time  reported  by  the  event  detection  program  will  not 
exactly  match  that  found  by  plotting  the  time  series.   The 
times  are  usually  within  3  seconds  o-f  each  other,  hut    hav,^ 
been  off  by  as  much  as  13  seconds  in  one  case.   The  time 
difference  between  the  two  methods  is  greater  at  the  end  of 
a  tape,  and  is  likely  to  occur  after  a  particularly  strong 
event  has  taken  place  (though  there  were  times  when  tims 
discrepancies  developed  without  strong  events  present-  and 
also  many  strong  events  existed  which  did  not  induce 
discrepancies).   Typically,  there  might  be  no  time 
difference  at  the  -^irst  part  of  the  tape,  then  after  a 
strong  event  a  Lhree  second  discrepancy  would  be  seen  And 
tliis  would  be  consistent  until  the  end  of  the  tape. 
because  the  errors  did  not  appear  randomly  throughout  the 
tape,  and  because  they  developed  impulsively,  I  believe 
that  the  problem  lies  in  the  time  counter  of  the  event 
detection  program  becoming  offset  from  the  time  of  the  raw 
data,  perhaps  because  of  short  records  in  the  raw  data. 
The  event  detection  praqram  was  not  written  to  adjust  the 
time  keeper  in  the  event  of  a  short  record.   Correcting 
this  may  eliminate  the  time  discrepancy  problem. 
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It  has  already  been  mentioned  that  the  accuracy  o-f  the 
time  delays  deteriorated  when  the  program  was  changed  to 
run  more  quickly.   This  also  made  the  detection  of  the 
event  peak  amplitude  less  likely.   As  a  result,  in  order  to 
get  accurate  locations  and  source  strengths,  both  time 
delays  and  peak  voltages  had  to  be  taken  manually  -from  time 
series  plots  o-f  each  event. 

The  other  major  problem  o-f  the  event  detector  program 
is  that  it  does  not  discriminate  between  an  Arctic  noise 
transient  and  an  artifact  such  as  an  air  gun  blast  or  a 
reverberation  shot.   Short,  strong  signals  Are    reported  as 
possible  events.   Adding  this  discrimination  to  the  program 
is  the  next  step  in  improving  its  usefulness. 

Visual  Confirmation 


Visual  confirmation  was  required  for  all  possible 
noise  events  in  order  to  eliminate  artifacts,  false  alarms 
and  multiple  detections  of  the  same  event.   In  addition,  in 
a  few  cases  visual  confirmation  revealed  two  events  where 
there  had  only  been  one  detection. 

The  event  detection  program  was  designed  to  preclude 
the  need  for  plotting  a  time  series  of  each  event.   The 
output  of  the  program  contains  time  delay  and  voltage 
amplitude  information  which  can  be  used  directly  in  the 
location  program.   However,  because  of  the  decreased 
accuracy  of  the  time  delay  and  amplitude  information,  and 


because  o-f  the  event  time  discrepancy  mentioned  previously, 
it  was  necessary  to  plot  the  time  series  of  each  event. 

The  first  step  of  the  visual  confirmation  is  to  review 
the  tape  log  for  any  artifacts  that  may  have  occurred 
during  the  recording.   The  times  are  noted,  and  these  Are 
compared  to  the  event  times  given  by  the  detection  program. 
Then  a  time  series  of  the  artifact  was  plotted  to  determine 
which  detections  were  associated  with  it.   In  general,  an 
artifact  such  as  an  air  gun  blast  did  not  affect  detections 
for  over  20  seconds. 

The  visual  confirmation  portion  of  the  procedure 
evolved  from  a  very  limited  look  only  at  events  which  could 
not  be  located  with  the  detection  program  generated  time 
delays,  to  a  three  step  plotting  procedure  for  each  event. 
During  the  early  period  of  this  work  the  hdetect    program 
output  was  used  directly  as  the  input  to  the  location 
program.   The  location  program  used  the  time  delays  to 
determine  the  event's  location  in  space.   Those  events 
which  could  not  be  located  needed  a  closer  look,  and  so 
their  time  series  were  plotted.   The  plots  were  made  of  a  2 
second  period  including  the  event  time  given  by  the  hdetect 
program.   Often  there  was  no  apparent  event  in  this  time 
series  plot,  and  the  detection  was  declared  a  false  alarm. 
When  an  event  did  plot,  manual  time  delays  were  taken  and 
used  to  locate  the  event.   These  manual  delays  located 
these  events  with  better  accuracy  then  the  hdetect 
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generated  time  delays.   It  soon  became  apparent  that  the 
best  answers  would  be  obtained  by  taking  manual  time  delays 
o-f  al  1  events.   Trying  to  localize  the  events  with  the 
program  generated  time  delays  was  dropped  -from  the 
procedure,  and  the  -first  step  a-fter  getting  the  hdetect 
program  output  became  doing  a  2  second  time  series  plot  o-f 
each  event. 

Arter  a  dozen  tapes  had  been  analyzed  in  this  tTiannei" 
the  discovery  of  the  event  time  discrepancy  was  made. 
A-fter  plotting  a  dozen  events  right  at  the  time  shown  by 
the  program,  the  -final  two  dozen  event  o-f  tape  4009  all 
appeared  to  be  -false  alarms.   The  quality  o-f  the  tape  was 
good  (low  background  noise) ,  so  this  seemed  highly 
suspicious.   A  broader  search  o-f  the  time  around  each  event 
showed  that  the  final  two  dozen  events  were  not  false 
alarms,  but  were  events  with  times  5  seconds  dif-ferent  than 
those  indicated  by  the  program,  so  that  none  of  those 
events  had  shown  up  on  the  2  seccnd  ti.Tie  series  plats.   The 
method  of  visually  confirming  events  was  changed  so  that  a 
waterfall  plot  was  made  around  the  time  of  each  event. 
This  showed  the  exact  time  of  the  event,  and  helped  discern 
the  pattern  of  time  discrepancy  between  time  series  and  the 
hdetect    program.   Once  the  pattern  was  found  events  were 
easy  to  -find  and  false  alarms  could  be  noted.   The 
confirmed  events  were  then  plotted  with  2  second  time 
ser i  es . 


-31- 

The  final  step  of  visual  confirmation  was  simply 
plotting  and  replotting  the  events  to  the  proper  gain  so 
that  the  higher  voltage  ampl  i  tL?des  would  not  be  cut  off. 
These  peak  voltages  were  manually  taken  from  the  time 
series  plots  and  used  to  determine  source  strength. 

The  final  method  used  for  visual  confirmation  was: 

1)  Check  tape  log  for  artifacts,  and  eliminate 
those  from  further  analysis. 

2)  Plot  a  waterfall  time  series  around  each 
possible  event,  separate  real  events  from  false  alarms,  and 
find   true  event  times. 

3)  Plot  2  second  time  series  of  each  real  events 
adjusting  gain  to  keep  from  clipping  higher  voltages. 

Using  this  technique  certainly  reduced  the  false  alarm 

rate.   A  breakdown  of  the  detection  statistics  for  tapes 

that  had  been  examined  by  both  methods  is  found  in 
Table  2-2. 


Table  2-2   Detection  Statistics  for  Two  Visual 
Confirmation  Methods 


Artifacts       Events       halse  Alarms 


Qrignai  Method  16.97.         37.57.  45.67. 

w/  2  sec  plots 

Ul  t  i  mate  Method  1  6 .  27.  65 .  47.  1 8 .  47. 

w/  waterfall  plot 

2  sec  clot 


(Percentages  of  detections  classed  in  each  category) 


The  human  interpreter  was  a  necessary  tool  in  this 
scheme.   There  was  not  necessarily  a  one-to-one 
correspondence  between  events  and  detections.   There  were 
cases  where  a  strong  event  would  cause  multiple  detections, 
and  cases  where  two  events  occurred  at  the  same  time  and 
caused  only  a  single  detection.   In  some  cases  a  series  o-f 
detections  seemed  to  be  an  event  and  an  echo,  or  perhaps 
stra-f  .   This  would  be  counted  as  a  single  event. 

The  method  o-f  determining  whether  a  detection  was  a 
■false  alarm  or  a  weak  event  was  sometimes  difficult-   In 
general,  if  the  detection  program  indicated  a  possible 
event,  "something"  could  be  seen  on  the  waterfall  plot. 
The  detection  was  dismissed  as  a  false  alarm  if  no  pattern 
for  taking  time  delays  could  be  seen.    (Because  of  the 
shape  of  the  hydrophone  array  there  were  consistent 
patterns  of  time  delays  depending  upon  the  direction  to  the 
event.)   Presumably  the  false  alarm  rate  depends  upon  the 
training  and  attention  of  the  human  interpreter. 

Manual  time  delays  were  taken  from  an  arbitrary 
reference  to  the  crossing  of  the  largest  peak  to  peak 
amplitudes,  as  shown  in  Figure  2—6.   For  most  events  this 
was  clear,  but  for  weak  or  complex  events  some  intuition 
was  needed. 

Voltage  amplitudes  were  taken  as  the  maximum  peak 
voltage  in  the  event  signal.   All  were  taken  as  magnitudes 
regardless  of  sign. 
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Figure  2—6   Sketch  showing  point  o-f  measurement  for  event 
time  delays.   The  event  is  timed  at  its  zero  crossing 
between  the  largest  pair  o-f  positive  and  negative  peaks. 
(The  measurement/analysis  system  has  a  polarity  o-f  negative 
voltage  -for  positive  pressure.) 
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Al 1  types  of  noise  event  signatures  previously 
observed  by  DyerC71  were  seen  in  the  ambient  noise  tapes  I 
evaluated.   The  majority  o-f  events  were  pops  and  extended 
pops.   There  were  also  a  -few  whines  and  stra-f  events. 
While  signature  types  were  noted  in  general,  the  signature 
type  o-f  each  individual  event  was  not  recorded. 


CHAPTER  3 
LOCATION  OF  NOISE  EVENTS 

Event.  Location  Program 

The  program  used  for  localization  was  based  on  the 
program  FQUAK    by  Peter  Stein  C13D.   This  program  places  the 
event  at  di-f-ferent  trial  locations  and  computes  the  slant 
range  to  each  hydrophone.   Figure  3—1  shows  the  coordinate 
system  used  for  these  calculations.   These  s^re    plotted 
against  the  time  delays  and  a  least  square  fit  is  done  to 
determine  slope  as  shown  in  Figure  3—2. 


slope  =  ^^A:^ ^^  ^ =  A  (3-1) 


ZAt  -  (A^R) 

y  intercept  =  -±=l!=± 'n =  B 

N 


The  standard  deviation  of  the  time  delays  from  the  sloce 
line  is  fi cured. 


(At  -  AR  -  B)2 
Sigma  =^/ =   ^  (j.- o., 


N 


The  location  having  the  lowest  standard  deviation  is  the 
location  of  the  event.   The  inverse  of  the  slope  is  the 

group  speed  of  the  signal.   The  y  intercept  of  the  plat  ■ 
added  to  the  reference  time  of  the  manual  time  deiavs  tc 
get  the  time  the  event  actually  occurred  (as  opposed  to 
>jhen  it  reached  the  hydrophone  array). 
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Figure  3-1   Coordinate  system  used  -for  calculation  of  event 
location. 
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Figure  3-2   Least  squares  -fit  o-f  time  delays  and  test 
location  slant  range. 
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I  have  aBsumed  that  the  strongest  peak  pressure  sensed 
at  the  hydrophone  is  due  to  a  waterborne  acoustic 
propagation  path  -from  a  source  located  in  the  ice  sheet. 
The  signal  enters  the  water  near  the  source  and  propagates 
directly  toward  the  hydrophone. 

I  have  assumed  that  the  signal  does  not  bounce  o-F-f  the 
ocean  bottom  or  the  ice  canopy  be-fore  reaching  the 
hydrophone.   Paths  bouncing  o-f-f  the  bottom  would  produce 
signals  with  much  lower  energy  than  the  direct  path  signal, 
and  can  be  ignored.   Signals  bouncing  o-f-f  the  ice  canopy 
ArEi    too  energetic  to  ignore  but,  as  I  show  subsequently, 
they  do  not  a-f-fect  the  time  delay  computations  significantly, 

The  location  program  is  based  on  arrival  times  being 
related  to  slant  range,  R  ,  and  does  not  take  the  upward 
re-fraction  o-f  the  acoustic  path  into  account.   The  impact 
this  has  on  the  results  is  discussed  in  the  next  section. 

The  location  program  takes  as  input  a  file  of  time 
delays  and  voltage  amplitudes,  and  outputs  a  file 
containing  the  best  event  location,  sound  speed,  and 
standard  deviation.   It  also  computes  source  strength  based 
on  the  voltage  amplitude  inputs,  the  event  location  and  a 
spherical  spreading  loss.   This  feature  was  originally 
included  so  that  the  source  strp^.igth  could  be  computed 
directly  from  the  event  Jetection  program  outputs.   Since 
the  peak  vol  Lages  recorded  by  the  detection  program  a.re    not 
as  accurpLe  as  those  done  by  hand,  and  since  the 
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transmisBion  loss  does  not  -follow  simple  spherical 
spreading,  these  computed  source  strengths  were  not  used 
■for  any  part  of  this  study. 

The  FQUAK    program  set  up  a  grid  of  points  around  a 
specified  center  position.   The  grid  consisted  of  a  point 
every  100  meters  from  -5000  to  +5000  meters  in  both  the  k 
and  y  directions.   This  resulted  in  100  x  100  test 
locations.   When  the  best  test  location  was  found  the 
interval  spacing  was  reduced  to  every  10  meters,  and 
another  10,000  test  locations  were  generated  using  the  best 
location  of  the  first  round  as  the  new  center.   The  process 
was  repeated  with  a  1  meter  interval  to  get  the  final 
answer*   The  scheme  evaluated  a  total  of  30,000  test 
locations,  covered  a  range  out  to  5000  meters,  and  tack 
about  20  minutes  to  run. 

I  noted  that  a  significant  number  of  events  found  ^-jith 
FQUAK    were  at  the  range  limit  of  5000  meters.   The  program 
locatiori     '>jas  (c-iritten  to  search  a  larger  area  faster.   The 
fineness  of  the  grid  was  decreased  to  20  x  20  vice  100  x 
100.   A  1000  meter  interval  was  added  to  enable  the  program 
to  search  out  to  10,000  meters.   This  reduced  the  total 
number  of  test  locations  to  1600  (20  x  20  x4) ,  and  the  time 
to  one  minute.   locat iori    gave  answers  which  were  very 
ccnsi  stent  wi  th  FQUAK  ,  ax  oept  i  n  one  parti  cui  3.r    si  tuat  i  on  , 

The  wider  grid  size  led  to  one  problem.   The  location 
program  sometimes  found  the  lowest  standard  deviation  for  a 
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poip.t  in  the  quadrant  directly  oppasite  the  true  location. 
This  was  suggested  by  the  sound  speed  being  reported  as 
approximately  -1440  m/sec,  as  illustrated  in  Figure  3-3. 
This  problem  was  solved  by  modifying  location    to  make  the 
program  t  ine locate .        This  program  used  the  grid  size  and 
spacing  o-f  FQUAK,     and  centered  the  search  so  that  the  user 
could  designate  which  o-f  the  quadrants  would  be  searched. 
A  casual  leak  at  the  manual  time  delays  o-f  an  event  easily 
reveals  the  appropriate  quadrant.   This  program  works  well, 
but  is  as  slow  as  the  original  FQUAK.        It  was  used  rarely. 

As  with  the  original  FQUAK ^     I  began  to  notice  that 
some  events  were  located  at  the  range  limit  o-F  the  locat ion 
program.   This  led  to  the  modi -f  i  cation  o-f  the  location 
program  to  -Form  the  program  tar  locate .        This  program  uses 
the  locat ion    grid  size  and  fineness,  but  allows  the  center 
point  to  be  any  of  the  far  ccix-r',er'=.    of  the  original  location 
grid,  or    at  the  limit  range  at  each  of  the  cardinal  points. 
This  is  shown  in  Figure  3—4.   This  allowed  events  to  be 
located  out  to  20,000  meters. 

The  locat  ion    program  source  code  and  a  brief  user's 
manual  e^rs    found  in  Appendix  B.   This  program  was  written 
in  the  c  programming  language  for  the  UNIX  operating 
system.   This  program  was  developed  to  the  point  of 
usefulness,  and  then  used  to  locate  events.   No  further 
program  development  was  done  (except  the  very  minor  changes 
to  produce  t  ine  locate    and  t'arlo'zate)   ,  so  there  &r^    surely 


-41- 


(0 
V 

c 

0 
0 
4> 
0) 

0 

« 

D 

E 


4.2       4.4       4.6 


4.8         5         5.2 

(Thousands) 

Slant  Range  (meters) 


5.4       5.6        5. 


Figure    3— 
away . 


3-3       Least    squares    fit    to    a    -false    location     ISO 
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Figure  3—4   Areas  covered  by  locatiori,    tar  locate    and 
f  irie locate . 


-43- 

improvements  to  be  made. 

The  location  program  is  quite  interactive.   One 
hydrophone  with  a  bad  time  delay  can  change  the  slope  and 
location  a  great  deal.   An  event  is  located  by  eliminating 
bad  time  delays  and  checking  the  sound  speed  and  standard 
deviation  o-f  the  location.   In  some  cases  no  hydrophones 
needed  to  be  removed,  but  in  most  cases  at  least  one 
hydrophone  was  removed  be-fore  an  event  was  considered 
located.   The  sound  speed  was  the  major  indicator  of 
whether  an  event  had  been  located.   I-f  the  sound  speed  was 
between  1380  and  1500  m/s  the  event  was  considered  located, 
O-f  course  an  attempt  was  made  to  get  close  to  1440  m/sec. 
This  had  to  be  balanced  with  reducing  the  standard 
deviation.   A  standard  deviation  below  0.01  seconds  was 
considered  good. 

A  table  summarizing  all  o-f  the  events  and  their 
location  parameters  is  -found  in  Appendix  C.   The  standard 
deviations  (sigma)  are  given  in  two  sets  o-f  units.   The 
■first  is  the  sigma  calculated  by  the  location    program,  and 
it  is  in  seconds.   The  second  sigma  is  a  translation  of 
that  standard  deviation  to  meters  using  the  sound  speed 
calculated  for  each  particular  event.   The  standard 
deviations  ranged  from  0.0010  to  0.0327  sec,  with  0.0077 
sec  being  ths  average.   The  significance  of  this  standard 
deviation  will  be  discussed  in  the  next  section. 

In  some  cases  just  removing  suspect  time  delays  did 
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not  lead  to  a  localization.   A  reexamination  o-f  the  event 
time  series  was  done  to  see  i -F  any  of  the  manual  time 
delays  was  incorrect.   Often  a  reexamination  of  the  time 
series  produced  a  change  of  1  to  4  of  the  time  delays. 
These  corrected  values  plus  values  from  the  other  channels 
would  then  be  used  to  locate  the  event.   About  15>1  of  the 
events  required  reexamination.   Most  of  those  were 
subsequently  located. 

Despite  the  above  efforts,  there  were  a  few  events 
that  could  not  be  located  within  the  1330  to  1500  m/sec 
sound  speed  limits.   These  events  may  be  from  propagation 
paths  other  than  the  assumed  direct  acoustic  path.   Events 
arriving  primarily  through  the  ice  longitudinal  wave  or 
the  ice  flexural  wave  would  have  phase  speeds  above  and 
below  my  sound  speed  limits.   These  non— locatabl e  evvents 
Ar&    indicated  in  the  event  location  summary  of  Appendix  C, 
and  they  were  not  used  for  any  analysis  which  required 
Adcurats    location. 

Figure  3—5  shows  the  position  of  the  events  located 
within  a  2  km  square  centered  on  the  array  origin.   Figure 
3-6  shows  the  position  of  all  events  located. 
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Figure  3-5   Noise  events  located  within  a  2  km  squar* 
surrounding  the  array  origin. 
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Figure  3-6   Position  of  all  noise  events  located 
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Ef-fects  o-f  Retraction  on  Location 

A  sound  speed  pro-fiie  was  used  to  get  the  re-fractive 
paths  -for  various  ranges.   This  was  simpli-Fied  by  the  fact 
that  all  the  hydrophones  were  at  a  depth  o-f  93  meters. 
Assuming  that  only  the  "direct"  path  is  involved  means  that 
each  horizontal  range  has  to  have  a  unique  launch  angle  in 
order  to  reach  the  hydrophone  at  its  speci-fic  depth.   Rays 
were  launchad  into  the  layers  at     the  sound  velocity  prcfile 
and  the  horizontal  range  to  the  hydrophone  was  calculated. 
The  time  required  to  travel  the  refractive  path  can  be 
calculated  and  compared  to  that  o-f  the  slant  range.   This 
time  error  can  then  be  related  to  the  error  of  the  location 
program. 

Figure  3—7  shows  the  linearized  sound  velocity  profile 
that  was  used.   It  is  based  on  the  sound  velocity  profile 
reported  for  the  eastern  Arctic  ocean  by  ChenCll.   Figure 
3-3  helps  to  illustrate  the  scheme  used  to  calculate  the 
ray  paths.   Equations  3.4,  3.5  and  3.6  were  used  to 
calculate  angles,  ranges,  depths,  and  propagation  time. 


V, 


g  cos  8 


cos  Q„  —  cos  Q 


( 3-4 ) 
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(1  +  sin  Qj)  (1  -  sin  8^^) 
(1  -  sin  Bj) (1  +  sin  Q^) 
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Figure    3-7       Sound    velocity    pro-file    used    in    predicting 
re-fractive    paths. 
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Figure  3—8  Scheme  used  to  compute  re-fractive  propagation 
paths.  Source  is  located  at  z  =  Om  and  hydrophone  at  z  = 
93m.   Sound  speed  gradients  change  at  SOm,  254m  and  362m. 


-so- 
using the  known  depths  (z)  and  estimated  speed 
gradients  (g)  o-f  Figure  3-7,  and  choosing  a  particular 
launch  angle  O^^)  ,  all  o-f  the  subsequent  angles  o-f 
intersection  o-f  the  layer  interfaces  iB^,    8^; ,  84)  can  be 
-found  -from  Equation  3-4.   The  angle  which  intercepts  the 
hydrophone  at  its  depth  O^)  can  also  be  found.   Uiith  the 
angles  known,  the  horizontal  range  and  propagation  time  far 
each  layer  can  be  determined  with  Equations  3-5  and  3-6, 
These  &r&    combined  to  get  the  total  horizontal  range  and 
propagation  time.   It  should  be  noted  that  there  ar&    two 
ranges  and  times  for  each  launch  angle.   The  first  path  is 
that  which  intercepts  the  hydrophone  on  the  way  down,  while 
the  other  intercepts  the  hydrophone  as  it  is  refracted  back 
toward  the  surface.   The  maximum  depth  reached  by  the 
propagation  path  was  also  found,  and  those  paths  that  went 
below  754  meters,  resulting  in  a  range  greater  than  30600 
meters  were  not  reported.   A  tabular  summary  was  made  cf 
launch  angle,  horizontal  range  from  the  hydrophone,  maximum 
depth  and  propagation  time,  and  this  may  be  found  in 
Appendix  D,  along  with  more  detailed  tables  listing  Qj^_4, 
r^_5  and  t^.^. 

Rays  connecting  source  and  hydrophone  with  one  or  more 
bounces  from  the  ice  were  not  considered  here.   The  effect 
of  those  rays  will  be  taken  into  account  in  Chapter  ^. 
The  refractive  propagation  time  was  less  than  the 
slant  range  propagation  time  because  the  refracted  path 
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travels  through  faster  water.   The  slant  range  propagatian 
time  was  calculated  by  dividing  the  slant  range  by  1438.48 
m/sec,  the  average  sound  speed  between  O  and  93  meters 
depth.   The  time  di-f-ference  between  the  slant  range  path 
and  the  re-fracted  path  Ars    shown  in  Figure  3-9  as  a 
-function  o-f  horizontal  range. 

This  time  difference  is  greater  than  the  average 
standard  deviation  of  the  I ocat i oti    program  only  after 
13,000  meters,  and  the  time  difference  at  20,000  meters  is 
only  about  twice  that  average.   The  standard  deviation  does 
not  reflect  the  time  difference  due  to  refraction  because 
all  of  the  hydrophone  time,  delays  a.re    adjusted  in  the  same 
manner  and  direction.   Figure  3-10  shows  that  sigma  does 
not  grow  with  horizontal  range.   Refraction  effects  dc  not 
influence  the  standard  deviation  greatly.   Closer  than 
13,000  meters  the  range  errcr    caused  by  other  factors  rr-iasks 
any  error  from  ignoring  refraction. 

There  is  a  better  point  of  focus  for  examining  the 
effect  of  refraction,  and  that  is  the  change  in  time  delay, 
not  the  change  in  the  propagation  time  itself.   A  point  was 
chosen  at  approximately  5  km  from  the  origin  of  the 
hydrophone  array,  and  another  chosen  at  approximately  6  km. 
The  slant  range  propagation  times  and  refractive 
prcpagation  tiines  were  calculated  for  each  point.   The  time 
delay  between  these  two  points  was  0.6857  sec  for  the 
refractive  oath  and  0.6930  sec  for  the  slant  range  oath. 
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Figure  3-9   Time  di -f -f  erence  between  the  slant  range 
acoustic  path  and  the  re-fractive  path  as  a  -function  o-f 
horizontal  range. 
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Figure  3—10   Scatter  plat  o-f  standard  deviation  and  range 
■for  all  located  noise  events. 


-54- 

The  di-fference  was  0.0073  (about  IX),  or  a  deviation  o-f 
0.OO37  5ec  -far    each  hydrophone.   Two  points  at 
approximately  18  and  20  km  were  also  evaluated.   The 
difference  between  their  time  delays  was  0.0182  sec  (about 
0.  IX)  ,  or  0.0091  sec  per  hydrophone.   These  numbers  Bre:    the 
same  order  as  the  total  error  of  the  location  program. 

The  location  program  may  compensate  for  some  of  this 
error  by  raising  the  sound  speed.    If  just  the  points  above 
were  used,  the  sound  speed  would  go  from  1438.5  to  1453.3 
at  5  km  and  to  1458.6  at  20  km.   With  24  time  delays  being 
used  in  the  location  program  the  effect  may  not  be  as 
great. 

The  main  source  of  error  in  the  location  program  is 
the  quality  of  the  manual  time  delays.   When  the  signal-to— 
noise  ratio  was  low,  picking  the  correct  peak  was  often 
difficult.   The  standard  deviation  will  reflect  the 
judgement  of  the  person  picking  off  the  time  delays.   The 
time  delays  were  only  measured  to  the  closest  0.003  sec. 
It  is  interesting  to  note  that  0.006  seconds  equates  to  the 
width  of  a  pencil  tip  on  the  time  series  plot  scale. 

The  final  question  to  be  answered  is  "How  do  the 
standard  deviation  and  refraction  errors  equate  to  the 
range  and  bearing  acourAczy    of  the  location    program."   Two 
hypothetical  noise  events  were  investigated,  one  at  5000 
meters  (2845,  4136)  and  the  other  at  20,000  meters 
(-8253,  18896)   The  time  delays  for  slant  range  propagation 


and  refractive  propagation  were  calculated.   The  location 
program  was  run  for  each  set  of  time  delays,  and  for  each 
set  partially  contaminated  with  0.016  sec  errors.    (Zero, 
+0.016  and  -0.016  were  each  added  to  one-third  of  the  time 
delays.)   The  results  are    summarized  in  Table  3-1. 

Table  3-1   Results  of  locat ion    Program  Accuracy  Test 


R 

Ar 

(m) 

(m) 

lant  range 

5  km 

4725 

-295 

20  km 

18477 

-2143 

ef raction 

k 

5  km 

4724 

-296 

20  km 

18477 

-2143 

1  ant  range  w/ 

5  km 

4191 

-829 

20  km 

20526 

-94 

©fraction  w/ 

5  km 

4209 

-811 

20  km 

20526 

-94 

6    Aj*      CT 


c 


(dieg)   (deg)     (sec)   (m/s) 

55.6   +0.1    0.0003    1440 
113.7   +0.1    0.0002    1437 

55.6   +0.1    0.0003    1454 
113.7   +0.1    0.0002    1458 

55.8   +0.3    0.0127    1453 
113.9   +0.3    0.0130    1443 

55.8   +0.3    0.0130    1466 
113.9   +0.3    0.0130    1463 

The  refraction  contaminated  by  errors  case  is  closest 
to  what  was  input  into  the  locat  ion    program  for  the  field 
events.   This  table  gives  an  estimate  of  the  accuracy  of 
the  locat  ion    program  as  800  m  at  5  km,  and  2000  m  at  20  km. 
The  bearina  accuracv  is  excellent. 
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CHAPTER  4 
STRENGTH  OF  NOISE  EVENTS 

Acoustic  Source  Model 

The  dipole  is  considered  a  possible  source  model. 
Peak  values  -for  the  source  parameter  o-f  -force,  F  ,  a.re 
used. 

The  acoustic  pressure  due  to  a  non— convecting  compact 
diools  source,  in  a  nonref  r  act  i  nq  in-finite  medium.  IeCSj: 


P  = 


►  in  9 

47TR 


i    5^  -.  ^ 
c  ^t    R 


(4-1) 


where  R  =  slant  range. 

Figure  4—1  shows  the  orientation  o-f  the  presumed  dipole. 
The  angle   a   is  the  launch  angle  -from  the  horizontal  plane 
down  into  the  water. 

Assuming  that  F  may  be  expressed  as  a  harmonic. 
\0   P'''Ot{  =GJF  =  zTf-fF  ..     the  pressure  may  then  be  expressed 


IpI  = 


sin  8 
47fR~ 


2  7T-fF    F 

+  _ 

c      R 


(4-2) 


Solving  -for  F  gives 


p,3  471  R^c 
°    sin  Q  L27f-fR 


KH^  —  -j.'. 


+  c3 


In    the    -far    -field    the    2  7T"fR    term    dominates    the    sum    in 
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Figure  4-1   Assumed  dipole  s 
de-finition  o-f  launch  angle,   8 


ource  orientation  and 


-58- 

the  denominator.   For  the  lowest  -frequency  considered  in 
this  study  (20  Hs )  ,  c  (1440  m/sec)  is  10/1  o-f  ZTT+R  at  115 
meters  and  only  1.17.  o-f  2Tri^R  at  1000  meters.   When  the  c 
in  the  denominator  is  neqlected  the  -force  can  be  written 


-PoP-X 


°    sin  8 


where   ^  is  the  naveiength 


The  peak  pressure,  p^.^  ,  should  lead  to  the  peak  -f 


orce 


F   .   This  de-finition  o-f  -force  was  used  as  the  parameter 
-for  dipole  strength.   Event  signatures  that  were  recorded 
-from  a  source  within  300  m  o-f  a  hydrophone  were  not  used  to 
calculate  dipole  strength,  F   ,  ^rom  peak  pressure,  p   . 

For  this  model  the  peak  acoustic  pressure  must  be 
-found.   The  hydroohone  sensitivity  o-f  -15^'  dB  re  t  volt  per 
1  //Pa  --jas  used  to  convert  vcltage  to  pressureL  1 7  j  . 

1  volt  =>  89  h4/m^  =  89  Pa   .  (4-5) 

The  dicole  strength  -formula  requires  wavelength,  \ 
Frequency  was  taken  -from  the  time  series  plots  -for  each 
event  via  axis  crossing  rate,  and  \      was  determined  by 
dividing  c  (1440  m/sec)  by  the  -frequency. 

Launch  ar=gle  needed  -for  the  dipole  model  can  be  -fauna 
as  in  Chapter  3  by  assuming  a  sound  velocity  pro-file  and 
computing  the  re^^r active  path.   There  is  a  uniqsje  launch 


angle  iar    each  horizontal  range  when  the  path  is  purely 
re-fractive,  but  a  range  o-f  launch  angles  when  sur-Face 
reflection  paths  are  included. 

The  final  parameter  in  the  dipole  strength  -formula  is 
slant  range.   The  answers  obtained  using  slant  range  Are; 
the  strengths  based  on  spherical  spreading  in  a  non- 
absorptive  medium,  equation  4—4.   Because  the  spherical 
spreading  assumption  is  a  poor  one,  '.:  r  e-f  ract  i  ve  and  surface 
reflective  propagation  paths  Are    caused  by  the  Arctic  sound 
velocity  profile),  equation  4—4  must  be  modified.   The 
effect  of  refraction  on  spreading  loss  will  be  discussed  in 
the  next  section. 

Volumetric  absorption  was  found  by  using  the 
absorption  formulas  of  DyerCSD.   Assuming  a  pH  of  3.2,  a 
salinity  of  33.5  °/oo,  a  temperature  of  0°  C  and  a  pressure 
of  40  atmospheres,  I  calculated  the  total  volumetric 
absorption  to  be  1.3  x  10  ""  dB/km  for  an  SO  Hz  signal.   For 
my  maximum  horizontal  range  of  20  km,  the  absorption  weald 
be  0.026  dB.   This  is  not  significant,  and  I  therefore  did 
not  include  a  volumetric  absorption  correction  in  the 
strength  calculations. 

Effects  of  Refraction  on  Transmission  Loss 

Spherical  spreading  loss  in  a  nonref ract i no  medium  is 
illustrated  in  Figure  4-2C143.   The  sound  pressure  is 
pressumed  to  soread  radial  1/.   The  sound  oressure  squared 
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Figure  4-2   Spherical  spreading.    (From  UrickClAU) 
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i=>    proportional  to  intensity,  and  intensity  is  the  power 
per  unit  area.   Since  the  power  -from  a  source  is  constant 


P  =  I.aTTRi^  =  Ir^^TT^rT       .  (4-6) 


The  intensity  at  the  reference  range  of  1  meter  can  be 
related  to  other  intensities  by 


If.  =  =  .  (4-7) 

47TR^      R^ 


Since  I  =  p^/Qo     ,  this  can  be  expressed  in  terms  o-f 
transmission  loss,  H  . 


H  =  -10  log =    10  log  R-^  =  20  log  R   ,  (4-8) 

P__ r ^ ' 

r  eT 
in  dB  re  the  distance  re-ference,  taken  as  1  m. 


The  spreading  scheme  -for  a  refractive  medium  is  shown 
in  Figure  4— 3C2D.   This  is  based  on  ray  theory  which 
assumes  that  acoustic  energy  does  not  cross  the  rays,  with 
energy  contained  between  two  rays  being  conserved.   The 
intensity  at  the  re-ference  range  between  the  two  rays  shown 
i  s: 


P  P 

I  =  7^ —  =  ^= A —   -        (4-9) 

27rRcos  8   R^8    27rcos  a   ^Q 


At  a  horizontal  distance   r  meters  from  the  source,  the 
intensity  is: 
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Figure  4-3   Spreading  loss  in  re-fraction. 
MedwinC2:) 


(From  Clav  and 


-ib.j. 


P  P  P 

I       =    =    =    .  (4—10) 


The    relation    between    intensities    becomes: 

Ve-F^TTcos    9^   A3  I^^^    ABcos    9^ 

I       =    — ^— , t=V   =    x-^ r      •  (4-11) 

^  ^TTr  ^rjsin    Q^  |  r^r|5in    Q^l 

The    loss    due    to    spreading    is: 


AQcos 


V  AQ^°^  ^o 


"re 


-f  r^r  I  sin 


(4-12) 


o^o 


In    terms    o-f    pressure,     I^^    =    Pt-'^/piCi?     -^f^d     1^-^^    =    Pret^'Pc 
and    there-fore 


Pr"  Pl^l         A®    '^°^    ^o 

^-^  .  (4-13) 


Ve^^  /^°"°    r^r|sxn    B, 


bince    b V    Snel  1  's    lew,        cos    B_/c_.    =    cos    S 


I-  1 


Pr"       _     Pi      If    '=°^    Ql  ^  /\q 

^;^    ~    P^   r^r|sin    Q^l     "    r  ^r  |  tan    sj        ' 

since  Pi^'Pa    =  1  i^  seawater  to  an 

excellent  approximation. 


Applying  the  dipole  model  to  this  spreading  los; 
SGuation  gives 


z  -  p.   ^ 


o   ^ 


P   ^  =   , r-  *— * 1  ',4—15,' 


wher=  A  =  F^_j/2X 
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Equation  4—15  assumes  a  unique  refractive  path  between 
source  and  hydrophone.   In  the  Arctic  there  may  be  other 
paths  due  to  the  non— specular  scattering  o-f  rays  o-f-f  the 
ice  canopy.   Figure  4-4  illustrates  how  rays  normally 
trapped  in  a  sur-face  duct,  may  be  deflected  down  to  a 
hydrophone.   The  minimum  verte>;ing  angle  calculated  from 
the  linearized  sound  velocity  profile  of  Chapter  3  was 
0.064  radians.   For  a  ray  to  stay  in  a  surface  duct  above 
the  hydrophone  it  must  be  reflected  from  a  slope  of  less 
than  0.032  radians  or  about  2°.   It  is  reasonable  to  assume 
that  the  ice  canopy  lacks  local  levelness  to  this  order,  so 
that  non— specular  rays  must  be  accounted  for. 

The  rays  which  rebound  from  the  ice  canopy  experience 
some  loss.   The  attenuation  for  the  FRAM  IV'  e>jperiment  has 
been  reported  at  0.1  dB/km  at  80  Hz C 1 1 3 .   This  attentuation 
mav  be  converted  to  a  loss  oer  bounce. 


P     =    O. 1    dB/km    =    -       ,  (4-16) 


where  b  =  loss  per  bounce,  and 
X  =  cycle  distance. 


The  cycle  distance  depends  on  the  launch  angle  and  the 
sound  speed  gradient.   For  a  launch  angle  of  0.032  radians 
and  the  assumed  sound  velocity  profile  of  Chapter  3,  the 
cycle  distance  is  3.7  km.   Therefore,  the  loss  per  bounce 
is  about  0.4  dB.   This  loss  is  lov^a  enouch  that  even  a  ray 
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Figure  4-4   Specular  and  non— specular  reflections  -from  the 
ice/water  inter-face. 
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which  has  bounced  several  times  may  contribute  a 
signi-ficant  amount  of  energy  at  the  hydrophone.   The  non- 
specular  rays  cannot  be  ignored.   Bounce  loss  at 
-frequencies  less  than  80  Hz  Ars    even  smaller,  since  the 
data  show  a  roughly  linear  dependence  on  -frequency. 

To  account  -for  the  non— specular  rays  the  spreading 
loss  is  calculated  using  the  ray  averaging  techni queC5] . 
The  pressure  -from  a  particular  ray  at  a  given  depth  and 
horizontal  range,  assuming  a  dipole  source  model,  is 


rt  o   da  dr 

p'"(r,z)  =  i , .  (4-17) 

r I  tan  Q^ |  dr  X/2 

dr 

The  term  represents  the  probability  that  a  rav  bundle 

X/2 

will  cross  a  certain  depth,  as  shown  in  Figure  4—5.   For  a 
single  linear  sound  speed  gradient  the  cycle  distance  can 

be  written  as 


X  =  sin  a^.^  =  2r^  sin  S^  ^  ^'"c^q   '       •;4-lS: 


where  r   is  the  radius  o-f 


curvature,  to  a  good  approximation  constant 
-for  all  small  angle  rays  in  a  linear  sound 
speed  gradient. 


Applying  equation  4—18  to  equation  4—17,  and  using  the 
small  angle  approximation,  gives 


A^  Iq^I  d8„ 
r>  '  o  I     o 

p'-(r,z)  =  .    I     .  (4-19) 

r^  r  (8  J 
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Figure  4-5  Probability  o-f  a  ray  bundle  croBsing  a  certain 
depth  at  a  given  horizontal  range. 
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In  order  to  average  the  contributions  of  the  possible 
rays,  this  pressure  is  integrated  over  all  possible  angles 
-for  a  given  receiver  depth,  and  then  averaged  over  depth 
doMn  to  the  hydrophone  at  z 


o' 
z. 


Q. 


p^(r)     = 


2    A^    1 


Q, 


r    r 


c       o 


dz 


O  8 


a 


dS 


o        ' 


(4-20) 


m 


where    Q,,    is    the    maKimuiTi    launch 


v 


angle  of  a  ray  that  will  hit  the  hydrophone 
at  a  given  range,  and  Q   is  the  minimum 
launch  angle. 


The  angle  S^  is  a  function  of  8   and  z 


■   Ccos  Sh  —  cos  Qq] 


(4-21) 


Using  the  small  angle  approximation  for  cosine  leads  to 


®1  =-^/9o" 


( 4—22 ) 


Substituting  this  into  equation  4-20  and  evaluating  the 
integral  over  angle  gives 


p'^  (r )  = 


2  A^ 


r  r 


c   o 


OL 


'^v^  -  .-   -  V^"'   -  -. 


dz   .   (4-23) 


a  '^  =  —  for  all  z  so  the  second  term  within  the  integral 

^c 
IS  always  zero.   Evaluating  the  first  term  o\'&r    depth  gives 

an  expression  for  pressure  in  terms  of  r  and  8^. 


D^ir)     = 


A'=^ 


^  r  z, 


-69- 


Q  ^'  - 

V 


W 


,,  n3/2 


r 


c  J 


(4-24) 


This  expression  can  be  used  to  find  the  source  strength. 


o  r  z 


A^  =  p-^ 


8. 


,  t3/2  , 


Q  ^  - 

V 


c  J 


=  G(r)p'^   .     (4-25) 


The  spreading  -function,  G,  is  presented  as  a  -function  ot  r 
alone  since  8^^  depends  on  r.   For  each  r  there  is  a  unique 
Q^ ,  and  there-fore,  a  unique  G.   The  spreading  function  was 
calculated  for  horizontal  ranges  from  300  m  to  20,000  m, 
and  tablulated  in  Appendix  D.  -The  spreading  function  is 
shown  in  a  log— log  plot  in  Figure  4—6.   For  comparison  the 


equivalent  spherical  spreading  for  a  dipole  source  is 


I  1  •— li  I 


shown.   From  this  one  can  see  that  source  strengths 
calculated  using  the  spherical  spreading  law  lead  to  an 
unrealistic  dependence  an  range. 

In  order  to  get  8.^^  and  r   a  linear  sound  speed 
gradient  of  0.054  sec    was  chosen.   This  gradient  gives 
the  same  S   at  r  =  3  km  as  the  multiple  step  profile  used 
in  Chapter  3.   Three  kilometers  was  chosen  since  it  was  the 
median  horizontal  range  for  the  noise  events. 

The  spreading  loss  function  and  measured  peak  pressure 
iTiaonitudes  were  used  to  calculate  dioole  strenath. 


F^  =  2Xa  =  2Xp^^G(r) 


(4-26) 
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Figure    4-6       Spherical      (heavy    line)     and    re-fractive    spreading 
loss,     E(r)  ,     as    a    -function    of    horizontal    range. 
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For  a  particular  event  the  source  strength  was 
calculated  from  each  hydrophone,  and  then  the  average  taken 
as  the  source  stren'^th  -for  the  event.   The  standard 
deviat'on  within  each  event  ran  from  9/C  to  125X  of  the  mean 
valu^.   An  event  strength  summary  in  Appendix  C  lists  the 
mean  measured  peak  pressure  and  the  mean  dipole  strength, 
along  with  the  standard  deviation  on  these  values  for  each 
event . 

Strength  of  Background  Noise 

I  was  interested  in  the  effect  that  environmental 
loading  might  have  had  on  the  temporal,  spatial  and 
strength  statistics.   It  has  been  shown  by  Makris  and 
Dyer C ion  that  low  frequency  (10—20  Hz  band)  ambient  noise 
rms  pressure,  averaged  over  a  long  time,  correlates  well 
with  environmental  stresses  and  moments.   Since  I  had 
ambient  nci se  pressure  for  most  of  the  period  of  the  FRAM 
IV  e>:  peri  men  t,  and  since  I  had  environmental  stresses  and 
moments  available  for  only  a  part  of  the  time,  I  chose  to 
use  the  20-BO  Hz  1 ong—ti me— average  rms  pressure  as  my 
environmental  indicator. 

The  10  to  20  Hz  band  ambient  noise  pressure  was 
converted  to  20  to  80  Hz  band  pressure  in  the  following 
manner.   Figure  4-7  shows  the  typical  spectrum  for  central 
Arctic  pack  ice  noise.   The  portion  of  the  spectrum  between 
10  and  100  Hz  can  be  approximated  by  a  straight  line. 
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FREQUENCY,  Hz 


Figure  4-7   Composite  central  Arctic  ambient  noise  spectrum 
observed  durinq  the  FRAM  IV  experiment.    (From  DYerC6:) 
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log  S  =  A  C  lag  -f  1   +   B   ,  (4-27) 

where   A  =  slope      =  —1.7273  Pa'^/Hz^   , 
B  =  intercept  =  -1.0909  Pa^/Hz   , 


or 


S  =  10^-f^ 


(4-28) 


The  band  rfrts  oressure  relates  to  the  Boectral  level 


P  rms,b 


/ 


S  d-f 


(4—29) 


I  have  asBLimed  that  as  the  sound  pressure  level  changes 
-from  time  to  time  the  intercept  'B  changes,  but  the  slope 
remains  the  same.   By  substituting  equation  4-28  into  4—29 
and  using  the  known  10  to  20  Hz  ambient  noise  band,  B  can 

be  fe-jritten  in  terms  of  the  known  oressure. 


B  =   1  OG 


P^rms. 10-20 


K. 


(4-30) 


where 


•^'l  = 


^q(A  +  1)  _  10 (A  +  1) 


A  +  1 


=0.1 020  Hz 


The  ambient  noise  rms  pressure  -for  the  20  to  80  Hz  band  may 
now  be  -found. 


80 


P^rms, 20-80 


-  f 


10^-F^  d-f 


(4-31 ) 
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rms, 10-20 


where 


gQ(A  +  1)  _  20^*^  "^  ^^ 


A  +  1 


=  0.0988  Hz    , 


or  final Iv, 


Prms , 20-80     Prms , 1 0-20 


(4-32) 


Prms, 20-80   ~  0.98  P^ms, 10-20 


Thus  the  band  -from  20  to  80  Hz  is  virtually  identical  to 
the  one  -from  10  to  20  Hz  in  rms  pressure,  tor  long-time- 
averages,  and  in  turn,  is  an  acceptable  surrogate  for 

envi  roRmental  forcing  (applied  stresses  and  moiTients )  .   The 
20  to  30  Hz  band  ambient  noise  rms  pressure  for  each  of  th« 
tapes  investigated  is  found  in  Table  4—1. 
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Tabla  4-1   20  to  SO  Hz  Band  Ambient  Noise  rms  PreBsure 


Tape  #             Date  Recorded  Prms  '^'O-BO 

'(Pa) 

4001  3-27-32  Not  Available 

2001  3-29-82  Not  Available 

2009  3-30-82  0.022 

3001  3-31-82  0.019 

4003  4-01-82  0.044 

4005  4-01-82  0.035 

4007  4-01-82  0.022 

4009  4-02-82  0.010 

4011  4-02-82  0.010 

4013  4-03-82  0.013 

2023  4-08-82  0.037 

4015  4-09-82  0.040 

3047  4- 1 3-82  0.010 

4016  4-15-82  0.017 
4019  4-15-82  •  0.016 
4021  4-19-82  0.013 

4023  4-19-82  0.011 

4024  4-19-82  0.011 

4027  4-20-82  0.012 

4029  4-20-82  0.012 

4031  4-20-82  0.012 

4033  4-20-82  0.012 

4040  4-21-82  0.034 

4047  4-21-82  0.114 

4049  4-21-82  O. 140 

4051  4-21-82  0.140 

4053  4-22-82  O . 080 

4055  4-22-82  0.082 

4057  4-22-82  0.053 

4059  4-22-82  0.065 

4061  4-22-82  0.034 

4063  4-22-82  0.023 

4065  4-22-82  0.027 

4067  4-22-82  0.018 
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CHAPTER  5 
ANALYSIS  OF  NOISE  EVENTS 

Detection  Analysis 

A  total  o-f  34  tapes  was  examined,  -for  a  total  time  of 
662  minutes.    (For  a  few  o-f  these  tapes  the  entire  20 
minutes  was  not  used. ) 

There  was  a  total  o-f  499  detections  o-f  events  flagged 
on  at  least  50'/.  of  the  hydrophone  channels.   Of  these,  139 
were  man-made  artifacts,  and  125  were  false  alarms 
(detections  which  were  so  weak  that  no  pattern  for  taking 
time  delays  could  be  discerned).   There  were  199  unique 
events,  and  36  multiple  dtections  of  those  events.   Stated 
in  another  way,  of  the  detections  which  were  not  artifacts, 
65.3/^  were  strong  enough  to  support  analysis  and  34.7yC  were 
too  weak  to  reasonably  analyze,  and  hence  labeled  false 
al arms. 

Since  the  detection  process  depends  on  si gnal — to— noi se 
ratio,  the  level  of  background  ambient  noise  should  affect 
the  event  detection  rate.   Figure  5—1  shows  normalized 
ambient  noise  pressure,  number  of  false  alarms  per  tape, 
and  number  of  unique  events  per  tape  for  each  tape 
examined.   There    is    some    trend    for    more    euents    being    round 
Mhen     the    ambient    pr essure    is    Ion,     and    more    false    alarms 
dec  1  ar  ed    t<jhen     the    amb  ien  t    pr  essur  e     is    high. 

This  is  more  clearly  seen  in  Figure  5—2,  which  shows 
the  average  number  of  false  alarms  and  unique  events  found 
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Figure  5-1   Normalized  ambient  noise  pressure,  number  o-f 
unique  events  per  tape,  and  number  o-f  -false  alarms  per 
tape  -for  each  data  tape  examined. 
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Figure  5-2  Average  number  o-f  -false  alarms  and  uniaue 
events  per  data  tape  -for  -four  ranges  o-f  ambient  noise 
pressure. 
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per  tape  in  each  o-f  -four  background  noise  pressure  ranges. 
The  0.01-0.02  Pa  range  used  15  tapes  to  compute  its 
average,  the  0.02-0.03  Pa  range  4  tapes,  the  0.03-0.04  Pa 
range  5  tapes,  and  the  over  0.04  Pai    range  8  tapes.   Two 
tapes  were  recorded  during  the  first  few  days  of  the  FRAM 
IV  experiemnt,  before  the  10-20  Hz  band  ambient  noise 
recordings  were  started. 

A  breakdown  of  detections  for  each  tape  is  found  in 
Appendix  C. 

Temporal  Analysis 

The  interarrival  time  between  events  ranged  from  1  to 
1064  seconds.   Each  event  time  was  taken  to  the  nearest 
second,  and  no  events  were  taken  as  having  the  same  event 
time.   If  two  events  happened  in  the  same  second,  one  was 
judged  to  be  earlier,  and  the  two  events  were  given  event 
times  one  second  apart.   The  interarrival  time  for  a 
particular  event  was  measured  from  the  previous  event, 
except  for  the  first  event  of  a  tape,  which  was  measured 
from  the  start  of  the  tape. 

The  interarrival'  times  were  divided  into  bins  of  20 
seconds.   The  first  bin  ("0")  contained  events  which  had 
interarrival  times  from  O  to  19  seconds,  the  second  bin 
from  20  to  39  seconds,  and  so  on.   The  number  of  events  per 
bin  is  presented  in  Table  5—1  and  shown  graphically  in 
Figure  5-3.   A  complete  listing  of  interarrival  times  for 
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each event  is  found  in  Appendix  C.   The  mean  of  the 
interarrival  times  </i.  )  is  100  seconds,  and  the  standard 
deviation  <  0")  166  seconds.   In  terms  o-f  bins,  the  mean  i 
5  and  the  standard  deviation  8.   The  standard  deviation  i 
1.66  times  the  mean. 

Table  5—1   Number  o-f  Events  per  Interarrival  Time  Bin 


Bin 


Events 


0 

67 

1 

36 

2 

19 

3 

15 

4 

12 

5 

7 

6 

4 

7 

3 

8 

6 

9 

1 

10 

3 

11 

4 

12 

2 

13 

1 

14 

3 

IS 

2 

16 

0 

17 

0 

18 

0 

19 

0 

20 

1 

21 

0 

22 

0 

23 

2 

24 

25 

26 

27 

2 

28 

29 

30 

0 

in 

Events 

31 

0 

32 

0 

33 

1 

34 

0 

35 

0 

36 

0 

37 

0 

38 

0 

39 

0 

40 

0 

41 

1 

42 

0 

43 

0 

44 

0 

45 

0 

46 

0 

47 

0 

48 

0 

49 

1 

50 

0 

51 

0 

52 

0 

53 

1 

54 

0 

55 

0 

56 

0 

57 

0 

58 

0 

59 

0 
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Three  dif-ferent  probability  density  -functions  were 
investigated  to  -find  an  appropriate  -fit  -for  Figure  5-3. 
They  were  1)  a  hal f -gaussi an  distribution,  2)  an 
exponential  distribution  and  3)  a  J  shaped  distribution. 

The  hal -f  — gaussi  an  probability  density  -function  isC3!]: 


p(t)  = 


-t-/2t^- 


V^^o 


(5-1) 


The  general  equations  -for  mean,  mean  square  value  and 
variance  {  (J    )     can  be  used  to  solve  -for  the  unkown 
constant,  t  : 


/CO 
t  p(t)  dt   . 
u 


n  square  value  =  | 


mea 


CO 
t"^    p(t)     dt 


2=1  ■■      '^ 


(7^=1        (t     -    yU     )-     p(t)      dt 

o 


=    m< 


fan    square    value    —    U 


(5-2) 


(5-3) 


(5-4) 


Substituting  equation  5-1  into  equations  5—2,  5-3  and  5-4 
leads  to  the  -following  relations: 


°  2~   ' 


mean  square  value  =  t 


(5-5) 


G 


=  [  '^  ^   "   ^o^   ?    ^=  O.l^hjJi 
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This  value  -for  t^  was  used  in  equation  5-1,  and  the 
probability  density  -function  integrated  over  appropriate  . 
limits  to  get  the  number  o-f  events  in  each  20  second  bin. 
The  result  is  plotted  against  the  experimental  distribution 
in  Figure  5—4. 

The  second  distribution  (the  exponential)  belongs  to 
the  family  o-f  gamma  distribution  -functions  C15!]: 


1  a  -t/t^ 

p  =  = t   e        .  (5-6) 

t   O:-^  1  r(Q/  +  1) 


When  Ql     =0,  this  becomes  the  exponential  probability 
density  -function 


1    -t/to 
p  =  —  e        .  (5-7) 


Again  using  equations  5—2,  5—3  and  5—4  leads  to: 

/^  =  t    ;   mean  square  value  =  211  .  (5—8) 

The  exponential  probability  density  -function  was  integrated 
over  the  bins,  and  the  results  are  shown  in  Figure  5-5. 

Another  demonstration  o-f  the  -fit  o-f  the  exponential 
probabilty  distribution  is  shown  in  Figure  5—6.   Taking  the 
natural  log  o-f  the  -function  should  lead  to  a  straight  line 
when  plotted  against  time  or  bin  number.   The  straight  line 
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Figure  5-4   Hal -f  — gaussi  an  distribution  compared  with 
experimental  values. 
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Figure  5-5   Exponential  distribution  compared  with 
experimental  values. 
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Figure  5-6   Semi  log  plot  o-f  exponential  distribution  and 
experimental  values  against  bin  number. 
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in  Figure  5-6  is  a  plot  o-f  the  natural  log  o-f  the  paints 
calculated  using  the  exponential  probability  density 
-function.   The  experimental  points  seem  to  curve  rather 
than  lie  on  a  straight  line. 

The  last  distribution  (J  shaped)  is  also  a  gamma 
distribution.   The  J  shaped  distributions  are  characterized 
by  Oi    <  O  .   I  chose  a  -fairly  common  distribution  with 
Of  =  -0.5  .   The  probability  density  function  is: 


1     -1/2   -t/t„ 
p(t)  =  t      e         ,  (5-9) 


and  the  key  parameters  are: 

'—> 
t   =  2/^    ;   mean  square  value  =  ,  (5— lu) 

n-  o  I — 

This  distribution  is  plotted  against  the  experimental 
values  in  Figure  5—7.   The  natural  log  o-f  both  calculated 
and  experimental  points  are  plotted  against  bin  number  in 
Figure  5—8.   This  distribution  seems  to  -fit  the 
experimental  points  best  o-f  all.   The  J  shaped  probabilty 
density  function  goes  to  in-finity  at  zero,  but  it  is 
i  ntegrabl e. 

A  Chi  square  goodness  o-f  -fit  test  v-jas  done  on  all 
three  distributions.   The  results  a.rEi    summarized  in  Table 
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Figure  5-7   J  shaped  distribution  compared  with 
experimental  values. 
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Figure  5-6   Semi  log  plot  of  J  shaped  distribution  and 
experimental  values  against  bin  number. 
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5-2.   Also  presented  in  Table  5-2  aire    the  ratios  of 
standard  deviation  to  mean. 

Table  5—2   Comparison  o-f  Distribution  Functions 

Chi  square  (J  /  fJi 

Experimental  1.66 

Hal-f-gaussian  127.17  0.76 

Exponential  45.69  1.00 

J  shaped  10.65  1.41 

For  a  distribution  to  pass  a  goodness  o+  -fit  test  it 
must  have  a  Chi  square  less  than  a  prescribed  limit.   The 
limit  for  my  test  (9  degress  of  freedom,  O^  =  0.005)  was 
23.6C161.   Only  the  J  shaped  distribution  passed  the  Chi 
square  test.   It  also  has  (J / LL     closest  to  the 
experimental  values.   In  summary,  the  interarrival  data 
reasonably  fit  a  J  shaped  distribution  given  by: 


^     -1/2  -t/1l± 
p(t)  =  ——-—   t      e         .  (5-11) 


^fW 


Since  event  detection  rate  depended  on  ambient  noise 
level,  interarrival  time  between  events  should  also  show 
environmental  dependence.   Table  5.3  gives  average  and 
standard  deviation  of  the  interarrival  time  for  different 
ambient  noise  pressure  levels. 


mean 

St 

an 

id 

ard  deviation 

(sec) 

(sec) 

67 

129 

190 

153 

147 

143 

1B3 

318 
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Table  5—3   Background  Noise  Level  Dependence  o-f 

Interarrival  Time 


Ambient  Noise  Interarrival  Time 

rms  pressure 
(20-80  Hz) 
(Pa) 

0.01-0.02 

0.02-0.03 

0.03-0.04 

over  0.04 

The  tapes  having  a  background  noise  level  o-f  0.01  to  0.02 
Pa  have  a  significantly  shorter  interarrival  time  than 
tapes  in  the  other  three  pressure  groups.   As  with 
detection  rate,  the    inter  arr  ival     time    does    depend    on 
ambient    noise    leve 1 . 

Spatial  Analysis 

A-fter  removing  nonlocatable  events  and  events  located 
outside  a  horizontal  range  of  20,000  meters,  164  events 
remained.   These  were  grouped  by  horizontal  range  into  42 
annul  i  o-f  equal  area  as  shown  in  Figure  5—9.   Each  annulus 
is  a  30  square  km  ring  centered  at  the  array  origin.   The 
-first  annulus  ("O")  went  -from  O  to  3090  meters,  the  second 
from  3090  to  4370  meters,  and  so  on. 

Table  5—4  shows  the  number  o-f  events  per  annulus  and 
Figure  5—10  shows  this  distribution  graphically. 
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Figure    5-9       Scheme    -for    annul  i    each    o-f    area    equal     to    30    km^. 
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Figure  5-10   Number  o-f  events  per  radius  annulus. 


-94- 


Table  5-4   Number  o-f  Events  per  Annul  us 


n  u  1  u 

0 

1 

2 

3 

4 

5 

6 

7 

B 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

Events 

91 
19 

4 

3 

2 
4 
9 
3 
1 


3 
4 
1 
1 
4 
O 
O 
1 
0 
O 


Annul  us 


21 


24 
25 

26 
27 
28 

29 

3 '-- 
31 


34 
35 
36 
37 
38 
39 
40 
41 


Events 

O 
O 
1 

o 
o 
1 
o 
o 
o 
o 
o 
o 
1 
o 
1 
o 
1 
1 
o 
o 
1 


The  average  number  of  events  per  annul  us  is  3.93  and  the 
standard  deviation  is  14.15  events.   Figure  5-10  shai^s  that 
the    TiLimber    ot    euerits    found    is    highly    deperiderit    on     their 
range    from    the    array.        In  the  center  annulus  there  were 
over  20  times  the  mean  number  o-f  events. 

The  dependence  on  range  is  not  a  surprise,  since 
spreading  (and  possibly  scattering  and  other  losses)  will 
reduce  the  strength  o-f  weak  transients  down  to  the  ambient 
noise  level.   For  this  reason,  the  center  annulus  is 
prcbably  the  best  indicator  o-f  actual  event  density.    In 
this  ring  there  were  91  events  per    30  square  kilometers  per 
662  minutes  o-f  observation  or  aoor oximate Iv    0,3    events    oer 
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square    k i lometer    per    hour. 

The  average  number  o-f  events  per  annul  us  and  the 
number  o-f  events  per  square  kilometer  per  hour  should 
depend  on  background  noise  level.   The  average  number  o-f 
events  per  annuius  was  -found  -for  each  ambient  noise  rms 
pressure  range,  and  adjusted  to  re-flect  the  number  o-f 
events  in  a  662  minute  period.   The  results  ar&    seen  in 
Table  5-5.   The  number  o-f  events  per  square  kilometer  oer 
hour  for  the  center  annuius  Ar^    also  shown  in  Table  5-5. 


Table  5—5   Average  Number  o-f  Events  per  Annul  u; 
-for  4  Ambient  Noise  Levels 


Ambient  Noise    Events 
rms  pressure      per 
(20—80  Hz)       Annulu! 


0.01-0.02  Pa  2.60 

0.02-0.03  Pa  0.29 

0.03-0.04  r-B  0.45 

over  0.04  Pa  0.36 


Minutes 

Adjusted 

#  Events 

o-f  tape 

Events  per 

per  km"^ 

Exami  ned 

Annuius 

per  hr 

287.5 

5.98 

0 .  452 

77 

2.46 

0.260 

100 

2.99 

0.  ISO 

157.5 

1 .  50 

0.  076 

Ent i  re 

Population       3.93         662  3.93        0.27^ 


The    aver  age    namber    of    events    per    aririalus    and    the 
number    of    events    per    square    k i lometer    per    hour    both    reflect 
the    effect    of    signal—to—noise    ratio    on     the    detect  ion 
scheme . 

The  entire  population  of  events  was  investigated  for 
angular  dependence.   Figures  5-11  and  5—12  show  the  number 
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Figure  5-11   Number  o-f  events  per  30°  sector.   Angles  &r> 
measured  -frcm  the  northern  leg  0+  the  hydrophone  array. 
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Figure  5—12   Number  o-f  events  per  SO'-'  sector 


Radi  ui 


gives  the  number  of  events,  while  angle  indicates  the 
sector  measured  -from  the  northern  leg  o-f  the  a.rr3.y.        Each 
ring  represents  10  events. 
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o-f  events  -found  per    30°  sector.   In  the  polar  diagram 
(Figure  5-12)  the  radius  shows  the  number  of  events.   The 
angles  ars    measured  -from  the  northern  leg  o-f  the  airrB.y. 
Figure  5—13  is  a  polar  plot  showing  the  number  ai    events 
per  10°  sector.  There    Mas    no    pr edomiriarit    angular    di r ect ion 
found.       However,  some  preference  can  be  seen  -for  bearings 
o-f  330°  and  190°  -from  the  northern  leg  o-f  the  array. 

Strength  Analysis 

The  mean  hydrophone  peak  pressure  magnitude  -for  each 
event  -fell  within  a  fairly  narrow  band  of  values.   The  mean 
peak  pressures  ranged  from  1.32  to  0.16  Pa,  with  an  average 
of  0.36  Pa  and  a  standard  deviation  of  0.-20  Pa.   Figure  5—14 
shows  the  mean  hydrophone  peak  pressure  values  for  ail 
events  located  between  100  m  and  20,000  m  plotted  against 
range  from  the  array    origin. 

The  different  symbols  shown  in  Figure  5—14  represent 
events  during  each  of  the  four  ambient  pressure  categories. 
The    ei^snts    Nith    a    higher    mean    hydrophone    peak    pr essur e    have 
a  tendency    to    occur    during    higher    amb ient    noise    levels. 
This  can  be  seen  in  Table  5—6,  where  the  maximum,  minimum, 
average  and  standard  deviation  of  the  mean  hydrophone  peak 
pressure  values  ar^    given  for  each  of  the  four  ambient 
noise  levels- 
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Figure  5—14   Mean  hydrophone  peak  pressure  measured  -far 
events  between  100m  and  20,000  m,  plotted  against 
horizontal  range  -from  the  FRAM  IV  array  origin. 
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Table  5-6   Mean  Hydrophone  Peak  Pressure 
for  4  Ambient  Noise  Levels 


Ambient  Noise  Mean  Hydrophone  Peak  Pressure 
rms  pressure                       (Pa) 

(20-BO  Hz)  max  min  average  std  dev 

0.01-0.02  Pa  1.32  0.16  0.31  0.17 

0.02-0.03  Pa  0.56  0.25  0.36  0.09 

0.03-0.04  Pa  0.35  0.23  0.43  0.15 

over  0.04  Pa  1.28  0.16  0.49  0.34 


Enti  re 

PoDulation        1.32         0.16         0.36        0.20 


Source  strength  (F  )  was  -found  -for  the  events  which 
had  hydrophone  locations  between  300  m  and  20,000  m  -from 
the  event. 

The  dipole  strengths  ranged  -from  33  kN  to  4.9  MN, 
with  an  average  o-f  431  kN  and  a  standard  deviation  o-f  555 
kN.   The  distribution  o-f  strengths  -for  the  151  events 
evaluated  is  shown  in  Figure  5—15  and  in  Table  5—7. 

Figure  5—16  shows  the  dipole  strength  -for  all  events 
plotted  against  horizontal  range  -from  the  a.rra.y    origin. 
Again,  it    cari    be    seeri     that    th&    str origer    euerits    occur    Mheri 
the    ambient    pressur e    level     is    high.       Table  5—8  gives  the 
strenath  values  -for  the  di-fferent  ambient  noise  levels. 
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Table  5-7   Strength  Distribution  -for  a  Population  o-f  Events 


O  to   100  kN 

100  to   200  kN 

200  to   300  kN 

300  to   400  kN 

400  to   500  kN 

500  to   600  kN 

600  to   700  kN 

700  to   aOO  kN 

800  to   900  kN 

900  to  1000  kN 

1000  to  1100  kN 

1100  to  1200  kN 

1200  to  1300  kN 

1300  to  1400  kN 

1400  to  1500  kN 


1.5 
over 


to 


MN 
MN 


#  o-f  Events 

19 
29 
33 
28 

8 

5 

4 

5 

2 

5 

3 

2 

3 

1 

1 

1 

2 


Table  5—8   Dipole  Strength  versus  Ambient  Noise  Levels 


Ambient  Noise 
rms  pressure 
(20-80  Hz) 


0.01-0.02  Pa 
0.02-O.03  Pa 
0.03-0.04  Pa 
over  0.04  Pa 


Dipole  Strength  (F^) 
(kN) 


max 
1051 
2041 
1153 

4939 


n 

average 

std  dev 

33 

259 

155 

64 

649 

634 

79 

643 

359 

59 

860 

1449 

Enti  re 

Popui  ati on 


4939 


431 


555 


-103- 


0) 

c 

> 

Id 


(P 

E 
z 


34 

32- 

30- 

28- 

26- 

24 

22- 

20- 

18-. 

16-' 

14- 


12- 
10-/. 

8- 

6- 

4- 

2- 

0 


V, 


/. 


^/ 


V 


/. 


V 


mv 


0-100         300-400        600-700       900-1000     1200-1300 

Event  Strength  (kN) 


Figure  5-15   Distribution  o+  strength  tor  a  population  o+ 
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Figure  5-16  Dipole  strength  -For  events  between  300  m  and 
20,000  m,  plotted  against  horizontal  range  -from  the  array 
or  i  g  i  n . 
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iiost  of  the  events  evaluated  -for  strength  occurred 
during  the  lowest  ambient  noise  levels.   Figure  5—17  shows 
the  strength  of  events  that  occurred  when  the  ambient  noise 
was  0.01  to  0.02  Pa.   The  log  o-f  the  dipole  strength  is 
plotted  against  the  log  o-f  the  horizontal  range  -from  the 
center  o-f  the  array.   The  points  scatter  more  so  to  the 
upper  le-ft  rather  than  lower  right,  because  distance  itsel-f 
•filters  out  weak  events.   A  weak  signal  -from  -far  away  would 
not  reach  the  hydrophone  array  with  enough  amplitude  to  be 
distinguished  -from  the  background  noise.   And  events 
located  farther  away  would  tend  to  be  strong  events. 
However,  events  located  close  to  the  Array    should  have  the 
entire  range  o-f  source  strength  levels.   This  would  produce 
a  wedged  shaped  plot  o-f  weaker  events  close  to  the  Array. 
Indeed,  Figure  5—17  shows  a  general  scattering  with  perhaps 
a  wedge  o-f  weaker  events  near  the  array    oriqirs. 

Nonetheless  the  trend  shown  in  Figure  5—17  suggests 
that  the  ray  average  model  used  to  estimate  re-fractive— 
sur-face  reflective  spreading  may  need  to  be  replaced  with  a 
more  refined  model.   For  example,  horizontal  ranges  less 
than  about  1000  m  may  include  too  small  a  loss,  and 
therefore  lead  to  too  small  a  strength,  because  the 
reflective  contributions  may  not  be  as  large  as  imputed. 
Such  a  cri'iicism  is  supported  by  the  notion  that  for  a 
given  slope,  £  ,  reflective  rays  and  hence  ray  averaging 
occurs  onlv  bevond  a  critical  horizontal  ranae. 


-106- 


z 


L 

CD 
T3 

£ 

c 
(P 

L 

\l) 


0 

a 

Q 

0 


0 


61  -r 

60- 
59- 
58- 
57- 
56- 
55- 
54- 
53- 
52- 
51  - 
50- 
49- 
48- 
47- 
46- 
45- 


+    + 


+ 


+ 


+    + 


+ 


.  t 


24      26       28      30       32      34      35       38       40       42       44 
10  log  Horizontal  Range  (dB  re  1  m) 


Figure  5-17   10  log  o-f  Dipole  strength  (dB  re  1  N)  versus 
10  log  o-f  horizontal  range  -from  the  Array    origin  (dB  re  1  m) 
-for  events  occurring  during  an  ambient  noise  level  o-f  O  ni 
to  0.02  Pa. 
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The  foregoing  speculation  suggests  that  the  average 
dipole  strength  -for  the  lowest  ambient  noise  case  is  best 
■found  -from  the  events  farther  from  the  array,  and  is 

Fq  ^  10^-^  N  ^  320  kN   ,  (5-12) 

with  a  much  smaller  standard  deviation  than  in  Table  5—8. 
Presumably  corresponding  adjustments  could  be  made  for  the 
higher  ambient  noise  cases,  but  the  FRAM  IV  data  set 
contains  too  few  events  at  higher  ambient  noise  to  plot  as 
in  Figure  5—17. 

The  strength  analysis  is  a  somewhat  ambivalent  one 
because  of  spreading  model  uncertainty,  and  because  data  on 
ice  slopes  are  not  available-   But  the  dipole  picture  of  an. 
event  likely  has  some  validity,  and  at  least  rough 
estimates  of  its  strength  have  been  extracted  from  the 
data. 
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CHAPTER  6 
SUMMARY  AND  THOUGHTS 

Through  the  use  ot  a  detection  program,  visual 
con-f  i  rmati  on  and  a  location  program,  a  population  of  199 
Arctic  noise  transients  was  gathered.   There  Are    -four  major 
resLil  ts. 

First,  more  events  ArB    found  when  the  afribiant  pr^ss.ire 
is  low,  and  more  -false  alarms  when  the  ambient  pressure  is 
high.   The  inter  arrival  time  and  the  average  number  o-f 
events  per  unit  area  also  depend  on  ambient  noise  level. 
Since  more  events  a.rB    -found  when  the  ambient  noise  is  low, 
the  interarrival  time  decreases,  and  the  spatial  density 
i  ncreases. 

Second,  the  interarrival  times  were  -fit  to  several 
possible  Drobability  distributions.   The  interarrival  tiiTe 
distribution  best  -fits  a  J  shaped  gamma  distribution.   The 
mean  interarrival  time  is  100  seconds. 

Third,  the  number  o-f  events  per  unit  3.r&a.    is  highly 
dependent  on  range,  since  distance  -filters  out  weak 
transients.   The  event  density  in  the  annul  us  closest  to 
the  center  o-f  the  Arrsiy    was  0.3  events  per  square  kilometer 
per  hour  over  all  observations  and  0.5  events  per  square 
kilometer  per  hour  -for  cuiet  times.   There  is  no 
predominant  angular  dependence  to  the  spatial  distribution 
o-f  events. 


-109- 

Last ,  the  mean  dipole  strength  for  the  observecS  events 
is  430  kN  overall  and  260  kN  during  low  ambient  noi^e 
levels.   Stronger  events  occurred  during  high  ambient  noise 
levels.   A  re-finement  o-f  the  spreading  loss  model  used  to 
calculate  these  values  may  lead  to  values  which  a.re 
slightly  higher. 

Anal  ysi  s  o-f  Arctic  acoustic  events  is  far  from 
complete.   Several  areas  for  improvement  have  been 
mentioned  earlier  in  the  thesis.   The  detection  program 
needs  to  be  made  more  robust  to  eliminate  the  event  time 
error.   A  scheme  for  ignoring  artifacts  should  be  included. 
The  location  program  wastes  time  looking  in  the  wrong 
direction,  although  the  bearing  accuracy  of  the  program  is 
"■y^ry    good.   The  algorithm  should  be  changed  to  quickly  +ind 
the  right  bearing,  and  then  search  in  a  sector. 

The  type  of  each  event,  whether  it  was  a  pop  or  a 
whine,  >jas  not  recorded.   Collecting  this  infor'mation  and 
correlating  it  with  interarrival  time  and  range  still  needs 
to  be  done. 
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APPENDIX  A 


User's  Guide  -for  the  hdetect    Program 


Figure  A-1:   Flow  Chart  of  the  hdetect    Program 


Source  Coda  -For  the  hdetect    Program 
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USER'S  GUIDE  FOR  THE  hdetect    PROGRAM 

The  purpose  o-F  the  hdetect    program  is  to  detect 
ambient  noise  transients  amidst  the  background  ambient 
noise  recorded  on  a  FRAM  data  tape.   This  is  done  by 
comparing  the  short  average  o-f  data  points  to  the  long 
average  o-f  points  on  a  single  channel  in  order  to  flag  a 
possible  detection,  and  then  waiting  until  50X  o-f  the 
channels  a.r&    -flagged  to  declare  an  actual  detection. 

The  input  -for  the  hdetect    program  is  a  framread    output 
-file  without  headers.   A  FRAM  data  tape  is  read  into  the 
-file  by  the  command 

-•framread  -head       < RETURN > 
The  program  will  ask  -for  the  input  device  (tape  drive 
designation)  ,  the  output  -file,  and  the  number  o-f  data 
segments  to  skip  and  to  read.   Each  segment  represents  3.8 
seconds  o-f  data  on  24  channels.   The  framread    program  reads 
a  first  segment  which  contains  no  data  records,  so  you 
should  specify  skipping  one  more  segment  than  you  would 
normally  calculate.   For  example,  reading  the  entire  first 
half  of  a  20  minute  FRAM  IV  tape  would  reguire  the  response 
of 

1  160  <  return::: 

to  the  question  of  ''enter   #skip,  #seqments:". 

Once  this  input  file  has  been  created  the  hdetect 
program  can  be  used.   The  program  is  started  with  the 
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command 

hdetect  <RETURN> 

The  program  will  ask  for  the  FRAM  tape  number,  the  Julian 

date  o-f  the  tape,  and  the  start  time  o-f  the  tape  in  hours, 

minutes  and  seconds.   The  program  will  then  ask  you  to 

select  the  channels  you  wish  to  use.   In  most  cases  the 

FRAM  data  tapes  did  not  have  ambient  noise  hydrophones  tied 

into  all  channels,  and  the  speci-fic  channel  that  a 

hydrophone  was  recorded  on  changed  throughout  the 

experiment.   Which  channels  were  in  use  and  for  which 

hydrophones  can  be  found  in  the  experiment  logs.   The 

program  assumes  that  the  channel  number  is  equal  to  the 

hydrophone  number,  but  allows  you  to  change  t-his  by 

inputting  the  channel  number  and  the  proper  hydrophone 

number,  or  "O"  if  the  channel  is  not  in  use.   For  example, 

if  channel  3  was  not  used,  and  channel  7  was  used  for 

hydrophone  21,  the  input  would  be 

3,0  < RETURN > 

7,21  < RETURN > 

0,0  < RETURN > 

The  "0,0"  ends  the  changes  to  the  channel  selection.   You 

must  now  hit  any  key  to  continue  the  program. 

You  will  be  asked  to  enter  the  input  device  (the  input 

framread    file) ,  the  number  of  skips  and  segments,  and  the 

name  of  the  output  file.   The  output  file  does  not  have  to 

exist  before  the  program  is  started.   It  will  be  created  by 

the  program.   The  number  of  skips  and  segments  eire    those 
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that  you  would  calculate  using  3.8  seconds  per  segment. 
For  ex-ample,  to  proces  the  entire  first  hal-f  of  a  20  minute 
FRAM  IV  tape  the  number  of  skips  and  segments  would  be 
0,160  < RETURN > 

This  is  all  of  the  input  required  by  the  user.   The 
program  proceeds  from  this  point  without  user  interaction. 

The  output  of  the  hdetect    program  is  a  file  containing 

a  list  of  detections  in  the  following  format: 

tapenumber  Juliandate  hour  minute  seconds 
0  eventnumber  eventtime 

channel  hydrophone  timedelay  amplitude 
channel  hydrophone  timedelay  amplitude 
channel  hydrophone  timedelay  amplitude 


channel  hydrophone  timedelay  amplitude 

O  eventnumber  eventtime 

channel  hydrophone  timedelay  amplitude 


channel  hydrophone  timedelay  amplitude 
O  eventnumber  eventtime 


-1 
The  "O"  at  the  start  of  a  line  indicates  a  new  event 
detection,  and  the  "—1"  at  the  start  of  a  line  indicates  an 
end  of  file.   Each  channel  that  was  flagged  for  a 
particular  event  is  listed  with  its  hydrophone  number, 
timedelay  from  the  earliest  channel  signal  arrival,  and  its 
peak  voltage  amplitude.   This  outfile  can  be  used  as  the 
input  file  for  the  location  programs  without  modification. 
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A-1   Flow  chart  o-f  the  hdetect  program. 
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Figure  A-1   Flow  chart  o-f  the  hdetect  program, 
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Figure  A-1   Flow  chart  o-f  the  hdetect  program. 
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/*  startdoc 
hdetect.c 

program  to  read  whoi-segy  format  data  tapes  from  the  fram  IV  program. 

usage : 

hdetect 

program  is  interactive. 


by  Mary  Townsend-Manning 

enddoc 

*/ 


♦include  <3tdio.h> 

♦include  <math.h> 

♦define  NCHAN  2  5 

♦define  RECLN  950     /*  number  of  samples  per  trace  (4  bytes  per  sample)*/ 

♦define  OBYTES  3800   /*  number  of  bytes  per  record  output  */ 

♦define  ZERO  0 

♦define  LONGFILTLN  64  /*  length  of  long  average  filter   */ 

♦define  FLN  64 

♦define  MAXEVENTS  4 

♦define  RESET_DELAY  0.3 

♦define  RATIO  2.38 

♦define  SIG_OELAY  0.02 

♦define  THRESHOLD  0.5 

♦define  EVENT_DELAY  0.5 

♦define  END  -1 

double  vconv(x,y,n) 

register  float  *x,  *y; 
register  int  n; 

I 

register  double  sum  =  0 . ; 

if(n  >  0) 
( 

do   ( 

sum  +«  *x++  *  *y — ; 
)  while ( — n  >  0) ; 
} 
return (sum) ; 
1 


mainO 
C 


Source  code  -for  the  hdetect    program. 


1   "P*^ 


float  l_ave [NCHAN] ,  sh_ave,  timeof f lag [NCHAN] ; 

float  ampofflag [NCHAN] ,  firsttime,  t imede lay [NCHAN] ; 

float  chaneventtiine[MAXEVENTS]  [NCHAN]  ; 

float  chaneventamp[MAXEVENTS] [NCHAN] ,  event ime [MAXEVENTS] ; 

int  flag [NCHAN] ; 

int  miin_active_events; 

int  m,  n; 

int   event_flag[MAXEVENTS] [NCHAN] ,    flag_sum,    nuinber_of_events; 

int   event_niiniber  [MAXEVENTS  ]  ; 

int    nchan   -   0 ; 

int  Chan [NCHAN],  k,  j,  i,  tid,  date,  hour,  min,  sec,  1; 

int  channel,  data,  toggle; 

char  answer; 

float  h [FLN] ,  longf ilt [LONGFILTLN] ; 

int    nskip,nseg  ; 

char  oddobuf [NCHAN] [OBYTES] , evenobuf [NCHAN] [OBYTES]; 

int  count, error; 

char  f name [80] , iname [80] ; 

FILE     *iptr,  *ptr,  *fp,  *fopen(); 

float    tiineseries  [NCHAN]  [RECLN+2*FLN]  ; 

float   time   =   0  . ; 

num_active_event3  -  0; 

for  (i-l;i<25;i++)  { 
Chan [i]- (i) ; 
flag[i]-  0  ; 

for (m-l;m<MAXEVENTS;m++)   { 
event_flag[m] [i]  -  0; 
chaneventtime [m] [i]  =  0; 
chaneventamp [m] [i]  =  0; 
) 
) 

/*  Program  initialization  from  keyboard  */ 


fprintf (stderr, "Program  InitializationXn") ; 
fprintf (stderr, "enter  FRAM  tape  #\n") ; 
fscanf (stdin, "%d", &tid) ; 

fprintf (stderr, "enter  Julian  dateSn") ; 
fscanf (stdin, "%d", Sdate) ; 

fprintf (stderr, "enter  time  -  HR, MN, SC\n") ; 
fscanf (stdin, "%d, %d, %d", &hour, Smin, Ssec) ; 

fprintf (stderr, "default  values  for  channels  and  phonesNn") ; 
fprintf (stderr, "are  channel  #  =  phone  #.\n"); 

fprintf (stderr, "enter  channel,  phone  to  change. \n"); 

fprintf (stderr, "enter  '0'  for  phone,  to  eliminate  a  channel. \n") 

fprintf (stderr, "enter  '0,0'  to  quit.Xn"); 

fscanf (stdin, "%d, %d", S j,Sk) ; 


while    (j    !-  0   &s    j   <   25)     ( 
chan[j]    -   k; 


Source    code    -For    the    hdetect    program. 
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fscanf  (stdin,  "%d,  %d",  S  j,  &lc)  ; 
) 

fprintf (atderr, "FRAM  TAPE  %d    Julian  Date:  %d\n", tid, date) ; 

fprintf (stderr, "  Time:        %d: %d: %d\n",  hour, min, sec) 

for{i=l;i<13;i++) 

fprintf (stderr, "CH   %d  PH   %d       CH   %d  PH   %d\n", 
i,chan [i] , i+12, chan [i+12] ) ; 


/*  Checlc  to  make  sure  inputs  are  correct  —  Change  if  necessary  */ 


f scanf (stdin, "%c", &answer) ; 

fprintf  (stderr,  "Hit  any  Icey  and  RETURN,  when  ready.") 

f scanf (stdin, "%c",  Sanswer)  ; 


for(i=l;i<NCHAN;i++)  { 

if  (Chan [i]  !-  0)  nchan++; 
) 

fprintf (stderr, "enter  input  device:   ")  ; 
scanf ("%s", iname) ; 

if((iptr  -  fopen (iname,  "r") )  ~  NULL) 
( 

fprintf (stderr, "can' t  open  %3\n", iname) ; 
exit ( 1 ) ; 
1 


fprintf  (stderr,  "enter  #s)cip,  #segments  :  \n")  ; 

fprintf (stderr, "values  of  0  and  320  will  read  entire  tape\n") 

scanf  (''%d,  %d",  Snskip,  finseg)  ; 

/•  load  bandwidth  filter  */ 

if ( (fp-fopen ("PMfloat", "r") )  ==  NULL)  ( 

printf  ("cannot  open  bandwidth  filter  file\n") ; 

exit  ( 0 )  ; 

} 
for (i=0;i<64;i++) 

fscanf (fp,"%f",ih[i] )  ; 
f close (fp) ; 


/*  load  averaging  filter  */ 


f  or  (i-0;  KLONGFILTLN;  i++) 

longfiltCi]  -  1 . 0/ (float) LONGFILTLN; 


/•  Open  output  file  */ 
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fprintf (stderr, "enter  out-file:   ") ; 

scanf ("%3"/ fname) ; 

if  (  (ptr   -   fopen(fnaine,    "w"))    »-  NULL) 

{ 

fprintf (stderr, "can' t  open  %s\n", fname) ; 

exit (1) ; 

} 

fprintf (ptr, "%d  %d  %d  %d  %d\n",  tid,  date,  hour,  min,  sec) ; 

time  »  time  +3.8  *  (nskip-1)  ; 

/*  enter  first  record   */ 

if  (ns)cip%2  ==  1)  ( 

for ( j-1; j<NCHAN; j++)  { 

fread(&evenobuf [ j] [0] ,  sizeof (float) ,  RECLN,  iptr) ; 
} 

toggle  »  1; 
} 

else  { 
for{j=l; j<NCHAN; j++)  { 

f read(&oddobuf [ j] [0] ,  sizeof (float) ,  RECLN,  iptr) ; 
} 

toggle  -  0; 
) 

/*  ENTERING  RECORD  READING  MODULE  */ 
time  -  time  +0.50  4; 

/•      fprintf (stderr, "using  buffer  size  %d  bytes\n", sizeof (buf)  ) ;  */ 

for(i=l;  i  <  nseg;  i++) 

{ 
/*        fprintf (stderr, "time  =  %f\n",  time) ; 

fprintf  (stderr,  "processing  record  %d\n",  ns)cip+i)  ;  */ 

/*  read  next  record  into  appropriate  buffer  */ 

if  (toggle  ==»!){ 

for ( j=l; j<NCHAN; j++)  { 

fread(&oddobuf [ j] [0] ,  sizeof (float) ,  RECLN,  iptr) ; 
} 

toggle  =  0; 
) 

else  ( 

for ( j-1; j<NCHAN; j++)  { 

fread(Sevenobuf [ j]  [0] ,  sizeof (float)  ,  RECLN,  iptr); 

) 

toggle  -  1 ; 
) 

/*    filter  and  square  data    */ 
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if  (toggle  =-  1)  { 

for( j=l; j<NCHAN; j++)  { 
if(chan[j]  !-  0)   { 

3q_filt (Soddobuf [ j] [0] , Sevenobuf [ j ] [0] ,h, 
Stimeseries [ j] [0] ) ; 

l_ave[j]  -  vconv(Slongfilt [0] ,&timeseries[ j] [2* (FLN-1) ] , 
LONGFILTLN) ; 

) 

) 
else  { 

for( j=l; j<NCHAN; j++)  { 

if(chan[j]  !=  0)   { 
sq_filt (Sevenobuf [ j] [0] ,fioddobuf [ j] [0] ,h, 
Stimeseries [ j] [0] ) ; 

l_ave[j]  =■  vconv(slongfilt  [0]  ,  Stimeseries  [j]  [2-«  (FLN-1)  ]  , 
LONGFILTLN) ; 
} 
} 
) 


/*  ENTERING  EVENT  DETECTION  MODULE   */ 

for(k=0;)c<R£CLN;)c   +=4)     { 

for(l=l;l<NCHAN;l++)     ( 
if    (chan[l]     !=   0)     { 

l_ave[l]    =    (63.0*l_ave[l]    +   timeseries  [1]  (2*  (FLN-1) +lc]  ) /64  .  0; 
sh_ave    =    (timeseries  [1]  [2»  (FLN-1) +ic]    + 
timeseries  [1]  [2*  (FLN-1) +k-l]    +   timeseries  [1]  [2*  (FLN-1) +k.-2]    + 
timeseries  [1]  [2*  (FLN-1) +l<-3]  )  /4.0; 


/*  reset  old  flags  */ 


if  (flag[l]  ==  1  SS  (time  -  timeof f lag [1] )  >  R£SET_DELAY) 
flag[l]  =  0; 


/*  set  flag  if  RATIO  of  signals  is  reached  */ 


if  ( (sh_ave/l_ave[l] )>=RATIO)  { 
if  (flagCl]  —  1)  { 

if  (sh_ave  >  ampof f lag[l] )  ( 
timeof flag [1]  =  time; 
ampofflag[l]  -  3h_ave; 

) 
} 

else    ( 

if    (n\im_active_events    =»=    0)     { 
flag[l]    =   1; 
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timeoff lag[l] 
ainpofflag[l]    = 


■  time; 
sh  ave; 


( 


) 

else  { 

for  (m=l;m<  (num_active_events  +  1)  ;ni++) 
if  (event_flag[in]  [1]  ~  1)  { 

if  ((time  -  chaneventtime [m] [1] ) 
<-  SIG_DELlAY)  { 
if  (sh_ave  >  chaneventamp[m] [1] ) { 
chaneventtime [m] [1]  -  time; 
chaneventamp[m] [1]  -  sh_ave; 
} 
) 
elae  ( 

flag[l]  =-  1; 
timeof f lag[l]  =  time; 
ampofflag[l]  =  sh_ave; 
} 
} 
else  { 


event_f  lag(m]  [1]  =  1; 
chaneventtime [m] [1]  =  time; 
chaneventamp [m] [1]  =  sh_ave; 
m  "  num  active  events; 


} 


) 


/*  end  of  set  flag  module  */ 
/*  start  new  event  module  */ 
flag_sum  =  0; 


for  (1=1;1<NCHAN;1++)  { 

if  (flagCl]  =-  1)  flag_3um++; 
) 
if  (( (float) flag_3um/ (float) nchan)  >=  THRESHOLD)  { 

num_active_events++; 

eventime  [num_active_events]  =»  time; 

number_of_events++; 

event_number [num_active_events]  =  number_of_event3; 

for  (1-1;1<NCHAN;1++)  { 

event_f lag [num_active_events] [1]  =  flag[l]; 
chaneventtime [num_active_event3] [1]  "  timeof flag [1] 
chaneventamp[num_active_events] [1]  -  arapof f lag[l] ; 
flag[l]  =■  0; 
timeofflag[l]  «  0; 
aunpof  f  lag  [1]  =  0; 
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} 
} 

/*  end  of  new  event  module   */ 

/*  start  of  deactivate  old  event  module  */ 

if  (num_active_event3  >  0  S4  (time  -  eventime[l])  >  EVENT_DELAY)  { 
fprintf  (ptr,"%d  %d   %f\n",ZERO, 

event_n umber [1] , event ime [1]  )  ; 
/*    fprintf  (stderr,"%d   %f \n",event_number [1] ,eventime [1] ) ;   */ 


find  time  delays  by  finding  earliest  channel  event 
time,  and  subtracting  that  from  the  other  channel 
times 


firsttime  »  100  00.0; 

for (1=1;1<NCHAN;1++)  { 

if(chan[l]  !=■  0  &&   event_f  lag  [IHl]  !=  0  &S 
chaneventtime [1] [1]  <  firsttime) 
firsttime  =  chaneventtime [1] [1] ; 
) 

for(l-l;l<NCHAN;l++)   ( 

if(chan[l]  !=  0  &S  event_f lag [1] [1]  !=  0)  ( 

timedelay[l]  =  (  (chaneventtime [1]  [1] )  -  firsttime) 
fprintf (ptr, "%d   %d   %f   %f \n", 1, chan [1] , 
timedelay [1] ,  chaneventamp [1] [1]); 


} 


1 


print  to  file  to  indicate  end  of  event    */ 

for  (1=1; l<num_active_events; 1++)  ( 
for  (m=l;  m<NCHAN;m++)  ( 

event_f lag [1] [m]  =  event_f lag(l+l] [m] ; 
chaneventtime [1] [m]  =  chaneventtime [1+1] [m] 
chaneventamp [1] [m]  =  chaneventamp [1+1] [m] ; 
} 

eventiine[l]  =»  eventime  [1+1]  ; 
event_number [1]  =  event_number [1+1] ; 
} 

num_active_event3 — ; 
} 
time  +-  0.004*4.0; 
) 


/*  EXIT  MODULE 
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/*   fprintf(stderr, "processing  record  %d\n",  nskip+nseg) ;  */ 
if (toggle  ==1)   { 

for(j=l; j<NCHAN; j++)   { 
if(chan[j]  !=  0)   I 
f  or  (k'O  ;  )c<OBYTES ;  k++) 

oddobuf [j] [k]  =  0; 
sci_filt (Sevenobuf [ j] [0] , Soddobuf [ j] [0] ,h, Stimeseries [ j] [0] ) ; 

l_ave[j]  -  vconv(&longfilt [0] ,&tiroeseries [ j] [2* (FLN-1) ] , 
LONGFILTLN) ; 
) 
} 
) 
else   ( 
for ( j=l; j<NCHAN; j++)   { 
if(chan[j]  !=  0)  { 
for (k=0 ; k<OBYTES ; k++) 
evenobuf [ j] [k]  =  0; 
3q_f ilt (Soddobuf [ j] [0] ,&evenobuf [ j] [0] ,h, &timeseries [ j] [0] ) ; 

l_ave[j]  =  vconv(Slongfilt [0] , Stimeseries [j] [2* (FLN-1) ] , 
LONGFILTLN) ; 
1 
1 
} 


for  (k='0;k<  (RECLN-2*  (FLN-1)  )  ;k+=4)   { 

for (1=1;1<NCHAN;1++)  ( 

if  (chanCl]  !  =■  0)  { 

l_ave[l]  -  (63.0*l_ave[l]  +  timeseries [1] [2* (FLN-1) +k] ) /64 . 0; 
3h_ave  -  (timeseries [1] [2* (FLN-1) +k]  + 
timeseries  [1H2*  (FLN-1) +k-l]  +  timeseries  [1]  [2*  (FLN-1) +k-2]  + 
timeseries [1]  (2*  (FLN-1) +k-3] ) /4.0; 


/*  reset  old  flags  */ 


if  (flag[l]  ==  1  &&  (time  -  timeof f lag [1] )  >  RESET_DELAY) 
flagtl]  =  0; 


/*  set  flag  if  RATIO  of  signals  is  reached  */ 


if  ( (sh_ave/l_ave[l] )>=RATIO)  ( 
if  {flag[l]  -=  1)  ( 

if  (sh_ave  >  ampof f lag(l] ) 
timeof  f  lag  [1]  =•  time; 
ampofflag[l]  =  sh  ave; 


) 


] 


else  { 

if  (num_active_events  ==0)  { 
flag[l]  =  1; 
timeof f lag [1]  =  time; 
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ampoff lag[l]  =  sh_ave; 
} 
else  { 

for  (m=l;m< (niun_active_event3  +  1) ;m++)  { 
if  (event_flag[m)  [1]  =■=  1)  { 

if  ((time  -  chaneventtime [m] [!}) 
<-  SIG_DELAY)  ( 
if  (3h_ave  >  chaneventamp [m] [1] ) { 
chaneventtime [m] [1]  =  time; 
chaneventamp [m] [1]  "  3h_ave; 
1 
} 
else  { 

flag[l]  =  1; 
timeoff lag[l]  =  time; 
ampofflag[l]  =  sh_ave; 
) 
■  ) 
else  ( 


event_flag[m] [1]  =  1; 
chaneventtime [m] [1]  =  time; 
chaneventamp [m] [1]  =  sh_ave; 
m  »  nura  active  events; 


) 
} 

/*  end  of  set  flag  module  */ 

/*  start  new  event  module  */ 

flag_sum  =  0; 


for  (1=1;1<NCHAN;1++)  ( 

if  (flag[l]  ==  1)  flag_sum++; 
) 
if  (( (float) flag_sum/ (float) nchan)  >=  THRESHOLD)  { 

num_active_event3++; 

eventime [num_active_events]  =  time; 

numbe  r_o  f  _e  vent  s ++ ; 

event_number [nuin_active_event3]    =  numbe r_of_e vents; 

for  (1-1;1<NCHAN;1++)  ( 

event_f  lag [num_active_events]  [1]  =  flag[l]; 

chaneventtime  [num_active_events]  [1]  =  timeoff lag [1] 

chaneventamp [num_active_event3] [1]  =  ampof f lag[l] ; 

flag[l]  -  0; 

timeofflagtl]  -  0; 

ampofflag[l]  =0; 
) 
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1 

/*  end  of  new  event  module   */ 

/*  start  of  deactivate  old  event  module  */ 

if  (nuin_active_events  >  0  SS  (time  -  eventime[l])  >  EVENT_DELAy)  { 
fprintf  (ptr,"%d  %d   %f\n",  ZERO, 

event_nuinber  [  1  ]  /  eventime  [  1  ] )  ; 
/*    fprintf  (stderr,"%d   %f  \n"/event_nuinber  [1] , eventime  [1] )  ;  */ 


find  time  delays  by  finding  earliest  channel  event 
time,  and  subtracting  that  from  the  other  channel 
times 


*/ 


firsttime  -  10000.0; 

for(l=-l;l<NCHAN;l++)  { 

if(chan[l]  !-  0  &&   event_f lag [1] [1]  !-  0  && 
chaneventtime [1] [1]  <  firsttime) 
firsttime  "  chaneventtime [1] [1] ; 
1 

for(l=-l;l<NCHAN;l++)   { 

ifCchanfl]  !=  0  &&   event_f lag[l] [1]  !=  0)  { 

timedelayCl]  »  (  (chaneventtime [1]  [1] )  -  firsttime) 
fprintf (ptr, "%d   %d  %f   %f \n" , 1, chan [1] , 
timedelayCl],  chaneventamp [1] [1]); 
) 
) 


print  to  file  to  indicate  end  of  event    */ 

for  (1=1; l<num_active_events; 1++)  ( 
for  (m=l;  m<NCHAN;m++)  { 

event_flag[l] [m]  -  event_f lag [1+1] [m] ; 
chaneventtime [1] [m]  =  chaneventtime [1+1] [m] ; 
chaneventamp [1] [m]  -  chaneventamp [1+1] [m] ; 
} 

eventime [1]  -  eventime [1+1] ; 
event_nximber  [1]  -  event_number  [1+1]  ; 
) 

nuin_active_event3  — ; 
} 

time  +-  0.004*4.0; 


1 
/*   print   out   all   events    */ 

if (num  active  events   >   0)     [ 
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f  or  (k-l;lc<  (num_active_event3  +  1), •)«++)   { 

fprintf  (ptr,"%d  %d   %f\n'',  ZERO, 

event_n umber  [Jc]  ,  event ime  [k] )  ; 
/*       fprintf  (stderr,"%d   %f  \n",event_nuinber  [k]  ,  eventime  [k]  )  ;  */ 


find  time  delays  by  finding  earliest  channel  event 
time,  and  subtracting  that  from  the  other  channel 
times 


*/ 


firsttime  =  10000.0; 

for (1-1;1<NCHAN;1++)  { 

if(chan[l]  !-  0  £&  event_f lag [k] [1]  !=  0  && 
chaneventtime [k] [1]  <  firsttime) 

firsttime  "  chaneventtime [k] [1] ; 
} 

for (1-1;1<NCHAN;1++)   { 

if(chan(l]  !-  0  &&   event_f  lag  [k]  [1]  !=•  0)  ( 

timedelay[l]  »  ( (chaneventtime [k] [1] )  -  firsttime) 
fprintf (ptr, "%d   %d   %f   %f \n", 1, chan [1] , 
timedelay [1] ,  chaneventamp [k] [1]); 
) 
) 


) 


) 

/*  final  summary  to  screen  */ 


fprintf (ptr, "%d",  END) 
fclose(ptr)  ; 
f close (iptr)  ; 
exit(O)  ; 


sq_filt (first, second, filter, output) 

/*  filters  and  squares  two  data  arrays  */ 


float  first [RECLN] ,  second [RECLN] ,  filter [FLN] 
float  output [RECLN+2*FLN] ; 
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/*  put  zeros  in  first  (FLN  -1)  places  of  output  */ 

int  j,  i; 

float  transition [2*FLN] ; 

float  sum; 

for (i-0 ; i<FLN-l ; i++) 

output [i]  =  0.0  ; 

/*  put  in  first  data   */ 

for (i=FLN-l;i<RECLN;i++)  { 

sum  =  vconv(&filter [0] , &first [i] ,FLN) ; 
output  [i]  =  sqrt  (sxom  *  sum)  ; 
} 

/*  put  in  transition  from  first  to  second  data  */ 

for(j=0; j<FLN-l; j++)  { 

transitionC j]  -  f irst [RECLN- (FLN-1)  + j] ; 

transition [ j+FLN-1]  =  second[j]; 
) 

for(i-0;i<FLN-l;i++) { 

sum  «  vconv(Sfilter [0] ,Stransition[i+FLN-l] ,FLN) 
output [i+R£CLN]  =  sqrt{sum  *  sum)  ; 
) 

/*  put  in  second  data  */ 

for (i=FLN-l;i<2*FLN;i++)  { 

sum  -  vconv(Sf ilter [0] , Ssecond[i] , FLN) ; 

output  [RECLN+i]  =»  sqrt  (sum  *  sum)  ; 
) 
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USER'S  GUIDE  FOR  THE  location,    far  locate, 
and  finelocate  PROGRAMS 

The  purpose  o-f  the  location  programs  is  to  find  the 
spatial  location  o-f  an  event  -from  the  time  delays  between 
signal  arrival  at  different  hydrophones.   The  program 
assumes  a  test  location  and  computes  the  slant  range  to  the 
individual  hydrophones.   The  slant  ranges  are  plotted 
against  the  experimental  time  delays  and  a  least  squares 
fit  is  done.   The  test  location  with  the  best  least  squares 
fit  is  considered  the  location  of  the  event. 

The  input  to  the  location  programs  is  the  output  file 
of  the  detection  program.   Manual  time  delays  may  be 
substituted  -for  the  program  generated  time  delays  in  this 
file,  but  this  editting  must  be  done  before  the  location 
program  is  invoked. 

This  location  program  is  very  user  interactive.   The 
user  starts  the  location  program  with  the  command 

location  < RETURN > 

The  program  asks  for  the  input  and  output  file  names.   It 
then  reads  the  input  file  and  asks  whether  the  user  would 
like  to  locate  the  first  event.   This  allows  the  user  to 
skip  down  to  the  event  of  interest.   The  program  then 
allows  the  user  to  adjust  which  hydrophone  time  delays  will 
be  used  in  the  location  process.   This  is  very  handy  for 
removing  questionable  time  delays,  in  order  to  get  a  better 
location  -solution.   With  the  hvdroohone  channels  chosen. 
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the  pragram  proceeds  with  the  actual  location  algorithm. 
In  the  t  jn^-Icoate    and  far  locate    programs  the  user  is  asked 
to  speci-fy  which  quadrant  or    direction  is  to  be  searched. 

The  I ocat ion    program  tries  test  locations  in  a  large 
grid,  and  when  the  "best"  location  is  -found,  then  searches 
a  smaller  grid  around  this  "best"  location.   The  locat ior: 
and  tar  I  ocate    programs  have  4  levels  o-f  grids  and  the 
f iTi-e locate    program  has  3  levels.   Intermediate  answers  ar^ 
displayed  -for  each  level. 

The  intermediate  and  -final  answers  display  the  x  and  y 
coordinates  o-f  the  best  location,  the  standard  deviation  of 
the  least  squares  -fit  (sigma)  ,  the  group  speed  (which 
should  be  around  1440  m/s) ,  and  the  y  intercept  of  the  time 
delay  /  slant  range  plot.   After  the  final  answer  the  user 
is  asked  whether  or  not  he  would  like  to  remove  outlying 
paints.   If  this  option  is  selected  the  program  removes 
hydrophones  with  a  deviation  from  the  least  squares  fit  of 
more  than  3  times  the  standard  deviation,  and  the  group 
speed  and  the  standard  deviation  are    recalculated  and 
di  spl ayed. 

The  user  is  then  asked  if  he  would  like  to  locate  the 
event  with  different  hydrophones,  and  if  so  returns  the 
user  to  the  start  of  the  channel  selection  process.   The 
user  may  repeat  this  location  scheme  as  many  times  as 
necessary  for  a  particular  event.   Once  the  user  is 
satisfied  with  the  location  answer,  and  declines  to  locate 
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the  event  with  di-f-ferent  phones,  the  program  calculates  the 
event  strength  based  on  the  amplitudes  in  the  input  file, 
the  location  of  the  event  and  spherical  spreading  losses. 
The  location  parameters  and  the  strength  are    then  written 
to  an  oat-file- 

At  this  point  the  user  is  given  the  option  to  exit  the 
program  or  locate  the  next  event  in  the  input  file.   The 
location  process  continues  until  the  user  exits  or  until 
the  last  event  in  a  file  is  located. 
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/*  location. c 

program  to  locate  the  spatial  position  of  an  event  based  on 
time  delays  taken  from  "detection" . 

Source  strength  is  also  computed. 

*/ 

^include  <stdio.h> 
#include  <math.h> 

#define  PHONES  31 

#define  LEVEL  4 

#define  FINENESS  20 

♦define  DEPTH  91.0 

♦define  SENSITIVITY  0.0000000112202 

main  () 


float  amp[PHONES],  phonex [PHONES] ,  phoney [PHONES] ; 

float  timedelay (PHONES] ,  r[PHONES],  bestrange [PHONES] ; 

float  sumtime,  sumr,  sumrsq,  sumtimer,  slope,  yintr; 

float  sigma,  bestsigma,  bestslope,  bestyintr; 

float  bestamp,  N,  source [PHONES] ,  sumsource; 

float  xgs,  ygs,  xcntr,  ycntr,  a,  b,  bestx,  besty; 

float  level,  gridsize,  xfineness,  yfineness,  time,  gpspeed; 

int  i,  j,  tape,  date,  hour,  min,  sec,  event,  flag; 

int  phone,   num,  1,  n,  m,   bestflag; 

int  phoneflag [PHONES ] ; 

int  eventselect,  rerun,  change,  answer,  bye,  quadrant; 

char  iname[80],  oname[80]; 

FILE  *ptr,  openO,  *locptr,  *optr; 

/*  PROGRAM  INITIALIZATION  */ 

bye  =•  0; 

fprintf (stderr, "input  file  =  \n")  ; 
3canf("%s",  iname) ; 

fprintf (stderr, "output  file  =  \n")  ; 
scanf("%s",  oname) ; 

/*  open  files  */ 

if((ptr  -  f open (iname,  "r") )  =-  NULL)  { 

fprintf (stderr, "can' t  open  %3\n",  iname); 

exit  (1) ; 
) 
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if((optr  -  fopen(oname,  "w") )  ==  NULL)   { 
fprintf (stderr, "can' t  open  %s\n",  oname) ; 
exit (1)  ; 

) 

if((locptr  -   fopen ("array_loc",    "r") )    ==  NXHiL)       { 
fprintf (stderr, "can' t   open   array_loc    file\n") ; 
exit (1) ; 

} 


/*  read  hydrophone  locations  into  array  */ 

for  (i=l;  i<PHONES;  i++)   { 

f scanf (locptr, "%d  %f  %f",&phone,  &phonex[i],  &phoney[i]); 
) 


/*  read  input  file  header  */ 

fscanf(ptr,  "%d  %d  %d  %d  %d",  Stape,  Sdate,  Shour,  Smin,  Ssec) 
/*  read  event  header  */ 

f scanf (ptr,  "%d  ",  Sf lag) ; 

while (bye  !--!){ 


eventselect  =  0; 

while (eventselect  !"  1)   ( 

if(flag  <  0) 
exit (0) ; 

for  (i=l;  i<PHONES;  i++)   { 

phoneflag[i]  »  0; 
) 

f scanf (ptr,  "%d  %f",  Sevent,  Stime) ; 

for(i=l;  i<PHONES;  i++)  { 
fscanf (ptr, "%d", &f lag) ; 

if (flag  >  0)  { 

fscanf (ptr,  "%d",  Sj); 

phoneflagtj]  »  1; 

fscanf (ptr,  "%f  %f",  Stimedelay [ j] ,  4amp[j]); 
) 
else  { 

i  -  PHONES; 
) 
) 
fprintf (stderr, "event  =  %d,  time  =  %f\n",  event,  time) ; 
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fprintf (stderr, "Do  you  wish  to  locate  this  event?  (1  =  yes)\n"); 
scanf("%ci",  fieventselect)  ; 


1 


/*  channel  selection   */ 

rerun  •  1 ; 

while  (rerun  »-  1)  { 

change  -  1; 

while  (change  ==1)   { 

fprintf (stderr, "    phone  delay     \n") ; 

for(i=l;  i<PHONES;  i++)   { 
if (phoneflag[i]  !=  0)   { 

fprintf (stderr, "     %d  %f   \n", i, timedelay [i] ) 

) 
) 

fprintf  (stderr, "Do  you  wish  to  change  status?  (1  =  yes)\n"); 
3canf("%d",  Schange) ; 
if  (change  »=  I)   { 

fprintf (stderr, "change  status  by  typing  phone#\n") ; 
fprintf (stderr, "type  -1  to  quit\n") ; 
scanf ("%d", s j) ; 

while  (j  !=  -1)  ( 

if  (phoneflagC j]  !=  0)   { 

phoneflagCj]  =  0; 
} 
else   { 

phoneflagCj]  =  1; 
1 

scanf ("%d",  & j) ; 
} 
) 
) 


/*  locate  event  */ 

num  -  0 ; 
sumtime  -  0.0; 

for  (i-1;  i<PHONES;  i++)   { 

if  (phoneflag[i]  !=-  0)   ( 
num++ ; 
sumtime  +=  timedelay [i] ; 

) 
) 

fprintf (stderr, "select  quadrant  to  search  (1=NE,  2=NW,  3=SW,  4=SE\n") 
fprintf (stderr, "5=N,  6=S,  7=E,  8=W) \n") ; 
scanf ("%d", Squadrant) ; 
if  (quadrant  -=  1)  { 
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xcntr  =  10000.0; 

ycntr  -  10000 .0; 
} 
else  if  (quadrant  ==  2)  { 

xcntr  -  -10000 .0; 

ycntr  =-  10000.0; 
} 
else  if  (quadrant  "3)   { 

xcntr  =  -10000.0; 

ycntr  -  -10000.0; 
) 
else  if  (quadrant  —4)   { 

xcntr  =•  10000.0; 

ycntr  =  -10000 .0; 
} 
else  if  (quadrant  ==  5)    ( 

xcntr  =  0.0; 

ycntr  =  10000.0; 
} 
else  if  (quadrant  =-  6)    ( 

xcntr  =•  0.0; 

ycntr  -  -10000.0; 
) 
else  if  (quadrant  ==  "7)    { 

xcntr  -  10000.0; 

ycntr  =■  0.0; 
) 
else  if  (quadrant  ==8)    { 

xcntr  -  -10000.0; 

ycntr  =•  0.0; 
} 

else  { 

xcntr  =  0.0; 

ycntr  =  0.0; 
} 

bestflag  =  0; 

for  (1=0;  KLEVEL;  1++)  { 
level  =  LEVEL-1-1; 
gridsize  =  pow (10 . 0, level) ; 

for(m=0;  m<FINENESS;  m++)  ( 
yfineness  =  m  -  (FINENESS/2); 
ygs  -  ycntr  +  yfineness  *  gridsize; 

for(n-0;  n<FINENESS;  n++)  { 
xfineness  =■  n  -  (FINENESS/2); 
xgs  -  xcntr  +  xfineness  *  gridsize; 
sumr  =  0.0; 
sumrsq  "  0.0; 
sumtimer  «  0.0; 
for  (i-1;  i<PHONES;  i++)   { 
if (phoneflagCi]  !-  0)   { 
a  -  xg3-phonex[i] ; 
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b  =  ygs-phoney [i] ; 

r[i]  =  sqrt (pow(a,2.0)  +  pow(b,2.0)  +  pow (DEPTH, 2 . 0) ) 

sumr  +—  r [i] ; 
sumrsq  +-  pow (r [i] , 2 . 0) ; 
sumtimer  +»  (r[i]  *  timedelay [i] ) ; 
} 


N  =»  num; 

slope  »  (  (N  *  sumtimer)  -  (sumr  *  stjmtime)  )  / 

(  (N  *  sumrsq)  -  pow  (sxomr,  2  .  0)  )  ; 
yintr  =  (sxomtime  -  (slope  *  sumr)  ) /N; 


sigma  =  0.0; 

for(i=l;  KPHONES;  i++)   { 

if (phoneflag[i]  !=  0)   { 

sigma  +=  pow ( (timedelay [i] -yintr- (slope  *  r[i])),2.0) 

) 
) 
sigma  «  sqrt (sigma/N) ; 

if  (bestflag  =-  0)   ( 
bestx  -  xgs; 
besty  -  ygs; 
bestsigma  -  sigma; 
bestslope  »  slope; 
bestyintr  =  yintr; 
for  (i=l;  i<PHONES;  i++)   { 
if  (phoneflag  !=  0)   I 
bestrange[i]  =  r[i]; 
) 
} 

bestflag  =  1; 
1 

else   ( 

if  (sigma  <  bestsigma)   ( 
bestx  =  xgs; 
besty  =  ygs; 
bestsigma  =  sigma; 
bestslope  =  slope; 
bestyintr  -  yintr; 
for  (i=l;  KPHONES;  i++)   { 
if  (phoneflag  !-  0)   { 
bestrangeCi]  =  r[i]; 
} 
} 
) 
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gpspeed  «  1.0/bestslope; 


fprintf  (stderr,  "bestx  =■  %f,  besty  =  %f,  sigma  =  %f\n", 

bestx,  besty,  bestsigma) ; 
fprintf  (stderr,  "group  velocity  =•  %f\n'',  gpspeed)  ; 
fprintf (stderr, "y  intercept  =  %f\n",  bestyintr) ; 
xcntr  =  bestx; 
ycntr  =  besty; 


fprintf (stderr, "Do  you  wish  to  remove  outlying  points?  (l=yes)\n"); 
scanf ("%d", Sanswer) ; 
if  (answer  »=  1)  { 

num  -0; 

sumr  =0 . ; 

sumtime  -  0 . ; 

sumrsq  -  0 . ; 

suititimer   »    0  . ; 

for  (i-1;  i<PHONES;  i++)  { 
if (phoneflag[i]  !-  0)  ( 

if  (sqrt (pow ( (timedelay [i] -bestyintr- (bestslope*bestrange [i] ) 
2.0))  <  2 .5*bestsigma)  { 
nuir++; 

sumtinie  ■*-=  timedelay  [i]  ; 
sumr  +=•  bestrange  [i]  ; 
sumrsq  +=  bestrange [i] *bestrange [i] ; 
sumtimer  +=  bestrange [i] *timedelay [i] ; 
} 
else  ( 

fprintf (stderr, "outlying  phone  #  %d\n",i); 
1 
) 
) 

N  =  num; 
bestslope  =■  (  (N'sumtimer)  -  (sumr*sumtime)  )  / 

(  (N*3umrsq)  -  (sumr*sumr)  )  ; 
bestyintr  -  (sumtime- (bestslope'sumr) ) /N; 
gpspeed  -  1.0/bestslope; 
fprintf (stderr, "bestx  =  %f,  besty  =  %f,  sigma  =  %f\n", 

bestx,  besty,  bestsigma) ; 
fprintf (stderr, "group  velocity  »  %f\n",  gpspeed); 
fprintf (stderr, "y  intercept  =  %f\n",  bestyintr); 
) 

fprintf (stderr, 

"Do  you  wish  to  relocate  with  different  phones'?  (l=yes)\n"); 
scanf ("%d", irerun) ; 
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/*  finding  source  amplitude   */ 

for(i=l;  i<PHONES;  i++)   { 
if  (phoneflag[i]  !-  0)   { 

30urce[i]  -  (amp [i] /SENSITIVITY)  *  bestrange [i] ; 
sumsource  +=  source [i]; 
) 
} 
bestamp  =  sumsource/N; 

fprintf (optr, "%d  %d  %f  %f  %f  %f\n",  event,  date,  time, 

bestx,  besty,  bestamp) ; 
fprintf (optr, "%f  %f  %f\n",  bestsigma,  gpspeed,  bestyintr) 

fprintf (stderr, "Do  you  wish  to  exit?  (-1  =  exit)\n"); 
scanf("%d",  Sbye) ; 
) 

f close (ptr) ; 
f close  (optr) ; 
fclose  (locptr) ; 

exit  ( 0 )  ; 

) 
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APPENDIX  C 


Table  C-1:   Event  Location  Summary 


Table  C— 2:   Tape  Summary 


Table  C— 3:   Event  Inter arrival  Time  Summary 
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209 

0.0070 

10.29 

1473 

13 

92 

0.0036 

5.22 

1459 

c 

187 

0.0032 

4.62 

1452 

none 

12 

0.0143 

19.45 

1364 

13,20,21 

71 

0.0107 

15.01 

1409 

none 

74 

0.0065 

9.23 

1416 

I  d 

24 
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>5P,i 


jCbtion  buiTiiiiary 


rkse 

s; 

y 

r 

R 

phi 

siqiTiS 

SlQiTia 

c 

HydropfsGries 

(Tape  t. 

ReiiiQved 

event  Si 

(in! 

(a) 

i«) 

int) 

:  degrees) 

(sec) 

(iri) 

!  IT! /sec) 

4011, 33t» 

2390 

-699 

2973 

2975 

104 

0.0207 

27.41 

1  ^-L 

none 

4011,166 

-922 

-6001 

6071 

6072 

189 

0.0052 

7.ol 

1462 

nane 

4011,169 

-6201 

5340 

3183 

3134 

311 

0.0060 

3.75 

1443 

none 

4011,215 

-96 

-391 

403 

413 

194 

0.0042 

6.26 

1492 

24 

4011,241 

2390 

-792 

2997 

2993 

105 

0.0113 

16.33 

1431 

none 

4011,250 

-420 

890 

984 

939 

335 

0.003E 

5.&0 

1491 

-? 

4011,259 

-531 

TSQA 

7^03 

TSi'iS 

0,0035 

3.07 

1455 

r  -TiP 

4011,273 

-584 

-1001 

1159 

1163 

210 

0.0173 

24.27 

1401 

7.11 

4011,283 

2190 

-3011 

3723 

3724 

144 

0.0032 

4.62 

1451 

7 

4011,307 

3260 

5890 

6732 

6733 

29 

0.0060 

8.84 

1463 

5,6,7 

4011,344 

1990 

766 

2132 

2134 

69 

0.0069 

10.29 

1431 

15,13 

4011,362 

JJ/0 

-10007 

11454 

11455 

151 

0.0097 

14.17 

1457 

:3 

4011,335 

-329 

-4001 

4036 

4037 

192 

0.0077 

11.22 

1463 

none 

4011,510 

19870 

13790 

24136 

24137 

55 

0.0169 

23.52 

1393 

J,6,  /  .l-'.j.l', 

4011,532 

-561 

244 

612 

619 

294 

0.0063 

9.93 

1475 

none 

4011,541 

-253 

607 

660 

666 

337 

0.0034 

11.90 

1423 

15, IB 

4011,573* 

14471 

-16409 

21373 

21379 

139 

0.0056 

3.31 

1573 

6,15.16.19.: 

4011.595 

35S 

-il 

_"  w'  / 

--n 

103 

0.0i02 

14.37 

140S 

'" 

4011, ai5 

6212 

-4435 

7633 

/633 

126 

0.u020 

2.32 

1435 

!3,i3 

4011.616 

-630 

2890 

2953 

2959 

34S 

0.0232 

33.55 

1447 

none 

4011,663 

2990 

1566 

3375 

3377 

0.0053 

7.31 

1479 

18 

4011,664 

265 

256 

368 

380 

46 

0.0035 

4.37 

1393 

4011,693 

30 

6790 

6790 

o'91 

I 

0,0y27 

3.33 

1452 

21 

4011.709* 

3446 

1081 

3612 

■613 

3 

O.OOdb 

V.41 

1:75 

-.■■| 

4011.721 

-206 

-249 

-  J  L  ,> 

iji 

0.0i32 

i-b  .  l'j 

1444 

none 

4011,311 

-396 

390 

974 

979 

3  OS 

j.0079 

11.75 

1492 

7 

4011,361 

i  -J  0  C  d 

-15390 

20960 

209i0 

139 

0.0072 

1.  ■-.'  t  i  i. 

1411 

9,12,13.22 

4011,863 

734 

-701 

1  '<>  1  J 

1019 

134 

o.uii: 

15.67 

1400 

ncne 

4011,335 

990 

345 

1048 

1053 

71 

0.0053 

7.41 

1401 

18,21.22 

4011,965 

6341 

-8212 

10375 

10376 

142 

0.0039 

5.61 

1423 

none 

4011.992 

-1301 

1390 

2294 

2296 

0.0233 

33.52 

1406 

7,23,24 

4011,1046 

-655 

1390 

2000 

2002 

341 

0.0049 

7.26 

1487 

7,12.13,23,: 

4011,1061 

-52 

107 

119 

151 

334 

0.0020 

2.70 

1333 

3,4,5,6,20 

4011,1128 

-636 

452 

780 

786 

305 

0.0047 

6.72 

1441 

none 

4011,1137 

-3001 

-133 

3004 

3005 

267 

0.0097 

13.79 

1416 

23,24 
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Table  C-1    Event  Location  Susisiary 


Naie 

x 

y 

r 

R 

phi 

siggia 

signa 

c 

Hydrophones 

(Tape  i. 

Removed 

event  1) 

(il 

(!) 

(II) 

(a) 

(degrees) 

(sec) 

(a) 

(is/sec) 

4013,120 

-833 

-4970 

5039 

5040 

190 

0.0091 

13.11 

1447 

none 

4013,177 

370 

-1896 

1932 

1934 

169 

0.0035 

5.00 

1437 

2,8,9,10 

4013,277 

-393 

366 

537 

545 

313 

0.0045 

6.38 

1418  8,9,14,16,22,23, 

4013,317 

1326 

-12 

1326 

1329 

91 

0.0028 

3.99 

1445 

none 

4013,351 

-54 

-99 

m 

146 

209 

0.0070 

9.72 

1396 

none 

4013,354 

1534 

761 

1712 

1715 

64 

0.0016 

2.31 

1447 

21 

4013,373 

216 

683 

71o 

722 

13 

0.0028 

4.11 

1450 

none 

4013,381 

1248 

-2815 

3079 

3081 

156 

0.0051 

7.41 

1442 

none 

4013,403 

2602 

-5450 

6039 

6040 

154 

0.0067 

9.56 

1420 

none 

4013,464 

342 

2261 

2413 

2414 

20 

0.0083 

12.03 

1454 

none 

4013,478 

-52 

662 

664 

671 

356 

0.0076 

10.56 

1388 

none 

4013,518 

175 

492 

522 

530 

20 

0.0082 

12.05 

1472 

14,16,17 

4013,661 

6390 

3658 

7801 

7801 

62 

0.0030 

4.34 

1463 

none 

4013,664 

-1001 

1205 

1567 

1569 

320 

0.0114 

16.72 

1473 

19,20,22 

4013,671 

19866 

-8441 

21585 

21585 

113 

0.0104 

14.72 

1416 

none 

4013,694 

-675 

9290 

9314 

9315 

356 

0.0039 

5.62 

1437 

none 

4013,723 

3100 

-1533 

3458 

3460 

116 

0.0059 

8.51 

1447 

none 

4013.755* 

264 

168 

313 

326 

58 

0.0151 

23.00 

1520 

24 

4013,776 

7989 

-5590 

9750 

9751 

125 

0.0052 

7.48 

1446 

none 

4013,797 

-339 

7013 

7024 

7024 

357 

0.0052 

7.56 

1461 

none 

4013,801 

-400 

-133 

422 

432 

252 

0.0078 

11.58 

1479 

none 

4013,866 

38 

307 

309 

323 

7 

0.0062 

8.92 

1448 

none 

4013.875 

49 

320 

324 

337 

9 

0.0031 

4.27 

1397 

none 

4013.908 

9854 

17869 

20406 

20406 

29 

0.0066 

9.26 

1414 

7.12 

4013,922* 

13 

26 

29 

97 

27 

0.0090 

13.61 

1514  8,9,10,U,!2,18 

4013.950 

1165 

-2253 

2536 

2538 

1  C-T 

0.0037 

5.30 

1442 

none 

4013,979 

7890 

920 

7943 

7944 

83 

0.0084 

11.99 

1432 

none 

4013,1014 

-143 

2890 

2894 

2895 

357 

0.0091 

12.80 

1413 

none 

4013,1018 

-137 

574 

590 

597 

347 

0.0169 

24.97 

1476 

none 

4013,1107 

-701 

453 

835 

840 

303 

0.0081 

11.51 

1427 

none 

4013,1124 

-9001 

8142 

12137 

12137 

312 

0.0064 

9.43 

1478 

none 

4013,1166 

870 

-7001 

7055 

7055 

173 

0.0050 

7.39 

1468 

18 

4013,1187 

-696 

693 

982 

987 

315 

0.0063 

9.34 

1472 

17,18 

4013,1203 

7789 

-2315 

8126 

8126 

107 

0.0035 

5.11 

1480 

none 

4015,232 

-lOi 

422 

434 

444 

347 

0.0026 

3.62 

1406  21,22 

4016    no 

events 
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Table  C-l  .  Event  Location  SuinBiary 


Naae 

X 

y 

r 

R 

phi 

sigsa 

siqia 

c 

Hydrophones 

(Tape  t. 

Removed 

event  1) 

(») 

(i) 

(a) 

(9) 

(degrees) 

(sec) 

(in) 

(ffl/sec) 

4019,178 

-8416 

-18097 

19958 

19958 

205 

0.0038 

5.52 

1473 

none 

4019,679* 

-17789 

10549 

20682 

20682 

301 

0.0023 

3.50 

1533 

none 

4019,689 

-20109 

11680 

23255 

23255 

300 

0.0026 

3.70 

1428 

none 

4019,776 

-21087 

12642 

24586 

24586 

301 

0.0037 

5.13 

1394 

none 

4019,841* 

-3001 

-88 

3002 

3004 

268 

0.0061 

9.10 

1503 

none 

4021,61 

-7356 

15998 

17608 

17608 

0.0103 

14.68 

1430 

none 

4021,154 

-6508 

18999 

20083 

20033 

341 

0.0062 

8.84 

1436 

13 

4023    no 

events 

4024    no 

events 

4027,28         -1001  -371 


1068 


1072 


250      0.0145       20.81 


1434  24 


4029,666 
4029,897 
4029,1059 


-4250 
347 

-3313 


-10010 

-110 

-9110 


10875 

364 

9695 


10875 
376 

9696 


203  0.0075 
108  0.0244 
200   0.0040 


11.20 

33.82 

5.71 


1489  3,6,11,13,29 
1385  none 
1439  29 


4031,823 


159 


-373 


405 


416 


157   0.0076 


11, 


1493  none 


4033,490 

4033,1041 

4033,1043 

4033,1053 

4033,1126 

4033,1173 

4040,26 

4040,148 

4040,385 

4040,398 

4040,714 

4040,730 

4040,871 

4040,1114 

4040,1157 


4047,984 


-2532 
-21010 
-10309 
-285 
-21090 
-20905 

-20488 

-6001 

1462 

-9001 

-5432 

6468 

-3424 

-20973 

-7408 

-457 


4190 
767 

1017 
120 

1075 
B05 

19987 

4252 

-4010 

8110 

-14089 

■17109 

-8910 

16769 

7391 

200 


4896 
21024 
10359 
309 
21117 
20920 

28622 

7355 

4268 

12116 

15100 

18291 

9545 

26853 

10823 

499 


4897 
21024 
10359 
323 
21118 
20921 

28622 

7355 

4269 

12116 

15100 

18291 

9546 

26853 

10824 

507 


329 
272 
276 
293 
273 
272 

314 
305 
160 
312 
201 
159 
201 
309 
317 


0.0047 
0.0119 
0.0054 
0.0042 
0.0020 
0.0085 

0.0125 
0.0030 
0.0054 
0.0150 
0.0056 
0.0126 
0.0022 
0.0061 
0.0075 


6.61 
17.80 
91 


10 
91 
89 


17.63 
4.32 
8.03 

21.35 
8.38 

17.76 
3.20 
9.03 

11.26 


1408  13,18,29,30 

1493  9,11,13,19,20,21,22,23.24 

1451  16,22,23.24 

1451  13 

1465  13 


1410  6 

1426  10 
1486  none 

1427  13 
1497  none 

1411  6 
1446  none 
1469  11 
1493  none 


294   0.0230   31.93 


1387  none 


4049,1064   -2320   -8110    8435    3436     196   0.0068 
4049,1066   -1889   -8001    8221    3221     193   0.0150 


9.75 
21.01 


1426  none 
1401  9,11,13 
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Table  C-1  Event  Location  Sufflfflary 


Naie 

X 

y 

r 

R 

phi 

sigaa 

Sigma 

c 

Hydrophoi 

(Tape  t. 

Removed 

event  #) 

(91 

(IB) 

(a) 

(a) 

(degrees) 

(sec) 

(ffl) 

(m/sec) 

4051 

no  events 

4053,3 

-16878 

18965 

25388 

25388 

318 

0.0021 

3.11 

1461 

7 

4055 

no  events 

4057 

no  events 

4059 

no  events 

4061,571 

-8393 

-20923 

22544 

22544 

202 

0.0127 

18.29 

1445 

21,22,23 

4061,879 

-3355 

-2455 

4157 

4158 

234 

0.0027 

3.91 

1445 

5,6,23 

4063,296 

-10503 

17985 

20827 

20827 

330 

0.0090 

12.91 

1441 

none 

4063,824 

-16793 

14787 

22375 

22376 

311 

0.0043 

6.29 

1455 

none 

4063,864 

-2001 

1346 

2412 

2413 

304 

0.0043 

6.18 

1424 

18 

4063,933 

-9077 

7976 

12083 

12084 

311 

0.0111 

15.60 

1408 

none 

4063,936 

-1130 

912 

1452 

1455 

309 

0.0037 

5.42 

1450 

18 

4065 

no  events 

4067 

no  events 

2001,543   -3455   -5112 


6170 


6171 


214   0.0026 


3.77 


1452 


2009,189 

-1997 

-20987 

21082 

21082 

185 

0.0036 

5.40 

1487  none 

2009. i53 

S74 

-5201 

5274 

5275 

170 

0.0070 

9.  BO 

1402  18 

2009,455 

1441 

-2001 

2466 

2468 

144 

0.0045 

6.55 

1445  12 

2009,925 

-105 

-073 

681 

687 

189 

0.0193 

27.43 

1422  none 

2023,74 

15430 

-16460 

22561 

22562 

137 

0.0104 

15.13 

1458  7,11 

3001,11 

-3921 

75 

3922 

3923 

271 

0.0125 

17.92 

1434  none 

3001,16 

-20103 

-9471 

22222 

22222 

245 

0.0137 

19.92 

1450  none 

3001,301 

-3001 

698 

3081 

3083 

283 

0.0149 

21.13 

1418  4,5,6,7,20,21 

,  ii. ,  i.  J 

3047,169 

-19109 

-2181 

19233 

19233 

263 

0.0073 

10.56 

1453  23 

3047,751 

890 

-538 

1040 

1044 

N  = 

mean  = 
std  dev  = 

121 

0.0134 

199 

0.0077 
0.0055 

19.52 

199 

11.09 
7.82 

1460  none 

199 

1444 

34 

Table  C-2  Taoe  Sussiarv 


Tape  # 

HinuteE 

Julian  De 

lections 

Artiiacts 

Events 

iiult 

P'i 

e 

r  .  1  _  _ 

F 

^M- 

Es a mined 

Date 

Responst 

'5 

Alans 

Hz  Pa 

4001 

20 

36 

0 

13 

5 

ij 

4ot 

Avail 

2001 

20 

B3 

6 

0 

1 

0 

c 

Net 

R  'v  a  1 1 

2009 

20 

B9 

15 

0 

4 

1 

10 

0.022 

3001 

10 

90 

0 

3 

0 

2 

0.013 

4003 

20 

91 

15 

0 

13 

1 

i 

0  =  044 

4005 

20 

91 

24 

5 

11 

2 

6 

0.035 

4007 

;7 

91 

11 

I'i 

6 

1 

4 

0.022 

iOOv 

i-'i 

"i: 

55 

"7 

^7 

8 

- 

0  .  'i  i  0 

4011 

20 

92 

47 

7 

35 

1 

4 

0.010 

4013 

20 

93 

43 

1 
i 

34 

5 

3 

0.014 

2023 

20 

98 

1 

0 

1 

0 

0 

0.037 

4015 

"T' 

99 

14 

8 

1 

0 

5 

0,040 

3047 

17.5 

103 

0 

.■•i 

2 

0 

0 

0.010 

40  U 

20 

105 

3 

■? 

0 

0 

0 

0.C17 

4019 

20 

105 

14 

T 

5 

1 

5 

0.016 

4021 

20 

109 

6 

0 

2 

2 

n 

0,013 

4023 

20 

109 

22 

20 

0 

0 

2 

0.011 

4024 

20 

109 

14 

11 

0 

0 

3 

o.on 

4027 

20 

no 

La 

1 

0 

0 

0.012 

4029 

20 

no 

11 

8 

3 

0 

0 

0.012 

4031 

20 

110 

11 

10 

1 

0 

0 

0.012 

4033 

20 

no 

19 

12 

6 

0 

1 

0.0i2 

4040 

20 

111 

10 

9 

3 

0.034 

4047 

L'-J 

ill 

6 

0 

1 

0 

0.114 

40^9 

20 

{  ■  i 

J 

n 

3 

^J 

0.140 

4051 

20 

ill 

6 

0 

0 

0 

Q 

0.140 

4053 

9n 

112 

0 

1 

- 

le 

0.080 

i055 

i  :   1  -J 

112 

0 

0 

0 

0 

0 

0,032 

4057 

20 

112 

9 

1 

0 

0 

8 

0.053 

4059 

20 

112 

1 

0 

0 

0 

1 

0,065 

4061 

20 

112 

6 

0 

2 

0 

4 

0.034 

4063 

20 

112 

9 

0 

5 

7 
•J 

1 

0.023 

4065 

20 

112 

0 

2 

0 

0 

1 

0.027 

4067 

20 

112 

6 

6 

0 

0 

0 

0.013 

Total 


662 


499 


139 


199 


36 


12= 


Detscticn  is  asclarea  vsnen  -t  least  Wi.   3t  the  hyarophGnas  nave  tiigaeci. 

Artiiacts  are  iiian-Qade  noises  Uhots,  air  gun  blasts). 

Events  ar?  decEcticns  atrcnc  enough  tQ  analyze. 

Multiple  Responses  are  Fiultiple  detecticns  oi  an  alreadv  counted  event. 

False  alaras  are  detections  too  Neak  to  analyze. 


laoie  u-i  hv5 

'Ml   Interirnval  i 

bue  'dusifflary 

Nase 

Event  Tiae 

Interarrivsl 

Tiie   Hori: 

antal  Range 

(Taoe  1,  event  #) 

(sec) 

(sec)  (205 

bin) 

(?,) 

4001,63 

63 

63 

3 

25310 

4001,126 

126 

63 

3 

6901 

4001,175 

175 

49 

2 

9975 

4001,177 

177 

2 

0 

10241 

4001,195 

195 

18 

0 

2271 

4001,653 

653 

463 

23 

2109 

4001.818 

31E 

160 

3 

11615 

4001,393 

895 

77 

•j 

10642 

4001,1066 

1066 

171 

3 

20136 

4001,1131 

113! 

65 

3 

13391 

4001,1156 

1156 

25 

1 

3196 

4001, 1P4 

1194 

38 

1 

23204 

4001,1215 

1215 

21 

1 

19011 

4003,21 

21 

21 

1 

358 

4003,52 

en 
■Ji. 

31 

1 

9716 

4003,283 

233 

231 

11 

230 

4003,234 

284 

1 

0 

1276 

4003,285 

235 

1 

« 

0 

341 

4003,301 

301 

16 

() 

432 

4003,419 

419 

113 

5 

4954 

4003,516 

516 

97 

4 

3320 

4003,634 

634 

US 

5 

7631 

4003,635 

635 

1 

0 

3715 

40j3.703 

703 

73 

'\ 

3650 

4003,995 

995 

287 

14 

20076 

4003,1053 

1053 

58 

' 

272 

4005,219 

219 

21" 

10 

3902 

4005,24! 

241 

22 

1 

206 

4005,395 

395 

154 

7 

953 

4005.303 

308 

413 

20 

309 

4005,916 

916 

103 

J 

998 

4005,997 

997 

31 

4 

13006 

4005,1012 

1012 

15 

0 

2159 

4005,1097 

1097 

85 

4 

20083 

4005.1169 

1169 

72 

3 

1979 

4005,1173 

U73 

4 

0 

143 

4005,1183 

1183 

10 

0 

111 

4007,204 

204 

204 

iC 

364 

4007,444 

444 

240 

12 

798 

4007,582 

532 

133 

6 

218 

4007,794 

794 

212 

10 

210 

4007, S20 

52U 

-i 

1 

584 

4007. tS3 

936 

led 

3 

35 

■154- 


Event  Intsrarrival  Time  Susiinarv 


iNaie  Event  TiiBe 

(Tape  S,  event  S!   (sec) 


Interarnval  lite 
(sec)  (205  bin) 


Horizontal  Range 
(ffli 


4009,100 

4009,172 

4009,173 

4009,193 

4009,194 

4009,278 

4009,292 

4009=300 

4009,380 

4009,394 

4009,385 

4009,397 

4009,410 

4009,447 

4009,476 

4009,493 

4009,579 

4009,600 

4009,603 

4009,620 

4009,621 

4009.626 

4009,667 

4009.677 

4009,728 

4009,736 

4009,302 

4009,317 

4009,828 

4009,373 

4009,916 

4009,948 

4009,949 

4009,967 

4009,932 

4009,1041 

4009,1160 


100 

i  •  i 

178 
193 
194 
278 
292 
300 
380 
334 
335 
397 
410 
447 
476 
493 
579 
600 
603 
620 
621 
626 
667 
677 
72s 
736 
802 
iil/ 
823 
873 
916 
948 
949 
967 
982 
1041 
1160 


uu 

J 

72 

3 

6 

0 

15 

0 

1 

0 

34 

4 

14 

0 

3 

0 

80 

4 

4 

0 

1 

0 

12 

0 

13 

0 

37 

I 

29 

I 

17 

0 

36 

4 

21 

I 

3 

0 

12 

0 

1 

0 

5 

0 

41 
to 

0 

16 

0 

15 

0 

11 

0 

50 

38 

1 

32 

I 

1 

0 

18 

0 

15 

0 

59 

2 

119 

■J 

3955 

370 

987 

2785 

1160 

12253 

16051 

V  9  "^ 

11046 

979 

3339 

3174 

492 

7887 

883 

3651 
37 

1971 

3376 
693 
737 

2367 
176 

385 

3597 

7934 

966 

198 

229 

3419 

9004 

7565 

499 

1997 

573 


-155- 


)l5  C-3    Event  Intararrival  Tii?,e  Susiffiarv 


Hdae 

Event  Tiffie 

{Tape  S,  event  S; 

(sec) 

4011,33 

33 

4011,168 

163 

4011,169 

169 

4011,215 

215 

4011,241 

241 

4011,250 

250 

4011,259 

259 

4011,278 

273 

4011,283 

283 

4011,307 

307 

4011,344 

344 

4011,362 

362 

4011,335 

385 

4011,510 

510 

4011,532 

532 

4011,541 

541 

4011,573 

573 

4011,595 

595 

4011.615 

0  1  J 

4011.616 

616 

4011,663 

663 

4011,664 

664 

4011,693 

693 

4011 J09 

709 

4011.721 

?- ! 

4011,311 

311 

4011.861 

861 

4011,363 

363 

4011,885 

885 

4011,985 

935 

4011,992 

992 

4011.1046 

1046 

4011,1061 

1061 

4011,1128 

1123 

4011,1187 

1187 

nterarri 

val  Tiiiie 

Korizontal  Range 

i5ec) 

!205  bin) 

(ia) 

83 

4 

2973 

35 

4 

3071 

1 

0 

3133 

46 

2 

403 

26 

1 

2997 

9 

0 

934 

9 

0 

7908 

19 

■J 

.  \Xil 

5 

0 

3723 

24 

I 

6732 

37 

1 

2132 

18 

0 

11454 

23 

1 

4036 

125 

6 

24136 

22 

1 

612 

9 

0 

660 

32 

1 

21378 

22 

1 

1  <_-: 
j  "J  i 

20 

:'600 

1 

0 

2953 

47 

2 

3375 

1 

0 

363 

29 

t 

6790 

16 

i  -. 

0 

3612 

90 

974 

50 

20960 

2 

0 

1015 

22 

1 

1048 

100 

5 

1U375 

7 

/ 

0 

2294 

54 

2000 

15 

0 

119 

67 

3 

780 

59 

L. 

3004 

-15i- 


ibis  C-3  Event  Interarrival 


iufTiSisrv 


Naae  Event  TiiTie 

('Tape  #.  tvsnt  #*   (sec) 


4013,120 
4ui3,177 
4013, 
4013, 


77 


4013.354 

4013,373 

4C13,3dl 

4013,403 

4013,464 

4013,478 

4013,518 

4013,661 

4013,664 

4013,671 

4013,694 

4013,723 

4013,755 

4013,776 

4013,797 

4013.801 

4013, 36i) 

4013.375 

4013.903 

4013.922 

4013,950 

4013,979 

40L3.I014 

4013,1018 

4013,1107 

4013,1124 

4013,1166 

4013,1187 

4013.1203 


120 

177 

277 

317 

351 

354 

373 

331 

403 

464 

478 

513 

661 

664 

671 

694 

723 

755 

776 

797 

801 

366 

875 

908 

922 

950 

979 

1014 

i0l8 

1107 

1124 

1166 

1187 

1203 


Interarni 

val  Tiffie 

HGriicntal  Range 

(sec! 

(205  Din) 

■,iR) 

120 

6 

5039 

57 

1 

1932 

100 

5 

537 

40 

2 

1326 

34 

1 

113 

v' 

0 

1712 

19 

0 

716 

6 

0 

3079 

22 

1 

6039 

hi 

•J 

2413 

14 

0 

664 

4v 

2 

522 

143 

7 

7801 

3 

0 

1567 

7 

0 

21585 

23 

1 

9314 

29 

1 

3458 

1 

313 

21 

1 

9750 

1 

7024 

4 

0 

422 

65 

■J 

309 

n 

7 

0 

324 

•JO 

i 

20406 

14 

{} 

29 

23 

i 

iJOO 

i7 

1 

7943 

Tc; 

I 

2394 

4 

0 

590 

39 

4 

335 

17 

0 

12137 

42 

1 

7055 

21 

1 

982 

16 

0 

3126 

^u 


1=: 


?7,? 


11 


434 


4016 


none 


Tibia  C-3  Event  Interamval  Tiine  Sumffliry 


Nate 


tvent  iise 


ilaoe  i,  evant  I)   (sec) 


4019,178 
4019,679 
401?, 669 
4019,776 
4019,341 

4021,81 


178 
679 
689 
776 
341 

31 


Interamval 

Tiae 

Hg 

ri; 

^ontal  Rang 

(5SC) 

(20 

5  bin) 

(iT,) 

17B 

8 

19958 

501 

25 

20682 

10 

0 

23255 

97 

4 

24536 

65 

■J 

3002 

31 

4 

17608 

-r- 

.1 

70093 

4023 
4024 


none 
none 


4027,23 


1068 


4029,666 
4029,897 
4029,1059 


666 

897 

1059 


666 

33 

231 

11 

162 

3 

10375 

364 

9695 


4031.; 


823 


4033,490 
4033,1041 
4033,1043 
4033,1053 

4033,1173 


490 
1041 

1043 
1053 
■  ;2^ 
1173 


490 


24 


4396 
21024 
10359 


4040,26 

4040,148 

4040,335 

4040,396 

4040,714 

4040,730 

4040,371  , 

4040,1114' 

4040.1157 

4047,984 

4049,1064 

4049,iufe6 


148 
335 
398 

714 

730 

371 

1114 

1157 

934 


237 

11 

13 

0 

316 

15 

16 

0 

141 

7 

243 

12 

43 

2 

934 

49 

1064 

trr 

28622 

7355 

4263 

12116 

15100 

18291 

9545 

26853 

10823 

499 

8435 

01- i 
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Table  C-3  Event  Interarrival  Tise  SaiTifflary 

Naae  Event  Tise      Interarrival  TiiTie   Horizontal  Range 

ilape  t,  event  f)   (sec;         (seci  (205  bin)         (lai 


4051 

none 

4053,3 

3 

4055 

none 

4057 

none 

4059 

none 

4061.571 

57! 

4061,379 

379 

4063,296 

296 

4063,324 

324 

4063,364 

864 

4063,933 

933 

4063,936 

936 

40-!!  5 

none 

4067 

a  one 

2001,543 

543 

2009,139 

139 

2009,453 

453 

2009,455 

455 

2009.925 

925 

571 

23 

308 

15 

296 

14 

523 

26 

40 

i. 

69 

3 

3 

0 

11 

0 

5 

0 

235 

14 

169 

3 

582 

29 

4157 
20927 

i-i-J.'  u 

2412 

12033 

1452 


543  27  6170 

13^  ^  21092 

26i  13  5274 

2  0  2466 

470  23  681 

1023.74            74           74  3  225c 1 


3001,11  11 

3001,16  16 

3001.301  301  235     14  3081 


3047,169  169  169      3         19233 

3047.751  751  582     29  1040 


199 
100 


itd  dev  = 


•159- 


iiDiS  L-t 


vent  itranqiih  HuiTifliarv 


(Taos  ♦.  event  I) 


iiii) 


Po 
IPa) 


siqiTid 
(Paj 


Fc 
(N) 


SlQIfid 

(N) 


Naae  is  the  FRAH  IV  tape  nusber  toHoHed  bv  the  time-into  the  tape  the  event  occurred. 

r  ia  the  honicntai  range. 

Fo  is  the  average  peak  hydrophone  pressure  tor  an  event. 

sigisa  (Pa)  is  the  standard  deviation  of  the  peak  hydrophone  pressure  within  an  event. 

Fo  is  the  average  peak  dipole  strength  calculated  for  an  event. 

sigsa  (N)  is  the  standard  deviation  of  the  dipole  strength  within  an  event. 

t   is  the  trequencv,  eKalicitlv,  the  axis  crossmc  rate. 


4001 ,126 

4001,175 

4001,177 

4001.195 

i001.i53 

4001,318 

4001,895 

4001,1131 

4001,1156 

4001.1215 


4003,52 

4003,234 

4003,285 

4003,301 

^003.419 

4003,516 

4003,634 

4003,635 

4003,708 

4005,219 
4005.395 
4005,303 
4005,916 

4005,997 

4005,1012 

4005,1169 


6901 

0,487 

0.092 

707540 

134109 

JO 

9975 

0,600 

0.109 

966109 

175081 

46 

10241 

0.545 

0.194 

1002732 

357437 

42 

2271 

0.558 

0.117 

39670^ 

195120 

42 

2109 

0.665 

0.097 

1340377 

208247 

32 

11615 

1.117 

0.213 

1413009 

270370 

61 

10642 

0.524 

0.085 

1094608 

177193 

37 

13391 

0.482 

0.085 

779610 

136750 

48 

3196 

0.495 

0.069 

1222273 

166340 

30 

19011 

O.o27 

0.101 

1229059 

198513 

40 

353 

0.233 

0.103 

59368 

13452 

7  i 

9716 

0.324 

0.063 

430173 

33712 

58 

1276 

0.291 

0.131 

137664 

30580 

64 

341 

0.362 

0.136 

619^0 

14917 

30 

432 

0.274 

0.057 

62054 

15905 

74 

4954 

0.436 

0.450 

622783 

640604 

53 

3320 

0.246 

0.046 

342405 

63394 

53 

76B1 

0.163 

0.048 

235652 

69027 

CT 

3715 

0.242 

0.044 

369396 

a6290 

49 

3650 

0.225 

0.063 

264693 

73731 

64 

3902 

0.305 

0.271 

353094 

314253 

64 

953 

0.312 

0.170 

153324 

70515 

64 

309 

0.465 

0.193 

79016 

27448 

74 

998 

0.278 

0.084 

129213 

27205 

71 

18006 

0.243 

0.059 

910298 

215630 

21 

2159 

0.230 

0.060 

56416S 

147537 

27 

1979 

0.521 

0.253 

555870 

230311 

58 

4007,204 
4007,444 
40v7,S20 


0.106 
0.180 


64308 
122593 
103017 


'ibie  C-4  Event  Strenqth  Su^iTsarv 


■160- 


Masp 


:!ape  ?,  event 


4009,100 


4009,193 

4009,194 

4009,273 

4009,292 

4009, 3t0 

4009,334 

4009,335 

4009,397 

4009,410 

4009,447 

4009,476 

4009,493 

4009,579 

4009,603 

4009,620 

4009.621 

4009,667 

4009,723 

4009,736 

4009,302 

4009,317 

4009,323 

4009 , 948 

4009,949 

4009,967 

4009,1041 

4009,1160 


r 

Po 

sioiTia 

Fo 

sigiTia 

t 

iffl) 

(Pa; 

(Pa) 

(Ni 

(N) 

(Hz) 

3955 

0.217 

0.043 

237677 

47362 

69 

370 

0.449 

0.240 

34708 

34036 

71 

987 

0.256 

0.103 

150340 

52260 

53 

2785 

0.234 

0.055 

261773 

57023 

64 

1160 

0.321 

0.060 

191163 

28220 

65 

12253 

0.267 

0.045 

339643 

56951 

61 

16051 

0.243 

0.036 

315602 

46536 

60 

11046 

0.260 

'J  •  V  ■  J  1 

33S321 

4970V 

oO 

979 

0.311 

O.IOS 

I59v35 

50432 

63 

3339 

0.256 

0.049 

322932 

62179 

59 

3174 

0.223 

0.083 

250131 

9B245 

66 

492 

0.478 

0.179 

131094 

33345 

65 

7887 

0.293 

0.043 

461103 

75739 

49 

883 

0,400 

0.290 

139730 

61426 

73 

2333 

0.173 

0.028 

199141 

31577 

59 

3651 

0.343 

0.075 

373537 

31453 

69 

1971 

0.218 

0.074 

230200 

34135 

5^^ 

•3376 

0.237 

0.047 

315200 

61634 

56 

693 

0.236 

0,065 

53243 

18204 

'^P 

2367 

0.323 

0.195 

414294 

234724 

54 

2032 

0.375 

0.132 

333341 

138491 

36 

335 

0.981 

0.466 

224139 

82413 

64 

3597 

0.746 

0.6iS 

1050835 

355095 

0  j 

7934 

0.313 

0.039 

459710 

123412 

'iZ 

966 

0.179 

0.035 

105661 

22672 

58 

3419 

0.277 

0.039 

323295 

46413 

64 

9004 

0.244 

0.043 

294434 

57917 

6"^ 

7565 

0.339 

0.060 

403141 

72556 

a4 

1997 

0.224 

0.009 

221127 

70482 

64 

573 

0.425 

0.132 

137233 

31215 

63 

'icie  C-4    Event  Strength  Suminarv 


■161- 


Nas5 

r 

Po 

siqifia 

Fo 

siqiiia 

i 

(Tape  ?.  event  S) 

iff.) 

(Pa) 

(Pa) 

(N) 

(N) 

(Hz) 

4Cli.l£)8 

6071 

0.193 

0.040 

234221 

46793 

65 

4011,169 

31S3 

0.222 

0.052 

276205 

64443 

62 

4011,215 

403 

0.276 

0.079 

63100 

10454 

67 

4011,241 

2997 

0.302 

0.123 

323546 

133350 

68 

4011,250 

984 

0.353 

0.162 

193201 

65473 

59 

4011.259 

7908 

0.395 

0.055 

491277 

63119 

62 

4011,278 

1159 

0.273 

0.102 

163349 

53900 

66 

-011,233 

3723 

V.  ijj 

0.042 

239550 

42709 

74 

4011,30? 

7800 

0.244 

0.100 

307771 

126135 

61 

4011,344 

2132 

0.139 

0.041 

213914 

47390 

56 

4011,362 

11454 

0.170 

0.033 

202390 

33375 

65 

40! 1,335 

4086 

0.162 

0  Mb 

207993 

45S67 

59 

4011,532 

612 

0.316 

0.120 

95187 

24435 

66 

4011.541 

660 

0.133 

0.033 

53321 

14260 

70 

4011.595 

367 

0.708 

0.200 

157508 

32376 

64 

4011,615 

7633 

0.253 

0.056 

291097 

64413 

67 

4011,616 

2958 

0.314 

0.416 

486903 

607412 

4fc 

4011,663 

3375 

0.217 

0.037 

227709 

33947 

'1 

4011.664 

366 

0.377 

0.033 

3179^ 

14121 

A  A 

4011,693 

6790 

0.266 

0.051 

365550 

70710 

56 

4011,721 

323 

0.245 

0.067 

52222 

17894 

"J" 

4011,811 

974 

0.244 

0.056 

120519 

22495 

65 

4011.863 

1015 

0.166 

0.033 

167514 

46703 

7i 
•J  J 

4011.3S5 

1043 

0 .  362 

0.067 

244063 

4646.i 

54 

4001.985 

10375 

0.23S 

0.042 

...  i  J  0  I  u. 

56061 

53 

4011,992 

2294 

0.292 

0.116 

320173 

117435 

61 

4011,1046 

2000 

0.215 

0.055 

229653 

53113 

53 

40ii,li2S 

730 

■J  *  ^i  J  0' 

0.082 

i:S230 

20570 

70 

4011,1137 

3004 

0.267 

0.090 

299337 

985C4 

65 

•162- 


=die  i-^    tvsnt  Strenqtn  Suuiiiary 


Naas 

r 

Po 

iigiBa 

Fo 

siqiBS 

t 

(Tape  i,  event  If 

(n) 

(Pa; 

(Pa) 

(N/ 

(N) 

(H:) 

4013,120 

5039 

0.232 

0.046 

249882 

49666 

71 

4013,177 

1932 

0.224 

0.047 

246066 

59194 

58 

4013,277 

537 

0.291 

0.065 

74728 

14232 

64 

4013,317 

1326 

0.203 

0.037 

162763 

32443 

58 

4013,354 

1712 

0.229 

0.045 

201933 

38367 

64 

4013,373 

716 

0.354 

0.145 

136538 

31449 

61 

4013,381 

3079 

0.313 

0.065 

325157 

66566 

71 

tUI  j  ■■^v-l 

6039 

0.30c 

0.068 

36o369 

31069 

64 

4013, ^^64 

2413 

0.240 

0.05B 

306733 

70881 

53 

4013.473 

664 

0.247 

0.072 

70133 

14240 

75 

4013,518 

522 

0.362 

0.115 

107996 

•  28307 

53 

4013,461 

7801 

0.240 

0.078 

260305 

34007 

■^1 

4013,664 

1567 

0.246 

0.075 

193674 

52565 

64 

4013,694 

9314 

0.241 

0.053 

262115 

62856 

71 

4013,723 

3458 

0.230 

0.053 

268779 

62309 

64 

4013,776 

9750 

0.209 

0.040 

403424 

77526 

40 

4013,797 

7024 

0.208 

0.050 

250307 

59593 

64 

*013,30l 

422 

0.251 

0.077 

79331 

15347 

48 

4013,866 

309 

1.157 

0.556 

219930 

826^5 

67 

4013,375 

324 

0.519 

0.247 

122952 

5a24i 

J-.' 

4013,950 

2536 

0.251 

0.045 

350436 

63028 

51 

4013,979 

7943 

0.214 

0.042 

360405 

70664 

46 

4013,1014 

2894 

0.130 

0.043 

234500 

69660 

46 

4013.1018 

590 

0.295 

0.076 

37526 

16507 

7', 

4013,110? 

335 

0.435 

0.134 

206443 

51d43 

^4 

4013,1124 

12137 

0.290 

0.071 

316572 

77457 

71 

4013,1166 

7055 

0.222 

0.063 

320351 

90645 

53 

4013,1137 

982 

l"i  T  =  c. 

0.103 

139420 

31011 

30 

4013,1203 

8126 

0.194 

0.043 

234018 

5239<^ 

64 

4015,232 

434 

0.663 

0.223 

161500 

43620 

53 

4019.178 

19953 

0.312 

0.054 

382727 

66870 

64 

4021,81 

17603 

0.233 

0.0-19 

349159 

71641 

53 

4027,28 

1068 

0.263 

0.195 

145002 

51904 

64 

4029,666 

10375 

0.172 

0.027 

249799 

39740 

53 

4029,397 

364 

0.162 

0,052 

33400 

7374 

64 

4029,105'^ 

9695 

y  «  i  i  --' 

0.046 

307988 

66416 

■-)■-= 

4031.823 


40f 


0.246 


0.102 


66153 


16474 


64 
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Table  L-4  Event  Strenatn  SusiiTsarv 


Nase 

iTape  I,  event  i) 

4033,490 

4033,1043 

4033,1053 

4040,143 
4040,385 
4040,393 
4040.714 
4040,730 
4040,371 
4040,1157 

4047,984 

4049,1064 
4049,1066 

4061,379 

4063, 3o4 

4063,936 

2001 .543 


^ 

Po 

Bioiita 

Fo 

siqiTia 

f 

ffl) 

(Pa) 

(Pa) 

(N) 

!N) 

(Hz) 

4896 

C.232 

0.048 

386B11 

30418 

46 

10359 

0. 190 

0.027 

674418 

94838 

309 

0.365 

0.120 

68367 

19654 

69 

7355 

0.392 

0.055 

1119069 

156887 

27 

4268 

0.463 

0.082 

690510 

121747 

51 

12116 

0.438 

0.119 

848671 

229537 

40 

15100 

0.-47 

O.llS 

75o999 

19^581 

4  J 

13291 

0.491 

0.105 

962200 

204663 

40 

9545 

0.591 

0.093 

1061264 

167076 

43 

10823 

0.457 

0.103 

786054 

135924 

45 

499 

0.759 

0.252 

214017 

54162 

60 

8435 

1.032 

0.135 

3389218 

424279 

i.-J 

B221 

1.281 

0.242 

4938915 

931151 

20 

4157 


0.374 


0.353   1152659    10772/ 


2412 

0.313 

0.056 

545947 

96319 

40 

12033 

0.562 

0.063 

2040664 

245937 

21 

1452 

0.290 

0.027 

505945 

31353 

27 

all') 


U.^36 


0.055 


699233     9440 


200V, 4 jo 
2009,455 
2009.925 

3001,11 
3001,301 

.3047,169 
3047,751 


5274 

0.470 

0.057 

1294184 

1 0  f  3  5  i 

23 

2466 

0.405 

0.06^ 

973961 

156751 

29 

c51 

0.2^9 

0.041 

193006 

62324 

3o 

3922 

0.274 

0.045 

772736 

123159 

27 

3081 

0.214 

0.039 

395029 

71946 

40 

19233 

0.237 

0.037 

640540 

100804 

29 

1040 

0.230 

0.112 

150767 

56491 

64 

151 


Hean  = 


431151 


JEV  = 


Hax  = 

«in  = 


4938915 
33400 
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APPENDIX  D 


Table  D— 1:   Angles,  Ranges  and  Times  for  Re+ractiv« 

Propagation  Paths 

Table  D— 2:   Spreading  Loss  Function,  G(r) 
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Table  D-1  Angles,  Ranges  and  Tiies  for  Refractive  Propagation  Paths 

Theta  :ero      Theta  1  Theta  2  Theta  3  Theta  4     z 

Theta  zero  is  the  surface  launch  angle. 

Theta  1  is  the  angle  at  80  meters. 

Theta  2  is  the  angle  at  the  hydrophone  (93«). 

Theta  3  is  the  angle  at  254  aieters. 

Theta  4  is  the  angle  at  363  meters. 

:  IS  the  aaxiinuB  depth  of  the  ray. 

(See  Figure  3-8) 


93 

93 
93 
94 
96 
97 
98 
99 
100 
102 
103 
104 
105 
107 
108 
109 
111 
112 
113 
115 
116 
118 
119 
121 
122 
124 
125 
127 
128 
130 
131 
133 
135 
136 
138 
140 
141 
143 


0.0657 

0.039 

0.001 

0.0658 

0.039 

0.004 

0.0660 

0.040 

0.006 

0.0670 

0.041 

0.013 

0.0680 

0.043 

0.018 

0.0690 

0.045 

0.021 

0.0700 

0.046 

0.024 

0.0710 

0.048 

0.027 

0.0720 

0.049 

0.030 

0.0730 

0.051 

0.032 

0.0740 

0.052 

0.034 

0.0750 

0.053 

0.036 

0.0760 

0.055 

0.038 

0.0770 

0.056 

0.040 

0.0780 

0.058 

0.042 

0.0790 

0.059 

0.044 

0.0800 

0.060 

0.046 

0.0810 

0.062 

0.047 

0.0820 

0.063 

0.049 

0.0830 

0.064 

0.051 

0.0840 

0.065 

0.052 

0.0850 

0.067 

0.054 

0.0860 

0.068 

0.056 

0.0870 

0.069 

0.057 

0.0880 

0.070 

0.059 

0.0890 

0.072 

0.060 

0.0900 

0.073 

0.062 

0.0910 

0.074 

0.063 

0.0920 

0.075 

0,064 

0.0930 

0.077 

0.066 

0.0940 

0.078 

0.067 

0.0950 

0.079 

0.069 

0.0960 

0.080 

0.070 

0.0970 

0.081 

0.07t 

0.0980 

0.083 

0.073 

0.0990 

0.084 

0.074 

0.1000 

0.085 

0.075 

O.lOiO 

0.036 

0.077 
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Anqles  ^or  Refractive  FTOuaaatlon  Paths 


'leta  zer-j 


Theta  1    Theta  2    Theta  3    Theta  4 


0.1020 
0.1030 
0.1040 
0.1050 
0.1060 
0.1070 

o.ioao 

0. 1090 

O.lilO 
0.1120 
0. 1130 
0.1140 
0.1150 
0.1160 
0.1170 
0.1180 
0.1190 
0. 1200 
0.1210 
0.1220 
0.1230 
0.1240 
0. 1250 
0. 1260 
0.1270 
0.1260 
0. 1290 
0.1300 
0.1310 
0.1320 
0.1330 
0.1340 
0.1350 
0.1360 
0.1370 
0.13S0 
0.1390 
0.1400 
0.1410 

J.i-30 
0.1440 
0.1450 
0, 1460 
0. 1470 


0.087 

0.078 

0.089 

0.079 

0.090 

O.OSl 

0.091 

0.082 

0.092 

0.083 

0.093 

0.085 

0.094 

0.086 

0.095 

0.097 

0.097 

0.08S 

0.098 

0.090 

0.099 

0.091 

0.100 

0.092 

0.101 

0.093 

0.102 

0.094 

0.103 

0.096 

0.104 

0.097 

0.106 

0.098 

0.107 

0.099 

0.103 

0.100 

0.109 

0. 102 

0.110 

0.103 

0.111 

0.104 

0.112 

0.105 

0.113 

0.106 

0.114 

0.108 

0.11s 

0.109 

0.117 

0.110 

O.llS 

0.111 

O.IH 

0.112 

0.120 

0.113 

0.121 

0.115 

0.122 

0.116 

0.123 

0.117 

0.124 

0.118 

0.125 

0.119 

0.127 

0.120 

0.123 

0.121 

0.129 

0.123 

0.130 

0.124 

.■|     <  T  1 

l;    i-'=. 

0.132 

0. 134 
0.135 
0.136 
0.137 


0.12t! 
0.129 
0.130 


145 
146 
148 
150 
152 
154 
155 
157 
159 
161 
163 
165 
167 
169 
171 
173 
175 
177 
179 
iSl 
183 

137 
189 
19: 
194 
1% 
193 
200 
202 
205 
207 
209 
212 
214 
216 
219 

ix.  -■ 


236 
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Theta  :9rc 


Anqies  far  Retractive  PrGDaqation  Patns 
Theta  1  Theta  2  Theta  3  Theta  4 


0.1480 
0.1^90 
0.1500 
0.1510 
0.1520 
0.1521 
0.1522 
0.1523 
0.-524 
0.1525 
0.1526 
0.1527 
0.1523 
0.1529 
0.1530 
0.1540 
0.1550 
0.1560 
0.1570 
0.1530 
0.1590 
0.1600 
0.1610 
0.1620 
0.1630 
0.1640 
0.1650 
0.1660 
0. 1670 
0.1630 
0.1681 
0.1682 
0.1683 
0.1634 
0.1635 
0.1686 
0.1687 
0.1638 
C.ii39 
0.1690 
0.1691 

0.1693 
0. 1694 
0.1695 


0.138 

0.133 

243 

0.139 

0.134 

246 

0.140 

0.135 

248 

0.142 

0.136 

251 

0.143 

0.137 

•254 

0.143 

0.137 

254 

0.143 

0.137 

0.002 

254 

0.143 

0.138 

0.006 

255 

0.143 

0.133 

0.003 

iij  D 

0.143 

0.138 

0.010 

256 

0.143 

0.138 

0.011 

■"/C-7 

0.143 

0.138 

0.013 

257 

0.143 

0.133 

0.014 

253 

0.144 

0.138 

0.015 

258 

0.144 

0.133 

0.016 

259 

0.145 

0.139 

0.024 

265 

0.146 

0.140 

0.029 

272 

0.147 

0.142 

0.034 

273 

0.143 

0.143 

0.039 

284 

0.14= 

0.144 

0.043 

291 

0.150 

0.145 

0.046 

297 

0.151 

0.146 

0.050 

304 

0.152 

0.147 

0.053 

310 

0.153 

0.148 

0.056 

317 

0.154 

0.149 

0.059 

323 

'^' . !  55 

0.150 

0.061 

330 

0.156 

0.151 

0.064 

337 

0.157 

0.153 

0.067 

343 

0.159 

0.154 

0.069 

350 

0.160 

0.155 

0.071 

yc-i 

0. 160 

0.155 

0.072 

358 

0.160 

0.155 

0.072 

358 

0.160 

0.155 

0.072 

359 

0.160 

0.155 

0.072 

360 

0.160 

0.155 

0.073 

360 

0.160 

0.155 

0.073 

361 

0.160 

0.155 

0.073 

362 

0.160 

0.156 

0.073 

0.005 

364 

0.161 

0.15t 

0.074 

0.003 

363 

0.161 

0.156 

0.074 

0.010 

371 

0.161 

0.156 

0.074 

0.011 

374 

O.lti 

0.156 

0,074 

0.013 

378 

0.161 

0.156 

0.074 

0.014 

381 

0.161 

0.156 

0.075 

0.015 

334 

0.161 

0.156 

0.075 

0.016 

358 

0.161 

0. 156 

0.075 

0.017 

391 
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:ero 


Angles  for  Re-fractr-e  Propagatian  Path; 
Theta  1  Theta  2  Thets  3  Theta  4 


0.1697 
0,1698 
0.1699 
0. 1700 
0.1701 
0.1702 
0.1703 
0.1704 
0. 1705 
0. 1706 
0.1707 
0.1703 
0.1709 
0.1710 
0.1711 
0.1712 
0.1713 
0.1714 
0.1715 
0.1716 
0.1717 
0.1718 
0.1719 
0.1720 
0.1730 
0.1740 
0. 1750 
0. 1760 
0.1770 
0.1780 
0.1790 
0.1800 
0.1900 
0. 2000 
0.2100 
0.2200 
0.2300 
0. 2400 
0.2500 
0.2600 
0.2700 
0.28j0 
0.2900 
0.3000 
0.3200 


0.161 

0.157 

0.075 

0.013 

395 

O.lfel 

0.157 

0.076 

0.019 

393 

0.162 

0.157 

0.076 

0.020 

401 

0.162 

0.157 

0.076 

0.021 

405 

0.162 

0.157 

0.076 

0.022 

408 

0.162 

0.157 

0.076 

0.022 

411 

0.162 

0.157 

0.077 

0.023 

415 

0.162 

0.157 

0.077 

0.024 

413 

0.162 

0.157 

0. 077 

0.025 

422 

0.162 

0.158 

0.077 

0.025 

425 

0.162 

0.158 

0.073 

0,026 

428 

0.163 

0.158 

0.078 

0.027 

432 

0.163 

0.153 

0.073 

0,027 

435 

0.163 

0.158 

0.078 

0.028 

439 

0.163 

0.153 

0.073 

0.023 

442 

0.163 

0.158 

0.079 

0.029 

445 

0.163 

0.158 

0.079 

0,030 

449 

0.163 

0.15B 

0.079 

0.030 

452 

0.163 

0.159 

0.079 

0,031 

456 

0. 163 

0.159 

0.030 

0.031 

459 

0.163 

0.159 

0.030 

0.032 

462 

0.164 

0.159 

0.080 

0.032 

46c 

0.164 

0.159 

0.030 

0.033 

469 

0.164 

0.159 

0.030 

0.034 

473 

0.165 

0.160 

0.033 

0.033 

507 

0.166 

0.161 

0.085 

0.043 

542 

0.167 

0.162 

0.037 

0,047 

576 

0.160 

0.163 

0.089 

0,050 

611 

0.169 

0.164 

0,091 

0.054 

647 

0.170 

0.166 

0.093 

0,057 

632 

0.171 

0.167 

0.095 

0.060 

718 

0.172 

0.168 

0.096 

0.063 

753 

0.133 

0.178 

0.193 

0.139 

0.203 

0.200 

0.214 

0.210 

0.224 

0.221 

0.234 

0.231 

0.244 

0.241 

0.255 

0.252 

0.265 

0.262 

0.275 

0.272 

0.235 

0.283 

0.295 

0.293 

0.316 

0.313 

0.336 

0.334 
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"eti  LsrG 


Anglss  tor  Retractive  PropagotiGn  Path?, 
Iheta  1  Theta  2 


0.3600 
0.3800 
0.4000 
0.4500 
0.5000 
0.5500 
0.6000 
0.6500 


0.356  0.354 
0.37?  0.375 
0.397        0.395 


0.447 
0.497 
0.543 
0.593 
0.q4S 


0.446 
0.496 
0.546 
0.597 
0.647 


0.3000 
0.9000 
1 . 0000 
1.1000 
1 . 2000 
1.3000 
1 . 4000 
1 . 5000 
1 . 5708 


0.799 

0.796 

0.399 

0.398 

0.999 

0.999 

1.099 

1.099 

1.199 

1.199 

1.300 

1.299 

1.400 

1.400 

1.500 

1.500 

1.571 

1.571 
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Ran5S3  far  Retractive  ProDaQation  PathB 


Fr.eta  zero 


rl 


n 


r4 


rti 


Rl 


rl  15  the  horizontal  projection  o-f  the  path  froai  the  source  (Onsi  to  8Cni. 

r2  15  the  horizontal  projection  oi   the  path  -froiii  30m  to  the  hydropione  depth. 

r3  is  the  hori:ontai  projection  oi   the  path  troai  the  hydrophone  oepth  to  25451,  or  ue   vt-r 

r4  15  the  horizontal  projection  of  the  path  from  2j4.ti  to  362iti,  or  tne  yerte;:,  if  sooner. 

r5  15  the  horizontal  projection  of  the  path  frois  362a  to  the  vertex. 

Rl  15  the  horizontal  distance  to  the  hydrophone  on  the  downward  path  of  tne  rav. 

R2  15  the  horizontal  distance  to  the  hydrophone  on  the  upi^ard  swing  of  the  ray, 

iSee  Fioure  3-8! 


0.0657 
0.0653 
0.0660 
0.0670 
0.0660 
0.0690 
0.0700 
0.0710 
0.0720 
0.0730 
0.0740 
0.0750 
0.0760 
0.0770 
0.07B0 
0.0790 
0.0300 
0.0310 
0.')320 
0.0330 
0.0840 
0,0050 
0 . 0860 
0.0870 
0.0830 
0.0390 
0.0900 
0.0910 
0,0920 
0.C930 
0, 0940 
0.0950 
0.i9d0 
0.0970 
0.0980 
0 . 0990 
0, 1000 


523.6 

640.6 

663.7 

519.7 

600,9 

666.  J 

512.0 

563, 2 

672.1 

475.1 

476.3 

700.0 

440.4 

429.2 

727.1 

407.3 

395.3 

753.7 

377.0 

369.8 

779.7 

347.3 

348.6 

305.2 

320.1 

330.8 

830.3 

293.3 

315.4 

355.1 

268.7 

302,0 

379.4 

244.7 

290.0 

903.4 

221.0 

279,3 

927,1 

199,9 

269.6 

950,6 

173.9 

260.8 

9". 7 

158.8 

'-7  7 

?96.6 

1  'C  i 
1  j^  ■  -• 

2-^5.3 

1019.3 

120,7 

233,4 

1041.3 

102,8 

232.0 

1C64.1 

035,5 

226.0 

10S6.2 

ObS.3 

220.4 

1108.1 

052.7 

215.2 

1129.3 

037.1 

210.2 

1151.4 

022.0 

205,5 

1172.9 

007.4 

201.1 

1194.2 

993.3 

196.9 

1215.4 

979.6 

193,0 

1236.4 

966.3 

!39,2 

1257.3 

953,4 

135.0 

1270.2 

^40.^ 

132.1 

1298,9 

92S.7 

173.8 

1319.5 

916.9 

175.7 

1340.0 

905.4 

i  i'  2 .  u 

13dO,4 

85-4.2 

169.7 

1380.3 

383.3 

166.9 

1401.0 

872.7 

164.3 

1421.2 

Z  ■'  '  '^ 

i-:1.7 

i4M,3 

2164,2 

2210, 

7^ 

2120.6 

2251 

3 

2075.2 

2293, 

1 

1 

1951.4 

2393 

7 

1869.6 

24b5, 

5 

1303.5 

2519 

4 

1746.8 

2566, 

J 

1696.4 

2609 

6 

1650.9 

2650, 

0 

1609.2 

26BS 

5 

1570.6 

■--!--- 

- 

LI     4.  J 

-! 

1534.7 

27ol 

6 

1501.2 

2^96. 

8 

1469.6 

2831 

4 

143=^.7 

"  0  L  C 

6 

1411,5 

2399 

T 

1334.7 

2932 

3 

1359.2 

29n6 

0 

1334,8 

29QQ 

0 

1311.5 

3031 

3 

1239.2 

3064 

J 

1267.3 

3097 

- 

1247,3 

3129 

7 

1227,6 

3162 

i. 

1208,5 

3194 

6 

1190,2 

3227 

1 

1172,5 

men 

4 

1155,5 

3291 

3 

1139. C 

3324 

i. 

1123.0 

3356 

5 

'  in'i   '■ 

T-.CZ 

- 

i  1 U  ;  .  J 

-t-.'i-u 

"" 

1  ,-••..-.-,  c 

:42l 

-■ 

1078.1^ 

3453 

£ 

1^63.9 

3436 

0 

1050.2 

3518 

i 

t 

1036.9 

3550 

,8 

1024.0 

3583 

- 
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Ranass  tar  Ratrictivs  Propaqstian  Paths 


rl 


r2 


ro* 


r4 


Tj 


o.ioio 

0.1020 
0.1030 
0.1040 
0.1050 
0.1060 
0.1070 
0.1 080 
0..090 
O.ilCO 
0.1110 
0.1120 
0.1130 
■).1I40 
0.1150 
0.1160 
0.1170 
0.1180 
0.1190 
0.12OO 
0.1210 
0.1220 
0.1230 
0. 1240 

A.  -e-ri 


0.1 


2S0 
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Tises  '"or  Refractive  -topaqation  Paths 
tl     t2     t3     t4     t5 


1 1 


tl,  t2.  t3,  t4.  and  t5  correspcnd  to  the  times  required  to  traverse 

the  paths  associated  with  rl,  r2,  r3,  r4  and  r5, 

Tl  arid  12  correspond  to  the  tiiiies  required  to  traverse  the  distances  Rl  and  R2. 

(See  Figure  3-6) 
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Tisiss  tor  Retractive  Propaqstion  Paths 
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Table  D-2  Spreading  Loss  Function,  8(r) 

Theta  zero      Theta  1     z     Rl     R2    6(R1)   6iR2)   10  log   10  log 

G(Rl)   G(R2) 

Theta  zero  is  the  surface  launch  angle. 

Theta  1  is  the  angle  of  the  refractive  path  at  the  hydrophone  depth  (93ffl) 

z  is  the  laxiiua  depth  of  the  refractive  path. 

Rl  is  the  horizontal  range  to  the  hydrophone  intersected  on  the  downward  swing 

R2  is  the  horizontal  range  to  the  hydrophone  intersected  on  the  upward  swing 

6(R)  is  the  spreading  loss  function. 
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